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PEEFACE. 


The  present  publication  is  an  essay  that  was  sent  in  (De- 
cember, 1887)  to  compete  for  the  Smith's  Prizes  at  Cambridge. 

To  the  onlooker  it  is  always  a  mournful  thing  to  see  what 
he  considers  splendid  abilities  or  opportunities  wasted  for  lack  of 
knowledge  of  some  paltry  common-place  truth.  Such  is  in  the 
main  my  feeling  when  considering  the  neglect  of  the  study  of 
Quaternions  by  that  wonderful  corporation  the  University  of  Cam- 
bridge. To  the  alumnus  she  is  apt  to  appear  as  the  leader  in  all 
branches  of  Mathematics.  To  the  outsider  she  appears  rather  as 
the  leader  in  Applied  Mathematics  and  as  a  ready  welcomer  of 
other  branches. 

If  Quaternions  were  simply  a  branch  of  Pure  Mathematics  we 
could  understand  why  the  study  of  them  was  almost  confined 
to  the  University  which  gave  birth  to  them,  but  as  the  truth  is 
quite  otherwise  it  is  hard  to  shew  good  reason  why  they  have  not 
struck  root  also  in  Cambridge.  The  prophet  on  whom  Hamilton's 
mantle  has  fallen  is  more  than  a  mathematician  and  more  than  a 
natural  philosopher — he  is  both,  and  it  is  to  be  noted  also  that  he 
is  a  Cambridge  man.  He  has  preached  in  season  and  out  of 
season  (if  that  were  possible)  that  Quaternions  are  especially  useful 
in  Physical  applications.  Why  then  has  his  Alma  Mater  turned  a 
deaf  ear  ?  I  cannot  believe  that  she  is  in  her  dotage  and  has  lost 
her  hearing.     The  problem  is  beyond  me  and  I  give  it  up. 

But  I  wish  to  add  my  little  efforts  to  Prof.  Tait's  powerful 
advocacy  to  bring  about  another  state  of  affairs.  Cambridge  is  the 
prepared  ground  on  which  if  anywhere  the  study  of  the  Physical 
applications  of  Quaternions  ought  to  flourish. 
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When  I  sent  in  the  essay  I  had  a  faint  misgiving  that  per- 
chance there  was  not  a  single  man  in  Cambridge  who  could 
understand  it  without  much  labour — and  yet  it  is  a  straightfor- 
ward application  of  Hamilton's  principles.  I  cannot  say  what 
transformation  scene  has  taken  place  in  the  five  years  that  have 
elapsed,  but  an  encouraging  fact  is  that  one  professor  at  any  rate 
has  been  imported  from  Dublin. 

There  is  no  lack  in  Cambridge  of  the  cultivation  of  Quater- 
nions as  an  algebra,  but  this  cultivation  is  not  Hamiltonian, 
though  an  evidence  of  the  great  fecundity  of  Hamilton's  work. 
Hamilton  looked  upon  Quaternions  as  a  geometrical  method,  and  it 
is  in  this  respect  that  he  has  as  yet  failed  to  find  worthy  followers 
resident  in  Cambridge.  [The  chapter  contributed  by  Prof.  Cayley 
to  Prof.  Tait's  3rd  ed.  of  'Quaternions'  deals  with  quite  a  different 
subject  from  the  rest  of  the  treatise,  a  subject  that  deserves  a 
distinctive  name,  say,  Cayley  an  Quaternions.] 

I  have  delayed  for  a  considerable  time  the  present  publication 
in  order  at  the  last  if  possible  to  make  it  more  effective.  I  have 
waited  till  I  could  by  a  more  striking  example  than  any  in  the 
essay  shew  the  immense  utility  of  Quaternions  in  the  regions  in 
which  I  believe  them  to  be  especially  powerful.  This  I  believe 
has  been  done  in  the  '  Phil.  Trans.'  1892,  p.  685.  Certainly  on 
two  occasions  copious  extracts  have  been  published,  viz.  in  the 
P.  R.  S.  E.,  1890-1,  p.  98,  and  in  the  '  Phil.  Mag.'  June  1892,  p. 
477,  but  the  reasons  are  simple.  The  first  was  published  after  the 
subject  of  the  '  Phil.  Trans.'  paper  had  been  considered  sufficiently 
to  afford  clear  daylight  ahead  in  that  direction,  and  the  second 
after  that  paper  had  actually  been  despatched  for  publication. 

At  the  time  of  writing  the  essay  I  possessed  little  more  than 
faith  in  the  potentiality  of  Quaternions,  and  I  felt  that  something 
more  than  faith  was  needed  to  convince  scientists.  It  was  thought 
that  rather  than  publish  in  driblets  it  were  better  to  wait  for 
a  more  copious  shower  on  the  principle  that  a  well-directed  heavy 
blow  is  more  effective  than  a  long-continued  series  of  little 
pushes. 

Perhaps  more  harm  has  been  done  than  by  any  other  cause  to 
the  study  of  Quaternions  in  their  Physical  applications  by  a  silly 
superstition   with  which   the   nurses   of  Cambridge   are    wont   to 
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frighten  their  too  timorous  charges.  This  is  the  belief  that  the 
subject  of  Quaternions  is  difficult.  It  is  difficult  in  one  sense  and 
in  no  other,  in  that  sense  in  which  the  subject  of  analytical  conies 
is  difficult  to  the  schoolboy,  in  the  sense  in  which  every  subject 
is  difficult  whose  fundamental  ideas  and  methods  are  different 
from  any  the  student  has  hitherto  been  introduced  to.  The  only 
way  to  convince  the  nurses  that  Quaternions  form  a  healthy  diet 
for  the  young  mathematician  is  to  prove  to  them  that  they  will 
"  pay"  in  the  first  part  of  the  Tripos.  Of  course  this  is  an  impos- 
sible task  while  the  only  questions  set  in  the  Tripos  on  the  subject 
are  in  the  second  part  and  average  one  in  two  years.  [This 
solitary  biennial  question  is  rarely  if  ever  anything  but  an  exercise 
in  algebra.  The  very  form  in  which  candidates  are  invited,  or  at 
any  rate  were  in  my  day,  to  study  Quaternions  is  an  insult  to  the 
memory  of  Hamilton.  The  monstrosity  "  Quaternions  and  other 
non-commutative  algebras"  can  only  be  paralleled  by  "Cartesian 
Geometry  and  other  commutative  algebras."  When  I  was  in  Cam- 
bridge it  was  currently  reported  that  if  an  answer  to  a  Mathe- 
matical Tripos  question  were  couched  in  Hebrew  the  candidate 
would  or  would  not  get  credit  for  the  answer  according  as  one 
or  more  of  the  examiners  did  or  did  not  understand  Hebrew, 
and  that  in  this  respect  Hebrew  or  Quaternions  were  strictly 
analogous.] 

Is  it  hopeless  to  appeal  to  the  charges  ?  I  will  try.  Let  me 
suppose  that  some  budding  Cambridge  Mathematician  has  fol- 
lowed me  so  far.  I  now  address  myself  to  him.  Have  you  ever 
felt  a  joy  in  Mathematics  ?  Probably  you  have,  but  it  was  before 
your  schoolmasters  had  found  you  out  and  resolved  to  fashion  you 
into  an  examinee.  Even  now  you  occasionally  have  feelings  like 
the  dimly  remembered  ones.  Now  and  then  you  forget  that  you 
are  nerving  yourself  for  that  Juggernaut  the  Tripos.  Let  me 
implore  you  as  though  your  soul's  salvation  depended  on  it  to  let 
these  trances  run  their  utmost  course  in  spite  of  solemn  warnings 
from  your  nurse.  You  will  in  time  be  rewarded  by  a  soul-thrilling 
dream  whose  subject  is  the  Universe  and  whose  organ  to  look 
upon  the  Universe  withal  is  the  sense  called  Quaternions.  Steep 
yourself  in  the  delirious  pleasures.  When  you  wake  you  will  have 
forgotten  the  Tripos  and  in  the  fulness  of  time  will  develop  into  a 
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financial  wreck,  but  in  possession  of  the  memory  of  that  heaven- 
sent dream  you  will  be  a  far  happier  and  richer  man  than  the 
millionest  millionaire. 

To  pass  to  earth — from  the  few  papers  I  have  published  it 
will  be  evident  that  the  subject  treated  of  here  is  one  I  have  very 
much  at  heart,  and  I  think  that  the  publication  of  the  essay  is 
likely  to  conduce  to  an  acceptance  of  the  view  that  it  is  now  the 
duty  of  mathematical  physicists  to  study  Quaternions  seriously. 
I  have  been  told  by  more  than  one  of  the  few  who  have  read  some 
of  my  papers  that  they  prove  rather  stiff  reading.  The  reasons 
for  this  are  not  in  the  papers  I  believe  but  in  matters  which  have 
already  been  indicated.  Now  the  present  essay  reproduces  the 
order  in  which  the  subject  was  developed  in  my  own  mind.  The 
less  complete  treatment  of  a  subject,  especially  if  more  diffuse,  is 
often  easier  to  follow  than  the  finished  product.  It  is  therefore 
probable  that  the  present  essay  is  likely  to  prove  more  easy 
reading  than  my  other  papers. 

Moreover  I  wish  it  to  be  studied  by  a  class  of  readers  who  are 
not  in  the  habit  of  consulting  the  proceedings,  &c,  of  learned 
societies.  I  want  the  slaves  of  examination  to  be  arrested  and  to 
read,  for  it  is  apparently  to  the  rising  generation  that  we  must 
look  to  wipe  off  the  blot  from  the  escutcheon  of  Cambridge. 

And  now  as  to  the  essay  itself.  But  one  real  alteration  has 
been  made.  A  passage  has  been  suppressed  in  which  were  made 
uncomplimentary  remarks  concerning  a  certain  author  for  what 
the  writer  regards  as  his  abuse  of  Quaternion  methods.  The 
author  in  question  would  no  doubt  have  been  perfectly  well  able 
to  take  care  of  himself,  so  that  perhaps  there  was  no  very  good 
reason  for  suppressing  the  passage  as  it  still  represents  my  con- 
victions, but  I  did  not  want  a  side  issue  to  be  raised  that  would 
serve  to  distract  attention  from  the  main  one.  To  bring  the 
notation  into  harmony  with  my  later  papers  dv  and  V  which 
occur  in  the  manuscript  have  been  changed  throughout  to  d^i  and 
A  respectively.  To  facilitate  printing  the  solidus  has  been  freely 
introduced  and  the  vinculum  abjured.  Mere  slips  of  the  pen 
have  been  corrected.  A  formal  prefatory  note  required  by  the 
conditions  of  competition  has  been  omitted.  The  Table  of  Con- 
bents  was  not  prefixed  to  the  original  essay.     It  consists  of  little 
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more  than  a  collection  of  the  headings  scattered  through  the 
essay.  Several  notes  have  been  added,  all  indicated  by  square 
brackets  and  the  date  (1892  or  1893).  Otherwise  the  essay 
remains  absolutely  unaltered.  The  name  originally  given  to  the 
essay  is  at  the  head  of  p.  1  below.  The  name  on  the  title-page 
is  adopted  to  prevent  confusion  of  the  essay  with  the  '  Phil.  Mag.' 
paper  referred  to  above.  What  in  the  peculiar  caligraphy  of  the 
manuscript  was  meant  for  the  familiar  fff()d<>  has  been  consistently 
rendered  by  the  printer  as  jff()d%.  As  the  mental  operation  of 
substituting  the  former  for  the  latter  is  not  laborious  I  have  not 
thought  it  necessary  to  make  the  requisite  extensive  alterations  in 
the  proofs. 

I  wish  here  to  express  my  great  indebtedness  to  Prof.  Tait,  not 
only  for  having  through  his  published  works  given  me  such 
knowledge  of  Quaternions  as  I  possess  but  for  giving  me  private 
encouragement  at  a  time  I  sorely  needed  it.  There  was  a  time 
when  I  felt  tempted  to  throw  my  convictions  to  the  winds  and 
follow  the  line  of  least  resistance.  To  break  down  the  solid  and 
well-nigh  universal  scepticism  as  to  the  utility  of  Quaternions  in 
Physics  seemed  too  much  like  casting  one's  pearls — at  least  like 
crying  in  the  wilderness. 

But  though  I  recognise  that  I  am  fighting  under  Prof.  Tait's 
banner,  yet,  as  every  subaltern  could  have  conducted  a  campaign 
better  than  his  general,  so  in  some  details  I  feel  compelled  to 
differ  from  Professor  Tait.  Some  two  or  three  years  ago  he  was 
good  enough  to  read  the  present  essay.  He  somewhat  severely 
criticised  certain  points  but  did  not  convince  me  on  all. 

Among  other  things  he  pointed  out  that  I  sprung  on  the 
unsuspicious  reader  without  due  warning  and  explanation  what 
may  be  considered  as  a  peculiarity  in  symbolisation.  I  take  this 
opportunity  therefore  of  remedying  the  omission.  In  Quaternions 
on  account  of  the  non-commutative  nature  of  multiplication  we 
have  not  the  same  unlimited  choice  of  order  of  the  terms  in  a 
product  as  we  have  in  ordinary  algebra,  and  the  same  is  true  of 
certain  quaternion  operators.  It  is  thus  inconvenient  in  many 
cases  to  use  the  familiar  method  of  indicating  the  connection 
between  an  operator  and  its  operand  by  placing  the  former  im- 
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mediately  before  the  latter.  Another  method  is  adopted.  With 
this  other  method  the  operator  may  be  separated  from  the  operand, 
but  it  seems  that  there  has  been  a  tacit  convention  among  users 
of  this  method  that  the  separated  operator  is  still  to  be  restricted 
to  precedence  of  the  operand.  There  is  of  course  nothing  in  the 
nature  of  things  why  this  should  be  so,  though  its  violation  may 
seem  a  trifle  strange  at  first,  just  as  the  tyro  in  Latin  is  puzzled 
by  the  unexpected  corners  of  a  sentence  in  which  adjectives 
(operators)  and  their  nouns  (operands)  turn  up.  Indeed  a 
Roman  may  be  said  to  have  anticipated  in  every  detail  the 
method  of  indicating  the  connection  now  under  discussion,  for  he 
did  so  by  the  similarity  of  the  suffixes  of  his  operators  and 
operands.  In  this  essay  his  example  is  followed  and  therefore  no 
restrictions  except  such  as  result  from  the  genius  of  the  language 
(the  laws  of  Quaternions)  are  placed  on  the  relative  positions  in  a 
product  of  operators  and  operands.  With  this  warning  the  reader 
ought  to  find  no  difficulty. 

One  of  Prof.  Tait's  criticisms  already  alluded  to  appears  in  the 
third  edition  of  his  '  Quaternions.'  The  process  held  up  in  §  500 
of  this  edition  as  an  example  of  "  how  not  to  do  it "  is  contained 
in  §  6  below  and  was  first  given  in  the  '  Mess,  of  Math.,'  1884. 
He  implies  that  the  process  is  a  "  most  intensely  artificial  applica- 
tion of"  Quaternions.  If  this  were  true  I  should  consider  it  a 
perfectly  legitimate  criticism,  but  I  hold  that  it  is  the  exact 
reverse  of  the  truth.  In  the  course  of  Physical  investigations 
certain  volume  integrals  are  found  to  be  capable  of,  or  by  general 
considerations  are  obviously  capable  of  transformation  into  surface 
integrals.  We  are  led  to  seek  for  the  correct  expression  in  the 
latter  form.  Starting  from  this  we  can  by  a  long,  and  in  my 
opinion,  tedious  process  arrive  at  the  most  general  type  of  volume 
integral  which  is  capable  of  transformation  into  a  surface  integral. 
[I  may  remark  in  passing  that  Prof.  Tait  did  not  however  arrive 
at  quite  the  most  general  type.]  Does  it  follow  that  this  is  the 
most  natural  course  of  procedure  ?  Certainly  not,  as  I  think.  It 
would  be  the  most  natural  course  for  the  empiricist,  but  not  for 
the  scientist.  When  he  has  been  introduced  to  one  or  two  volume 
integrals  capable  of  the  transformation  the  natural  course  of  the 
mathematician  is  to  ask  himself  what  is  the  most  general  volume 
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integral  of  the  kind.  By  quite  elementary  considerations  he  sees 
that  while  only  such  volume  integrals  as  satisfy  certain  conditions 
are  transformable  into  surface  integrals,  yet  any  surface  integral 
which  is  continuous  and  applies  to  the  complete  boundary  of  any 
finite  volume  can  be  expressed  as  a  volume  integral  throughout 
that  volume.  He  is  thus  led  to  start  from  the  surface  integral 
and  deduces  by  the  briefest  of  processes  the  most  general  volume 
integral  of  the  type  required.  Needless  to  say,  when  giving  his 
demonstration  he  does  not  bare  his  soul  in  this  way.  He  thinks 
rightly  that  any  mathematician  can  at  once  divine  the  exact  road 
he  has  followed.     Where  is  the  artificiality  ? 

Let  me  in  conclusion  say  that  even  now  I  scarcely  dare  state 
what  I  believe  to  be  the  proper  place  of  Quaternions  in  a 
Physical  education,  for  fear  my  statements  be  regarded  as  the 
uninspired  babblings  of  a  misdirected  enthusiast,  but  I  cannot 
refrain  from  saying  that  I  look  forward  to  the  time  when 
Quaternions  will  appear  in  every  Physical  text-book  that  assumes 
the  knowledge  of  (say)  elementary  plane  trigonometry. 

I  am  much  indebted  to  Mr  G.  H.  A.  Wilson  of  Clare  College, 
Cambridge,  for  helping  me  in  the  revision  of  the  proofs,  and  take 
this  opportunity  of  thanking  him  for  the  time  and  trouble  he  has 
devoted  to  the  work. 

ALEX.  McAULAY. 


University  of  Tasmania, 
Hobart. 

March  26,  1893. 


CONTENTS. 


Section  I.    Introduction. 


PAGE 


General  remarks  on  the  place  of  Quaternions  in  Physics        ...  1 

Cartesian  form  of  some  of  the  results  to  follow 6 


Section  II.     Quaternion  Theorems. 

1.  Definitions 12 

2.  Properties  of  f ....  15 

5.  Fundamental  property  of  (J 18 

6.  Theorems  in  Integration 19 

9.     Potentials 22 


Section  III.     Elastic  Solids. 

11.  Brief  recapitulation  of  previous  work  in  this  branch       ....  25 

12.  Strain,  stress-force,  stress-couple 26 

14.  Stress  in  terms  of  strain 27 

16.  The  equations  of  equilibrium 32 

16a.  Variation  of  temperature         .......  .35 

17.  Small  strains 37 

20.  Isotropic  bodies 40 

22.  Particular  integral  of  the  equation  of  equilibrium 42 

24.  Orthogonal  coordinates 44 

27.  Saint-Venant's  torsion  problem 46 

29.  Wires 49 


Section  IV.     Electricity  and  Magnetism. 

34.  Electrostatics — general  problem 55 

41.  The  force  in  particular  cases 63 

43.  Nature  of  the  stress 65 

46.  Magnetism — magnetic  potential,  force,  induction 67 

49.  Magnetic  solenoids  and  shells 69 

54.  Electro-magnetism — general  theory       .......  71 

60.  Electro-magnetic  stress  ..........  74 


XIV 


CONTENTS. 


Section  V.     Hydrodynamics. 

PAGE 

61.  Preliminary 76 

62.  Notation 76 

63.  Euler's  equations 77 

68.  The  Lagrangian  equations       .........  80 

69.  Cauchy's  integrals  of  these  equations       .......  81 

71.     Flow,  circulation,  vortex-motion 82 

74.     Irrotational  motion          ..........  84 

76.     Motion  of  a  solid  through  a  liquid 85 

79.     The  velocity  in  terms  of  the  convergences  and  spins        ....  88 

83.     Viscosity 91 


Section  VI.    The  Vortex-atom  Theory. 


85.  Preliminary     ............  94 

86.  Statement  of  Sir  Win.  Thomson's  and  Prof.  Hicks's  theories  ...  94 

87.  General  considerations  concerning  these  theories    .....  95 

88.  Description  of  the  method  here  adopted  .......  96 

89.  Acceleration  in  terms  of  the  convergences,  their  time-Muxes  and  the 

spins          ............  97 

91.     Sir  Wm.  Thomson's  theory 98 

93.  Prof.  Hicks's  theory 100 

94.  Consideration  of  all  the  terms  except  -  v  {ff")l^       .....  101 
1)6.     Consideration  of  the  term  -  v  •  (c"')/2 103 


QUATERNIONS   AS   A   PRACTICAL   INSTRU- 
MENT  OF  PHYSICAL  RESEARCH. 

SECTION   I. 

Introduction. 

It  is  a  curious  phenomenon  in  the  History  of  Mathematics 
that  the  greatest  work  of  the  greatest  Mathematician  of  the 
century  which  prides  itself  upon  being  the  most  enlightened  the 
world  has  yet  seen,  has  suffered  the  most  chilling  neglect. 

The  cause  of  this  is  not  at  first  sight  obvious.  We  have  here 
little  to  do  with  the  benefit  provided  by  Quaternions  to  Pure 
Mathematics.  The  reason  for  the  neglect  here  may  be  that 
Hamilton  himself  has  developed  the  Science  to  such  an  extent  as 
to  make  successors  an  impossibility.  One  cannot  however  resist 
a  strong  suspicion  that  were  the  subject  even  studied  we  should 
hear  more  from  Pure  Mathematicians,  of  Hamilton's  valuable 
results.  This  reason  at  any  rate  cannot  be  assigned  for  the 
neglect  of  the  Physical  side  of  Quaternions.  Hamilton  has  done 
but  little  in  this  field,  and  yet  when  we  ask  what  Mathematical 
Physicists  have  been  tempted  by  the  bait  to  win  easy  laurels  (to 
put  the  incentive  on  no  higher  grounds),  the  answer  must  be 
scarcely  one.  Prof.  Tait  is  the  grand  exception  to  this.  But  well- 
known  Physicist  though  he  be,  his  fellow- workers  for  the  most 
part  render  themselves  incapable  of  appreciating  his  valuable 
services  by  studying  the  subject  if  at  all  only  as  dilettanti.  The 
number  who  read  a  small  amount  in  Quaternions  is  by  no  means 
small,  but  those  who  get  further  than  what  is  recommended  by 
Maxwell  as  imperatively  necessary  are  but  a  small  percentage  of 
the  whole. 

I  cannot  help  thinking  that  this  state  of  affairs  is  owing  chiefly 
to  a  prejudice.  This  prejudice  is  well  seen  in  Maxwell's  well- 
known  statement — "  I  am  convinced  that  the  introduction  of  the 
M.  1 
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ideas,  as  distinguished  from  the  operations  and  methods  of  Quater- 
nions, will  be  of  great  use  to  us  in  all  parts  of  our  subject."* 
Now  what  I  hold  and  what  the  main  object  of  this  essay  is  to 
prove  is  that  the  "operations  and  methods"  of  Quaternions  are  as 
much  better  qualified  to  deal  with  Physics  than  the  ordinary  ones 
as  are  the  "  ideas". 

But,  what  has  produced  this  notion,  that  the  subject  of 
Quaternions  is  only  a  pretty  toy  that  has  nothing  to  do  with  the 
serious  work  of  practical  Physics?  It  must  be  the  fact  that  it  has 
hitherto  produced  few  results  that  appeal  strongly  to  Physicists. 
This  I  acknowledge,  but  that  the  deduction  is  correct  I  strongly 
disbelieve.  As  well  might  an  instrument  of  which  nobody  has 
attempted  to  master  the  principles  be  blamed  for  not  being  of 
much  use.  Workers  naturally  find  themselves  while  still  inex- 
perienced in  the  use  of  Quaternions  incapable  of  clearly  thinking 
through  them  and  of  making  them  do  the  work  of  Cartesian  Geo- 
metry, and  they  conclude  that  Quaternions  do  not  provide  suitable 
treatment  for  what  they  have  in  hand.  The  fact  is  that  the  sub- 
ject requires  a  slight  development  in  order  readily  to  apply  to  the 
practical  consideration  of  most  physical  subjects.  The  first  steps 
of  this,  which  consist  chiefly  in  the  invention  of  new  symbols  of. 
operation  and  a  slight  examination  of  their  chief  properties,  I 
have  endeavoured  to  give  in  the  following  pages. 

I  may  now  state  what  I  hold  to  be  the  mission  of  Quaternions 
to  Physics.  I  believe  that  Physics  would  advance  with  both 
more  rapid  and  surer  strides  were  Quaternions  introduced  to 
serious  study  to  the  almost  total  exclusion  of  Cartesian  Geo- 
metry, except  in  an  insignificant  way  as  a  particular  case  of 
the  former.  All  the  geometrical  processes  occurring  in  Physical 
theories  and  general  Physical  problems  are  much  more  grace- 
ful in  their  Quaternion  than  in  their  Cartesian  garb.  To  illus- 
trate what  is  here  meant  by  "theory"  and  "general  problem" 
let  us  take  the  case  of  Elasticity  treated  below.  That  by  the 
methods  advocated  not  only  are  the  already  well-known  results  of 
the  general  theory  of  Elasticity  better  proved,  but  more  general 
results  are  obtained,  will  I  think  be  acknowledged  after  a  perusal 
of  §§  12  to  21  below.  That  Quaternions  are  superior  to  Cartesian 
Geometry  in  considering  the  general  problems  of  (1)  an  infinite 
isotropic  solid,  (2)  the  torsion  and  bending  of  prisms  and  cylinders, 
(3)  the  general  theory  of  wires,  I  have  endeavoured  to  shew  in 

*  Elect,  and  Mag.  Vol.  i.  §  10. 


INTRODUCTION.  3 

§§  22 — 33.  But  for  particular  problems  such  as  the  torsion 
problem  for  a  cylinder  of  given  shape,  we  require  of  course  the 
various  theories  specially  constructed  for  the  solution  of  particular 
problems  such  as  Fourier's  theories,  complex  variables,  spherical 
harmonics,  &c.  It  will  thus  be  seen  that  I  do  not  propose  to 
banish  these  theories  but  merely  Cartesian  Geometry. 

So  mistaken  are  the  common  notions  concerning  the  preten- 
sions of  advocates  of  Quaternions  that  I  was  asked  by  one  well- 
known  Mathematician  whether  Quaternions  furnished  methods  for 
the  solution  of  differential  equations,  as  he  asserted  that  this  was 
all  that  remained  for  Mathematics  in  the  domain  of  Physics ! 
Quaternions  can  no  more  solve  differential  equations  than  Carte- 
sian Geometry,  but  the  solution  of  such  equations  can  be  performed 
as  readily,  in  fact  generally  more  so,  in  the  Quaternion  shape  as  in 
the  Cartesian.  But  that  the  sole  work  of  Physical  Mathematics 
to-day  is  the  solution  of  differential  equations  I  beg  to  question. 
There  are  many  and  important  Physical  questions  and  extensions 
of  Physical  theories  that  have  little  or  nothing  to  do  with  such 
solutions.  As  witness  I  may  call  attention  to  the  new  Physical 
work  which  occurs  below. 

If  only  on  account  of  the  extreme  simplicity  of  Quaternion 
notation,  large  advances  in  the  parts  of  Physics  now  indicated,  are 
to  be  expected.  Expressions  which  are  far  too  cumbrous  to  be  of 
much  use  in  the  Cartesian  shape  become  so  simple  when  translated 
into  Quaternions,  that  they  admit  of  easy  interpretation  and,  what 
is  perhaps  of  more  importance,  of  easy  manipulation.  Compare  for 
instance  the  particular  case  of  equation  (15  in)  §  16  below  when 
F  =  0  with  the  same  thing  as  considered  in  Thomson  and  Tait's 
Nat.  Phil.,  App.  C.     The  Quaternion  equation  is 

pjSVrtQwb.  =  0. 

The  Cartesian  exact  equivalent  consists  of  Thomson  and  Tait's 
equations  (7),  viz. 

d    i^dw  /da        \      dw  da      dw  da) 
dx  {    dA  \dx        J      db  dz      dc  dy) 

d    [    dw  da     div  da     dw  (da        \) 
dy  I    dB  dy      da  dz      dc  \dx        J  \ 

div  da     dw  da     dw  (da        \) 
dC  dz      da  dy      db  \dx        J) 


and  two  similar  equations. 
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Many  of  the  equations  indeed  in  the  part  of  the  essay  where 
this  occurs,  although  quite  simple  enough  to  be  thoroughly  useful 
in  their  present  form,  lead  to  much  more  complicated  equations 
than  those  just  given  when  translated  into  Cartesian  notation. 

It  will  thus  be  seen  that  there  are  two  statements  to  make 
good : — (1)  that  Quaternions  are  in  such  a  stage  of  development 
as  already  to  justify  the  practically  complete  banishment  of  Carte- 
sian Geometry  from  Physical  questions  of  a  general  nature,  and 
(2)  that  Quaternions  will  in  Physics  produce  many  new  results  that 
cannot  be  produced  by  the  rival  and  older  theory. 

To  establish  completely  the  first  of  these  propositions  it  would 
be  necessary  to  go  over  all  the  ground  covered  by  Mathematical 
Physical  Theories,  by  means  of  our  present  subject,  and  compare 
the  proofs  with  the  ordinary  ones.  This  of  course  is  impossible  in 
an  essay.  It  would  require  a  treatise  of  no  small  dimensions. 
But  the  principle  can  be  followed  to  a  small  extent.  I  have  there- 
fore taken  three  typical  theories  and  applied  Quateimions  to  most 
of  the  general  propositions  in  them.  The  three  subjects  are  those 
of  Elastic  Solids,  with  the  thermodynamic  considerations  necessary, 
Electricity  and  Magnetism,  and  Hydrodynamics.  It  is  impossible 
without  greatly  exceeding  due  limits  of  space  to  consider  in  addi- 
tion, Conduction  of  Heat,  Acoustics,  Physical  Optics,  and  the 
Kinetic  Theory  of  Gases.  With  the  exception  of  the  first  of  these 
subjects  I  do  not  profess  even  to  have  attempted  hitherto  the 
desired  applications,  but  one  would  seem  almost  justified  in  argu- 
ing that,  since  Quaternions  have  been  found  so  applicable  to  the 
subjects  considered,  they  are  very  likely  to  prove  useful  to  about 
the  same  extent  in  similar  theories.  Again,  only  in  one  of  the 
subjects  chosen,  viz.,  Hydrodynamics,  have  I  given  the  whole  of 
the  general  theory  which  usually  appears  in  text-books.  For 
instance,  in  Electricity  and  Magnetism  I  have  not  considered  Elec- 
tric Conduction  in  three  dimensions  which,  as  Maxwell  remarks, 
lends  itself  very  readily  to  Quaternion  treatment,  nor  Magnetic 
Induction,  nor  the  Electro-Magnetic  Theory  of  Light.  Again,  I 
have  left  out  the  consideration  of  Poynting's  theories  of  Elec- 
tricity which  are  very  beautifully  treated  by  Quaternions,  and  I 
felt  much  tempted  to  introduce  some  considerations  in  connection 
with  the  Molecular  Current  theory  of  Magnetism.  With  similar 
reluctance  I  have  been  compelled  to  omit  many  applications  in 
the  Theory  of  Elastic  Solids,  but  the  already  too  large  size  of  the 
essay  admitted  of  no  other  course.     Notwithstanding  these  omis- 
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sions,  I  think  that  what  1  have  done  in  this  part  will  go  far  to 
bear  out  the  truth  of  the  first  proposition  I  have  stated  above. 

But  it  is  the  second  that  I  would  especially  lay  stress  upon. 
In  the  first  it  is  merely  stated  that  Cartesian  Geometry  is  an  anti- 
quated machine  that  ought  to  be  thrown  aside  to  make  room  for 
modern  improvements.  But  the  second  asserts  that  the  improved 
machinery  will  not  only  do  the  work  of  the  old  better,  but  will 
also  do  much  work  that  the  old  is  quite  incapable  of  doing  at  all. 
Should  this  be  satisfactorily  established  and  should  Physicists  in 
that  case  still  refuse  to  have  anything  to  do  with  Quaternions, 
they  would  place  themselves  in  the  position  of  the  traditional 
workmen  who  so  strongly  objected  to  the  introduction  of  ma- 
chinery to  supplant  manual  labour. 

But  in  a  few  months  and  synchronously  with  the  work  I  have 
already  described,  to  arrive  at  a  large  number  of  new  results  is  too 
much  to  expect  even  from  such  a  subject  as  that  now  under  dis- 
cussion.    There  are  however  some  few  such  results  to  shew.     I 
have  endeavoured  to  advance  each  of  the  theories  chosen  in  at 
least  one  direction.     In  the  subject  of  Elastic  Solids  I  have  ex- 
pressed the  stress  in  terms  of  the  strain,  in  the  most  general  case, 
i.e.  where  the  strain  is  not  small,  where  the  ordinary  assumption 
of  no  stress-couple  is  not  made  and  where  no  assumption  is  made 
as  to  homogeneity,  isotropy,  &c.     I  have  also  obtained  the  equations 
of  motion  when  there  is  given  an  external  force  and  couple  per 
unit  volume  of  the  unstrained  solid.     These  two  problems,  as  will 
be  seen,  are  by  no  means  identical.     In  Electrostatics  I  have  con- 
sidered the  most  general  mechanical  results  flowing  from  Maxwell's 
theory,  and  their   explanation  by  stress  in  the  dielectric.     These 
results  are  not  known,  as  might  be  inferred   from   this  mode  of 
statement,  for  to  solve  the  problem  we  require  to  know  forty-two 
independent  constants  to  express  the  properties  of  the  dielectric  at 
a  given  state  of  strain  at  each  point.    These  are  the  six  coefficients 
of  specific  inductive  capacity  and  their  thirty-six  differential  coeffi- 
cients with  regard  to  the  six  coordinates  of  pure  strain.     But,  as 
far  as  I  am  aware,  only  such  particular  cases  of  this  have  already 
been  considered  as  make  the  forty-two  constants  reduce  at  most  to 
three.     In  Hydrodynamics  I  have  endeavoured  to  deduce  certain 
general  phenomena  which   would   be  exhibited  by  vortex-atoms 
acting  upon  one  another.     This  has  been  done  by  examination  of 
an  equation  which  has  not,  I  believe,  been  hitherto  given.     The 
result  of  this  part  of  the  essay  is  to  lead  to  a  presumption  against 
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Sir  William  Thomson's  Vortex-Atom  Theory  and  in  favour  of 
Hicks's. 

As  one  of  the  objects  of  this  introduction  is  to  give  a  bird's-eye 
view  of  the  merits  of  Quaternions  as  opposed  to  Cartesian  Geome- 
try, it  will  not  be  out  of  place  to  give  side  by  side  the  Quaternion 
and  the  Cartesian  forms  of  most  of  the  new  results  I  have  been 
speaking  about.  It  must  be  premised,  as  already  hinted,  that  the 
usefulness  of  these  results  must  be  judged  not  by  the  Cartesian 
but  by  the  Quaternion  form. 


Elasticity. 

Let  the  point  (x,  y,  z)  of  an  elastic  solid  be  displaced  to 
(x,  y',  z).  The  strain  at  any  point  that  is  caused  may  be  supposed 
due  to  a  pure  strain  followed  by  a  rotation.  In  Section  III.  below, 
this  pure  strain  is  called  -v/r.  Let  its  coordinates  be  e,  /  g, 
a/2,  b/2,  c/2;  i.e.  if  the  vector  (£  V,  f)  becomes  (£',  77',  £')  by  means 
of  the  pure  strain,  then 

&c,  &c. 

Thus  when  the  strain  is  small  e,  f,  g  reduce  to  Thomson  and  Tait's 
\jr6)\  +f}  1  +  g  and  a,  b,  c  are  the  same  both  in  their  case  and 
the  present  one.  Now  let  the  coordinates  of  M*,  §  16  below,  be 
E,  F,  G,  A/2,  B/2,  C/2.     Equation  (15),  §  16  below,  viz. 

^ives  in  our  present  notation 

E  =  e2  +  c2/4  +  i2/4  =  {dx'jdxf  +  (dy'/dxf  +  (dz'/dxf, 

&c,     &c. 

A  =  a  (/+  g)  +  bc/2  =  2  {(dx'/dy)  (dx'/dz)  +  (dy'/dy)  (dy'/dz) 

+  (dz'jdy)  (dz'/dz)}, 
&e,     &c. 

which   shew  that  the  present  E,  F,   G,  A/ 2,  B/2,   C/2  are  the 
A,  B,  C,  a,  b,  c  of  Thomson  and  Tait's  Nat.  Phil,  App.  C. 

*  This  result  is  one  of  Tait's  (Quaternions  §  365  where  he  has  <p'cp=m2).    It  is 
given  lure  for  completeness. 
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Let  us  put 

\x  y  zl 

J  [       )   =Jn,  &c,  &c, 

\yz  ) 

J\yz)   =  J™'  J\zx)=  J*v  &C  '  &c''  &c"  &c< 

I  have  shewn  in  §  14  below  that  the  stress-couple  is  quite 
independent  of  the  strain.  Thus  we  may  consider  the  stress  to 
consist  of  two  parts — an  ordinary  stress  PQRSTU  as  in  Thomson 
and  Tait's  Nat.  Phil,  and  a  stress  which  causes  a  couple  per  unit 
volume  L'M'N'.  The  former  onry  of  these  will  depend  on  strain. 
The  result  of  the  two  will  be  to  cause  a  force  (as  indeed  can  be 
seen  from  the  expressions  in  §  13  below)  per  unit  area  on  the  x- 
interface  P,  U+N'j%  T-Mj%  and  so  for  the  other  interfaces.  If 
L,  M,  N  be  the  external  couple  per  unit  volume  of  the  unstrained 
solid  we  shall  have 

L'  =  -  LI  J,  M'  =  -  MJ  J,  N'  =  -  N/J, 

for  the  external  couple  and  the  stress  couple  are  always  equal  and 
opposite.     Thus  the  force  on  the  ^-interface  becomes 

P,  U-N/2J,  T+M/2J 

and  similarly  for  the  other  interfaces. 

To  express  the  part  of  the  stress  (P  &c.)  which  depends  on 
the  strain  in  terms  of  that  strain,  consider  w  the  potential  energy 
per  unit  volume  of  the  unstrained  solid  as  a  function  of  E  &c. 
In  the  general  thermodynamic  case  iv  may  be  defined  by  saying 
that  w  x  (the  element  of  volume)  =  (the  intrinsic  energy  of  the 
element)  —  (the  entropy  of  the  element  x  its  absolute  temperature 
x  Joule's  coefficient).  Of  course  10  may  be,  and  indeed  is  in 
§§  14,  15  below,  regarded  as  a  function  of  e  &c. 

The  equation  for  stress  is  (15  b)  §  10  below,  viz., 
JipG)  =  2x*<Jivx'o>  =  2p1/Sf/0s'a»SV1*QwVa. 

The  second  of  the  expressions  is  in  terms  of  the  strain  and  the 
third  in  terms  of  the  displacement  and  its  derivatives.  In  our 
present  notation  this  last  is 
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JP  _  fdx'V dw      fdx'V dw      (dxV dw 
T  ~  \fa)  dE+  \dy)  dF+  [dz  )  dG 

dx'  dx'  dw      9  dx'  dx  dw         dx'  dx  dw 
"  dy  dz  dA         dz  dx  dB         dx  dy  dG ' 
&c,  &c. 

JS     dy'  dz'  dw      dy'  dz  dw     dy'  dz  dw 
2    ~  dx  dx  dE     dy  dy  dF     dz  dz  dG 

(dy'  dz'      dy'  dz'\dw      (dy'  dz'     dy'  dz'\  dw 
\dy  dz       dz  dy)  dA      \  dz  dx      dx  dz  J  dB 

/dy'  dz'     dy'  dz'\  dw 
\dx  dy      dy  dx)  dC 

&c,  &c. 

In  §  14  I  also  obtain  this  part  of  the  stress  explicitly  in  terms 
of  e,  f,  g,  a,  b,  c,  of  w  as  a  function  of  these  quantities  and  of  the 
axis  and  amount  of  rotation.  But  these  results  are  so  very  com- 
plicated in  their  Cartesian  shape  that  it  is  quite  useless  to  give 
them. 

To  put  down  the  equations  of  motion  let  Xz,  Yx,  Zx  be  the 
force  due  to  stress  on  what  before  strain  was  unit  area  per- 
pendicular to  the  axis  of  x.  Similarly  for  Xy,  &c.  Next  suppose 
that  X,  Y,  Z  is  the  external  force  per  unit  volume  of  the  un- 
strained solid  and  let  D  be  the  original  density  of  the  solid.  Then 
the  equation  of  motion  (15  n)  §  16  a  below,  viz. 

Dp'  =  F  +  tA, 
gives  in  our  present  notation 

X  +  dXJdx  +  dXJdy  +  dXJdz  =  isb'D, 
&c,  &c. 

It  remains  to  express  Xx  &c.  in  terms  of  the  displacement  and 
LMN.     This  is  done  in  equation  (15  I)  §  10  below,  viz. 

™  =  -  2/>1'£V1*ciwa>  +  3FMF/Jl>2'^V1V2/2,SV1V.2V^p1>2>;.* 

In  our  present  notation  this  consists  of  the  following  nine 
equations : 

*  The  second  term  on  the  right  contains  in  full  the  nine  terms  corresponding  to 
(JV,N  -JV]M)I'2.J.  Quaternion  notation  is  therefore  here,  as  in  nearly  all  cases 
which  occur  in  Physics,  considerably  more  compact  even  than  the  notations  of 
determinants  or  Jaeobiuns. 


dE  dx      dC  dy      dB  dz  J  2  J" 

Y      9  (dw  dtf     dm  dy^     dw  d_i/\      J1SL  -  JnN 
x        \dE  dx  "*"  dC  dy^dBdz)^'       2J 

„         (dw  dz       dw  dz       dw  dz'\      JnM  —  JV,L 

^x  —  <M  jp  J"    Tj7iT:+TnT-      + 
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(dw  dx      dw  dx      dw  dx'\     JnN  —  JWM 


,rf^  cfe    '  dC  dij      dBdzJ^         2J 
and  six  similar  equations. 

We  thus  see  that  in  the  case  where  LMN  are  zero,  our  present 
Xx,  Xv,  X3  are  the  PQR  of  Thomson  and  Tait's  Nat  Phil.  App.  C 
(d),  and  therefore  equations  (7)  of  that  article  agree  with  our 
equations  of  motion  when  we  put  both  the  external  force  and  the 
acceleration  zero. 

These  are  some  of  the  new  results  in  Elasticity,  but,  as  I  have 
hinted,  there  are  others  in  §§  14,  15  which  it  would  be  waste  of 
time  to  give  in  their  Cartesian  form. 


Electricity. 

In  Section  IV.  below  I  have  considered,  as  already  stated,  the 
most  general  mechanical  results  flowing  from  Maxwell's  theory  of 
Electrostatics.  I  have  shewn  that  here,  as  in  the  particular  cases 
considered  by  others,  the  forces,  whether  per  unit  volume  or  per 
unit  surface,  can  be  explained  by  a  stress  in  the  dielectric.  It  is 
easiest  to  describe  these  forces  by  means  of  the  stress. 

Let  the  coordinates  of  the  stress  be  PQRSTU.  Then  FXF2FZ 
the  mechanical  force,  due  to  the  field  per  unit  volume,  exerted 
upon  the  dielectric  where  there  is  no  discontinuity  in  the  stress,  is 
given  by 

Ft  =  dP/dx  +  d  U/dy  +  dT/dz,  &c,  &c. 

and  (I,  m,  n)  being   the  direction  cosines    of  the  normal  to  any 
surface,  pointing  away  from  the  region  considered 

F;  =-[lP  +  m  U  +  nT]a  -  [  I,  &c,  &c, 

where  a,  b  indicate  the  two  sides  of  the  surface  and  F(,  F2',  F3'  is 
the  force  due  to  the  field  per  unit  surface. 

It  remains  to  find  P  &c.  Let  X,  Y,  Z  be  the  electro-motive 
force,  a,  ft,  7  the  displacement,  w  the  potential  energy  per  unit 
volume  and  Kxx,  Kuy,  K„,  Kyz,  Kzz,  Kxy  the  coefficients  of  specific 
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inductive  capacity.  Let  1  +e,  1  +/,  1+g,  a/2,  b/2,  c/2  denote  the 
pure  part  of  the  strain  of  the  medium.  The  K's  will  then  be 
functions  of  e  &c.  and  we  must  suppose  these  functions  known,  or 
at  any  rate  we  must  assume  the  knowledge  of  both  the  values  of 
the  K's  and  their  differential  coefficients  at  the  particular  state  of 
strain  in  which  the  medium  is  when  under  consideration.  The 
relations  between  the  above  quantities  are 

4ma  =  KxxX  +  KscyY  +  KexZ,  &c,  &c. 

w  =  (Xa  +  Y/3  +  Zry)/2 
=  (KXXX*  +  KyJ>  +  KJF  +  2KnTZ  +  2KZXZX 

+  2KxuXY)/^. 

It  is  the  second  of  these  expressions  for  w  which  is  assumed 
below,  and  the  differentiations  of  course  refer  only  to  the  K's. 
The  equation  expressing  P  &c.  in  terms  of  the  field  is  (21)  §  40 
below,  viz. 

0o)  =  —  ^VDcoE  —  ^qww, 

which  in  our  present  notation  gives  the  following  six  equations 

p= _  i  (  _  aX  +  /3F+  yZ)  -  dw/de,  &c,  &c, 
S  =  \  (/3Z  +  yY)  —  dw/da,  &c,  &c. 

I  have  shewn  in  §§  41 — 45  below  that  these  results  agree  with 
particular  results  obtained  by  others. 


Hydrodynamics. 

The  new  work  in  this  subject  is  given  in  Section  VI. — "  The 
Vortex-Atom  Theory."  It  is  quite  unnecessary  to  translate  the 
various  expressions  there  used  into  the  Cartesian  form.  I  give 
here  only  the  principal  equation  in  its  two  chief  forms,  equation 
(9)  §  89  and  equation  (11)  §  90,  viz. 

P  +  v  -  a2/ 2  +  ^Try'JjfiSarVa  +  udmjdt)  d%  =  H, 

P  +  V-  ol/2  +  (47r)-yj"/  [dsSVu  (  Var  -  ma)  +  ud  (md%)/dt]  =  H. 

In  Cartesian  notation  these  are 
fdp/p+V  +  q*/2 

-  (4tt)-1///{2  [{x'  -  x)  (wV  -  0  + . . .  +. .  .]/r:!  +  (dc/dt)/r}dx'dy'dz'  =  H. 

fdp/p+V+q*/2 

-  (**Tllff  {(x  -  x)  [2  (wV  -  v0  -  cu]  +...  +  .  ..}lr3 .  dx'dy'dz 

-  (^n//  [d  {cdxdy'dz')/dt}/r  =  H. 
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The  fluid  here  considered  is  one  whose  motion  is  continuous 
from  point  to  point  and  which  extends  to  infinity.  The  volume 
integral  extends  throughout  space.  The  notation  is  as  usual. 
It  is  only  necessary  to  say  that  H  is  a  function  of  the  time  only, 
r  is  the  distance  between  the  points  x',  y',  z'  and  x,  y,  z\ 

c  =  du/dx  +  dv/dy  +  dwjdz ; 

djdt  is  put  for  differentiation  which  follows  a  particle  of  the  fluid, 
and  d/dt  for  that  which  refers  to  a  fixed  point. 

The  explanation  of  the  unusual  length  of  this  essay,  which  I 
feel  is  called  for,  is  contained  in  the  foregoing  description  of  its 
objects.  If  the  objects  be  justifiable,  so  must  also  be  the  length 
which  is  a  necessary  outcome  of  those  objects. 


SECTION    II. 

Quaternion  Theorems. 
Definitions. 

1.  As  there  are  two  or  three  symbols  and  terms  which  will  be 
in  constant  use  in  the  following  pages  that  are  new  or  more  general 
in  their  signification  than  is  usual,  it  is  necessary  to  be  perhaps 
somewhat  tediously  minute  in  a  few  preliminary  definitions  and 
explanations. 

A  function  of  a  variable  in  the  following  essay  is  to  be  under- 
stood to  mean  anything  which  depends  on  the  variable  and  which 
can  be  subjected  to  mathematical  operations ;  the  variable  itself 
being  anything  capable  of  being  represented  by  a  mathematical 
symbol.     In  Cartesian  Geometry  the  variable  is  generally  a  single 
scalar.     In  Quaternions  on  the  other  hand  a  general  quaternion 
variable  is  not  infrequent,  a  variable  which  requires  4  scalars  for 
its  specification,  and  similarly  for  the  function.     In  both,  however, 
either  the   variable  or  the  function   may  be  a  mere  symbol   of 
operation.     In  the   following  essay  we   shall   frequently  have   to 
speak  of  variables  and  functions  which  are  neither  quaternions  nor 
mere  symbols  of  operation.    For  instance  K  in  §  40  below  requires 
6  scalars  to  specify  it,  and  it  is  a  function  of  yfr  which  requires 
G  scalars  and  p  which  requires  3  scalars.      When  in  future   the 
expression  "any  function"  is  used  it  is  always  to  be  understood  in 
the  general  sense  just  explained. 

We  shall  frequently  have  to  deal  with  functions  of  many  inde- 
pendent vectors,  and  especially  with  functions  which  are  linear  in 
each  of  the  constituent  vectors.  These  functions  merely  require 
to  be  noticed  but  not  defined. 

Hamilton  has  defined  the  meaning  of  the  symbolic  vector  V 
thus  : — 

_      .  d       .  d      .  <! 
dx     J  dy      "  dz' 
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where  i,  j,  k  are  unit  vectors  in  the  directions  of  the  mutually 
perpendicular  axes  x,  y,  z.  I  have  found  it  necessary  somewhat  to 
expand  the  meaning  of  this  symbol.  When  a  numerical  suffix 
1,  2,... is  attached  to  a  V  in  any  expression  it  is  to  indicate  that  the 
differentiations  implied  in  the  V  are  to  refer  to  and  only  to  other 
symbols  in  the  same  expression  which  have  the  same  suffix.  After 
the  implied  differentiations  have  been  performed  the  suffixes  are 
of  course  removed.  Thus  Q  (a,  /9,  <y,  8)  being  a  quaternion  func- 
tion of  any  four  vectors  a,  /3,  7,  8,  linear  in  each 


Q  (V2viv*) 

n  /dX  da 
\dx  dx 

a  j  + 

~  /dX  da  . .  \ 
Q\dydxll)  + 

■«(££" 

n  /dX  d/j, 
+  2\didy 

v)  + 

~.  (dX  da  ..  \ 
Qtdydy*>)  + 

«££*' 

\ax  dz 

m  + 

«(££*)  + 

«££»: 

and  again 

Q,  (\>  /V  ?„  V2)  =  Q3  (\,  a2,  Vt  +  V„  V2). 

It  is  convenient  to  reserve  the  symbol  A  for  a  special  meaning. 
It  is  to  be  regarded  as  a  particular  form  of  V,  but  its  differentia- 
tions are  to  refer  to  all  the  variables  in  the  term  in  which  it 
appears.     Thus  Q  being  as  before 

Q  (\,  a,  A,  V:)  =  Qa  (a,,  n.2,  V,  +  V2,  V,)*. 

If  in  a  linear  expression  or  function  Vx  and  p1  (p  being  as  usual 
=  ix  +jy  +  kz)  occur  once  each  they  can  be  interchanged.  Simi- 
larly for  V  and  pr  So  often  does  this  occur  that  I  have  thought 
it  advisable  to  use  a  separate  symbol  £,  for  each  of  the  two  V1  and 
Pl,  £2  for  each  of  the  two  V„  and  p2  and  so  for  £3>  &c.  If  only  one 
such  pair  occur  there  is  of  course  no  need  for  the  suffix  attached  to 
£  Thus  £  may  be  looked  upon  as  a  symbolic  vector  or  as  a  single 
term  put  down  instead  of  three.  '  For  Q  (a,  ft)  being  linear  in  each 
of  the  vectors  a,  /3 

Q  (6  t)  =  Q  (^  ft)  -  Q  ft  0  +  Q  (J J)  +  Q  (*. *) C1)- 

There  is  one  more  extension  of  the  meaning  of  V  to  be  given. 

*  These  meanings  for  Vd  V2— a  *  useA  iu  a  PaPer  on  "Some  General  Theorems 
in  Quaternion  Integration,"  in  the  Mess,  of  Math.  Vol.  xiv.  (1884),  p.  26.  The 
investigations  there  given  are  for  the  most  part  incorporated  below.  [Note  added, 
1892,  see  preface  as  to  the  alteration  of  y'  mto  A.] 
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u,  v,  w  being  the  rectangular  coordinates  of  any  vector  o-,  ^V  is 
defined  by  the  equation 

_      .  d        .  d      7    d 
du        dv        diu 

To  o-V  of  course  are  to  be  attached,  when  necessary,  the  suffixes 
above  explained  in  connection  with  V.  Moreover  just  as  for 
Vj,  px  we  may  put  £,  £  so  also  for  <7V1,  <t1  may  we  put  the  same. 

With  these  meanings  one  important  result  follows  at  once. 
The  Vj's,  V2's,  dkc,  obey  all  the  laws  of  ordinai^y  vectors  whether 
with  regard  to  multiplication  or  addition,  for  the  coordinates 
d/dx,  d/dy,  d\dz  of  any  V  obey  with  the  coordinates  of  any  vector 
or  any  other  V  all  the  laws  of  common  algebra. 

Just  as  o-V  may  be  defined  as  a  symbolic  vector  whose  coordi- 
nates are  djdu,  d/dv,  d/dw  so  cf>  being  a  linear  vector  function  of 
any  vector  whose  coordinates  are 

(ftAci  «AC2  «363c3)  (i-e-  &  =  axi  +  hJ  +  cik>  &c-)- 
0Q*  is  defined  as  a  symbolic  linear  vector  function  whose  coordi- 
nates are 

(d/dav  d/db1}  d/dcv  d/da2,  d/db2,  d/dc2,  d/da3,  d/db3,  djdc3), 

and  to  ^Q  is  to  be  applied  exactly  the  same  system  of  suffixes  as 
in  the  case  of  V.  Thus  q  being  any  quaternion  function  of  <£,  and 
to  any  vector 

^Qjft)  ■q1  =  —  (idq/da1  +jdqjdb1  +  Mq/dcJ  Si<o 

—  (idq/da2  +jdq/db2  +  kdq/dc2)  Sjco 

—  (idq/da3  +jdq/db3  +  hdqjdc3)  Ska. 

The  same  symbol  ^(j  is  used  without  any  inconvenience  with  a 
slightly  different  meaning.  If  the  independent  variable  </>  be  a 
self -conjugate  linear  vector  function  it  has  only  six  coordinates. 
If  these  are  PQRSTU  (i.e.  <j>i  =  Pi  +  Uj  +  Tk,  &c.)  0(i  is  defined 
as  a  self- conjugate  linear  vector  function  whose  coordinates  are 

(d/dP,  d/dQ,  d/dR,  %d/dS,  \djdT,  hd/dU). 

We  shall  frequently  have  to  compare  volume  integrals  with 
integrals  taken  over  the  bounding  surface  of  the  volume,  and 
again  surface  integrals  with  integrals  taken  round  the  boundary 
of  the  surface.  For  this  purpose  we  shall  use  the  following  nota- 
tions for  linear,  surface  and  volume  integrals  respectively  JQdp, 

*  I  have  used  an  inverted  I)  to  indicate  the  analogy  to  Hamilton's  inverted  A. 
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f$Qa"S,  ffJQdS  where  Q  is  any  function  of  the  position  of  a  point. 
Here  dp  is  a  vector  element  of  the  curve,  dX  a  vector  element  of 
the  surface,  and  d%  an  element  of  volume.  When  comparisons 
between  line  and  surface  integrals  are  made  we  take  d%  in  such  a 
direction  that  dp  is  in  the  direction  of  positive  rotation  round  the 
element  d%  close  to  it.  When  comparisons  between  surface  and 
volume  integrals  are  made  dX  is  always  taken  in  the  direction 
away  from  the  volume  which  it  bounds. 


Properties  of  £ 

2.  The  property  of  £  on  which  nearly  all  its  usefulness 
depends  is  that  if  a  be  any  vector 

o-  =  -  giSfgir, 

which  is  given  at  once  by  equation  (1)  of  last  section. 

This  gives  a  useful  expression  for  the  conjugate  of  a  linear 
vector  function  of  a  vector.  Let  <f>  be  the  function  and  co,  r  any 
two  vectors.  Then  <p'  denoting  as  usual  the  conjugate  of  $  we 
have 

SwcpT  =  Srcp'a), 

whence  putting  on  the  left  t  =  —  §S§t  we  have 

St  (-  {SaxpO  =  Srcp'co, 

or  since  t  is  quite  arbitrary 

<f)'cD  =  -ZSco<l>£ (2). 

From  this  we  at  once  deduce  expressions  for  the  pure  part  (pco 
and  the  rotational  part  Veto  of  (pco  by  putting 

(0  +  <f>' )  ft)  =  -  tfSnt  -  {ScocpZ  =  Zcp~co  | 

(cp  -  <p')  to  =  -  0$Sfft>f  +  ?&»><££=  VVty^.to  =  2  Few  | ' 

And  all  the  other  well-known  relations  between  </>  and  0'  are  at 
once  given  e.g.  Sty£  =  Sty '£,  i.e.  the  "convergence"  of  </>  =  the 
"  convergence  "  of  cp'. 

3.  Let  Q  (X,  jjl)  be  any  function  of  two  vectors  which  is  linear 
in  each.  Then  if  (pco  be  any  linear  vector  function  of  a  vector  w 
given  by 

</)&)  =  -  XfiScoa  \  /4x 

we  have  Q(£  «  =  SQ(-  «a,/3)  =  2Q(«,0)  f "'" 
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or  more  generally 

Q(£to£)  =  2Q(%'«,/3) W. 

To  prove,  it  is  only  necessary  to  observe  that 

and  that  —  %S£xa  =  Xa- 

As  a  particular  case  of  eq.  (4)  let  <f>  have  the  self-conjugate  value 

4>(0  —  —  \  X  (fiScoa  +  aS(o{3)               \  ,  , 

Then  Q(r,0C)  =  iS{Q(a,/3)  +  Q(Aa))J (5)' 

or  if  Q  (X,  fj,)  is  symmetrical  in  X  and  /x 

Q(£<K)  =  2Q(a,/3)    (6). 

The  application  we  shall  frequently  make  of  this  is  to  the  case 
when  for  a  we  put  V1  and  for  fi,  <tx,  where  er  is  any  vector  function 
of  the  position  of  a  point.  In  this  case  the  first  expression  for  cf> 
is  the  strain  function  and  the  second  expression  the  pure  strain 
function  resulting  from  a  small  displacement  cr  at  every  point. 
As  a  simple  particular  case  put  Q  (X,  /x)  —  S\/x  so  that  Q  is  sym- 
metrical in  X  and  fx.     Thus  <f)  being  either  of  these  functions 

Another  important  equation  is 

•Q(fc#)=Q(*'fcO  (6a). 

This  is  quite  independent  of  the  form  of  </>.  To  prove,  observe 
that  by  equation  (2) 

and  that  —  £7^$^  =  (££,.  Thus  we  get  rid  of  £  and  may  now 
drop  the  suffix  of  ^  and  so  get  eq.  (6a).  [Notice  that  by  means 
of  (6a),  (4a)  may  be  deduced  from  (4);  for  by  (6a) 

Q  (fc  fat)  =  Q  (tffc  W  =  SQ  (x'cc,  fi)  by  (4).] 

3a.  A  more  important  result  is  the  expression  for  (fflco  in 
ten i is  of  <f>.     We  assume  that 

yS<p\cf)fi(f)i>  =  niSXfiv, 

where  X,  /x,  v  are  any  three  vectors  and  m  is  a  scalar  independent 

*  [Note  added,  1892.  For  practice  it  is  convenient  to  remember  this  in  words  : — 
.!  term  in  which  f  and  (j>i;  occur  is  unaltered  in  value  by  changing  them  into  <p'£ 
and  £  respectively."] 
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of  these  vectors.     Substituting  £,  f2,  f3  for  X,,  //,,  y  and  multiplying 

stUMWM,  =  v™ (66), 

which  gives  w  in  terms  of  <f>.     That  /S^^^S'^^^,  =  6  is  seen  by 
getting  rid  of  each  pair  of  fs  in  succession  thus : — 

as  (^  Fr2o  &3 = -  sf?& .  r2^3  =  s  (xx  -  wsw  & 

=  -  »$&  =  6. 

Next  observe  that 

Multiplying  by  V££2  and  again  on  the  right  getting  rid  of  the 
£"s  we  have 

V^Scf)cocf)^2  =  -2ma> (6c), 

whence  from  equation  (66) 
or  changing  &>  into  <£_1&) 

^-%Wl«    W 

By  equation  (6a)  of  last  section  we  can  also  put  this  in  the  form 

so  that  ^>-1&)  is  obtained  explicitly  in  terms  of  <£  or  <£'. 

Equation  (6c)  or  (Qd)  can  be  put  in  another  useful  form  which 
is  more  analogous  to  the  ordinary  cubic  and  can  be  easily  deduced 
therefrom,  or  *less  easily  from  (6d),  viz. 

=  0 (6/). 

*  [Note  added,  1892.     The  cubic  may  be  obtained  in  a  more  useful  form  from 
the  equation  wSfIf2f3^|-10f2^f3=  -  3Ff1f2Sty>«0f1<#'2  thus 

=  <pu>S .  tff^fjFftfa  -  20frS' .  0U0&F&&. 
Again  0frSf .  0w0f2Ff1f2=0f15.  0uT'.  tfsFft-a=tfiS*«  (  ~  fiS&^ft  +  frWi^ft) 

=  -  0  (fi<Sffi0«)  Sft*ft  +  0  (frS  Jifli)  «■#« 

=  jPaSfyf-  (p-ZS&w  =  <f>-uS  ?<!>£+  <p3u>. 

Hence  03w  —  m"<p2w  +  ni'<pw  —mu= 0, 

where  6m = S&  ^V^i^a^s 

2m'=-SFf,f2F0fl(/.f, 

M.  2 
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As  a  useful  particular  case  of  equation  (6d)  we  may  notice 
that  by  equation  (4)  §  3  if 

(fray  =  —  cr1$ft>V] 

^a)  =  -  3  vvy^ajsvy^sSv^ (%), 

and  #-lo)  =  -  3  Va^SaVyjSVyy^a^ (6A). 

4.     Let   0,   -\jr   be   two   linear  vector   functions   of  a    vector. 
Then  if 

where  %  is  a  quite  arbitrary  linear  vector  function 

#  =  ^> 

for  we  may  put  ^=tS^(o  where  t  and  w  are  arbitrary  vectors,  so 
that 

&T(f)0)  =  Sryjr(i), 

or  (j)(o  =  yfrco. 

Similarly*  if  <£  and  yjr  are  both  self-conjugate  and  ^  is  a  quite 
arbitrary  self-conjugate  linear  vector  function  the  same  relation 
holds  as  can  be  seen  by  putting 

•  X%~  T^^Q)  +  0)$%T. 


Fundamental  property  of  (j. 

5.     Just  as  the  fundamental  property  of  ^V  is  that,  Q  being 
any  function  of  a 

ZQ  =  -QtS8a(TVi> 

so  we  have  a  similar  property  of  q.  Q  being  any  function  of  </>  a 
linear  vector  function 

BQ  =  -Ql884>balS  (7). 

The  property  is  proved  in  the  same  way  as  for  V,  viz.  by  ex- 
panding 8B<f>^([^  in  terms  of  the  coordinates  of  ^q.  First  let  <£ 
be  not  self-conjugate,  and  let  its  nine  coordinates  be 

(aAciaA°2aAc8)- 

*  [Note  added,  1892.  The  following  slightly  more  general  statement  is  a  practi- 
cally much  more  convenient  form  of  enunciation  :  if  Sx&£  =  Sxfff,  where  x  is  a 
perfectly  ^arbitrary  self-conjugate  and  </>  and  f  are  not  necessarily  self-conjugate 
then  0  —  i/7]. 
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Thus 

=  BaxdQldax  +  Bb^Qjdb^  +  Sc.dQ/dc,, 
+  Ba2dQ/da2  +  Bb^dQ/db^  +  Bc2dQ/dc2, 
+  BatdQ/daa  +  BbzdQ/db3  +  Bc3dQ/dc3  =  BQ. 

The  proposition  is  exactly  similarly  proved  when  <£  is  self-con- 
jugate. 


Theorems  in  Integration. 

6.  Referring  back  to  §  1  above  for  our  notation  for  linear 
surface  and  volume  integrals  we  will  now  prove  that  if  Q  be  any 
linear  function  of  a  vector* 

JQdp  =  ffQ (VdtA)  (8), 

ffQdZ  =  JJJQAds (9). 

To  prove  the  first  divide  the  surface  up  into  a  series  of 
elementary  parallelograms  by  two  families  of  lines — one  or  more 
members  of  one  family  coinciding  with  the  given  boundary, — 
apply  the  line  integral  to  the  boundary  of  each  parallelogram  and 
sum  for  the  whole.  The  result  Avill  be  the  linear  integral  given  in 
equation  (8).  Let  the  sides  of  one  such  parallelogram  taken  in 
order  in  the  positive  direction  be  a,  ft  +  ft',  —  a  —  a',  —  ft ;  so  that 
a'  and  ft'  are  infinitely  small  compared  with  a  and  ft,  and  we  have 
the  identical  relation 

0  =  a  +  ft  +  ft'-a-*-ft  =  ft'-a. 

The  terms  contributed  to  JQdp  by  the  sides  a  and  —a  —  a  will 
be  (neglecting  terms  of  the  third  and  higher  orders  of  small 
quantities) 

Qa-Qa-  Qa  +  QlaSftV1  =  -  Qa  +  Q.aS/SV,. 

Similarly  the  terms  given  by  the  other  two  sides  will  be 

Qft'  -  Q1ftSaVl 

*  These  two  propositions  are  generalisations  of  what  Tait  and  Hicks  have  from 
time  to  time  proved.  They  were  first  given  in  the  present  form  by  me  in  the 
article  already  referred  to  in  §  1  above.  In  that  paper  the  necessary  references  are 
given. 

2—2 
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so  that  remembering  that  /3'  -  a  =  0  and  therefore  Qff-Qa.'  =  0 
we  have  for  the  whole  boundary  of  the  parallelogram 

Q,*S(3V1  -  QfiSaV,  =  Q^Wafi .  V,)  =  QVdZA, 
where  dTS  is  put  for  Fa/3.     Adding  for  the  whole  surface  we  get 
equation  (8). 

Equation  (9)  is  proved  in  an  exactly  similar  way  by  splitting 
the  volume  up  into  elementary  parallelepipeda  by  three  families 
of  surfaces  one  or  more  members  of  one  of  the  families  coinciding 
with  the  given  boundary.  If  a,  /3,  y  be  the  vector  edges  of  one 
such  parallelepiped  we  get  a  term  corresponding  to  Q@'  —  Q/x  viz. 

Q  (vector  sum  of  surface  of  parallelepiped)  =  0, 

and  we  get  the  sum  of  three  terms  corresponding  to 

QiaS/3V1-Qi(3ScNl   . 
above,  viz. 

-  Qt  ( VPy)  S*V1  -  Qt  (  VyOL)  S/3V1  -  Q1  ( Fa/3)  8yVt  =  -  Q^Safiy, 

whence  putting  Safty  =  —  d%  we  get  equation  (9). 

7.  It  will  be  observed  that  the  above  theorems  have  been 
proved  only  for  cases  where  we  can  put  dQ  =  —  Q1  SdpV1  i.e.  when 
the  space  fluxes  of  Q  are  finite.  If  at  any  isolated  point  they  are 
not  finite  this  point  must  be  shut  off  from  the  rest  of  the  space  by 
a  small  closed  surface  or  curve  as  the  case  may  be  and  this  surface 
or  curve  must  be  reckoned  as  part  of  the  boundary  of  the  space. 
If  at  a  surface  (or  curve)  Q  has  a  discontinuous  value  so  that  its 
derivatives  are  there  infinite  whereas  on  each  side  they  are 
finite,  this  surface  (or  curve)  must  be  considered  as  part  of  the 
boundary  and  each  element  of  it  will  occur  twice,  i.e.  once  for  the 
part  of  the  space  on  each  side. 

In  the  case  of  the  isolated  points,  if  the  surface  integral  or  line 
integral  round  this  added  boundary  vanish,  we  can  of  course  cease 
to  consider  these  points  as  singular.  Suppose  Q  becomes  infinite 
at  the  point  p  =  a.  Draw  a  small  sphere  of  radius  a  and  also  a 
sphere  of  unit  radius  with  the  point  a  for  centre,  and  consider  the 
small  sphere  to  be  the  added  boundary.  Let  d%'  be  the  element 
of  the  unit  sphere  cut  off  by  the  cone  which  has  a.  for  vertex  and 
the  element  c£S  of  the  small  sphere  for  base.  Then  d?L  =  a2dX' 
and  we  get  for  the  part  of  the  surface  integral  considered  a2JjQ^d%' 
where  QdK  is  the  value  of  Q  at  the  element  d%.     If  then 

UT(p-a)^T*(j>-a)QdSU(p-a)  =  Q 
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the  point  may  be  regarded  as  not  singular.  If  the  limiting  ex- 
pression is  finite  the  added  surface  integral  will  be  finite.  If  the 
expression  is  infinite  the  added  surface  integral  will  be  generally 
but  not  always  infinite.  Similarly  in  the  case  of  an  added  line 
integral  if  Lt  T  (p_a)  =0  T(p  -  a.)  Q  U(p  -  a)  is  zero  or  finite,  the  added 
line  integral  will  be  zero  or  finite  respectively  (of  course  including 
in  the  term  finite  a  possibility  of  zero  value).  If  this  expression 
be  infinite,  the  added  line  integral  will  generally  also  be  infinite. 

This  leads  to  the  consideration  of  potentials  which  is  given 
in  §  9. 

8.  Some  particular  cases  of  equations  (8)  and  (9)  which 
(except  the  last)  have  been  proved  by  Tait,  are  very  useful. 
First  put  Q  =  a  simple  scalar  P.     Thus 

jPdP=jjVd%VP (10), 

jjPdt=jjjVPd% (ii). 

If  P  be  the  pressure  in  a  fluid  —JJPdX  is  the  force  resulting  from 
the  pressure  on  any  portion  and  equation  (11)  shews  that  —  VP  is 
the  force  per  unit  volume  due  to  the  same  cause.  Next  put 
Qa>  =  Sua  and  Vcoa.     Thus 

fSdpa=jfSdZVa (12), 

jVdpa  =  fJV(Vd%V  .  a)  =JJdSSVa  -  JJ^fltfSo^.flS), 

ffSdZ*  =  fffSV*d% (14), 

JjVdSa  =  JjfVVad% (15). 

Equations  (12)  and  (14)  are  well-known  theorems,  and  (13)  and 
(15)  will  receive  applications  in  the  following  pages.  Green's 
Theorem  with  Thomson's  extension  of  it  are,  as  indeed  has  been 
pointed  out  by  Tait  particular  cases  of  these  equations. 

Equations  (14)  and  (15)  applied  to  an  element  give  the  well- 
known  physical  meanings  for  S^a  and  FV<r*.  The  first  is 
obtained  by  applying  (14)  to  any  element,  and  the  second  (regard- 
ing irasa  velocity)  is  obtained  by  applying  (15)  to  the  element 
contained  by  the  following  six  planes  each  passing  infinitely  near 
to  the  point  considered — (1)  two  planes  containing  the  instanta- 
neous axis  of  rotation,  (2)  two  planes  at  right  angles  to  this  axis, 
and  (3)  two  planes  at  right  angles  to  these  four. 

*  [Note  added,  1892.  Let  me  disarm  criticism  by  confessing  that  what  follows 
concerning  Vya  is  nonsense.] 


22  QUATERNION  THEOREMS.  [§  8. 

One  very  frequent  application  of  equation  (9)  may  be  put  in 
the  following  form: — Q  being  any  linear  function  (varying  from 
point  to  point)  of  i2t  and  V1?  R  being  a  function  of  the  position  of 
a  point 

MQ(R1,V1)dz  =  -fffQ1(B,V1)d%  +  fJQ(R,d$)...(lG). 


Potentials. 

9.  We  proceed  at  once  to  the  application  of  these  theorems 
in  integration  to  Potentials.  Although  the  results  about  to  be 
obtained  are  well-known  ones  in  Cartesian  Geometry  or  are  easily 
deduced  from  such  results  it  is  well  to  give  this  quaternion  method 
if  only  for  the  collateral  considerations  which  on  account  of  their 
many  applications  in  what  follows  it  is  expedient  to  place  in  this 
preliminary  section. 

If  R  is  some  function  of  p  -  p  where  p  is  the  vector  coordinate 
of  some  point  under  consideration  and  p  the  vector  coordinate  of 
any  point  in  space,  we  have 

pVi2  =  -p,ViZ. 

Now  let  Q(R)  be  any  function  of  R,  the  coordinates  of  Q  being- 
functions  of  p  only.  Consider  the  integral  JjjQ  (R)  d%  the  variable 
of  integration  being  p  (p  being  a  constant  so  far  as  the  integral  is 
concerned).  It  does  not  matter  whether  the  integral  is  a  volume, 
surface  or  linear  one  but  for  conciseness  let  us  take  it  as  a  volume 
integral.     Thus  we  have 

Now  pV  operating  on  the  whole  integral  has  no  meaning  so  we 
may  drop  the  affix  to  the  V  uutside  and  always  understand  p. 
Under  the  integral  sign  however  we  must  retain  the  affix  p  or  p 
unless  a  convention  be  adopted.  It  is  convenient  to  adopt  such  a 
convention  and  since  Q  will  probably  contain  some  pV  but  cannot 
possibly  contain  a  P-V  we  must  assume  that  when  V  appears  with- 
out an  affix  under  the  integral  sign  the  affix  p  is  understood. 
With  this  understanding  we  see  that  when  V  crosses  the  integral 
sign  it  must  be  made  to  change  sign  and  refer  only  to  the  part  we 
have  called  R.     Thus 

WQ(&)ds  =  -fffVlQ(Rl)d%  (17). 
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Generally  speaking  R  can  and  will  be  put  as  a  function  of 
T~l  {p  —  p)  and  for  this  we  adopt  the  single  symbol  a.  Both  this 
symbol  and  the  convention  just  explained  will  be  constantly  re- 
quired in  all  the  applications  which  follow. 

10.  Now  let  Q  be  any  function  of  the  position  of  a  point. 
Then  the  potential  of  the  volume  distribution  of  Q,  say  q,  is  given 
by:- 

q  =  MuQd% (18), 

the  extent  of  the  volume  included  being  supposed  given.  We 
may  now  prove  the  following  two  important  propositions 

V22  =  4ttQ (19), 

tfSdSV.q^lnrfflQfo (20). 

The  latter  is  a  corollary  of  the  former  as  is  seen  from  equation  (9) 
§  6  above. 

Equation  (19)  may  be  proved  thus.  If  P  be  any  function  of 
the  position  of  a  point  which  is  finite  but  not  necessarily  continu- 
ous for  all  points 

VjOJuPcfe, 

is  always  finite  and  if  the  volume  over  which  the  integral  extends 
is  indefinitely  diminished,  so  also  is  the  expression  now  under 
consideration,  and  this  for  the  point  at  which  is  this  remnant  of 
volume.  This  in  itself  is  an  important  proposition.  The  expres- 
sion, by  equation  (17),  =  —  fff¥uPd%  and  both  statements  are 
obviously  true  except  for  the  point  p.  For  this  point  we  have 
merely  to  shew  that  the  part  of  the  volume  integral  just  given 
contributed  by  the  volume  indefinitely  near  to  p  vanishes  with 
this  volume.  Divide  this  near  volume  up  into  a  series  of  element- 
ary cones  with  p  for  vertex.  If  r  is  the  (small)  height  and  dco  the 
solid  vertical  angle  of  one  of  these  cones,  the  part  contributed  by 
this  cone  is  approximately  U(p  —  p')Pprdai/S  where  Pp  is  the 
value  of  P  at  the  point  p >'.     The  proposition  is  now  obvious. 

Now  since 

V2q  =  V2fff»Qd%  =  fffV*uQd*, 

we  see  that  the  only  part  of  the  volume  integral  JjfuQdS  which 
need  be  considered  is  that  given  by  the  volume  in  the  immediate 
neighbourhood  of  p,  for  at  all  points  except  p,  V-u  =  0.  Consider 
then  our  volume  and  surface  integrals  only  to   refer  to  a  small 
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sphere  with  p  for  centre  and  so  small  that  no  point  is  included  at 
which  Q  is  discontinuous  and  therefore  VQ  infinite.  (This  last 
assumes  that  Q  is  not  discontinuous  at  p.  In  the  case  when  Q  is 
discontinuous  at  p  no  definite  meaning  can  be  attached  to  the 
expression  V2</.)     We  now  have 

V2?  =  V*JJfuQd%  =  -  VffJVuQd*  [by  §  9] 
=  V//JuVQ<fe  -  VffudZQ, 

equation  (9)  §  6  being  applied  and  the  centre  of  the  sphere  being 
not  considered  as  a  singular  point  since  the  condition  of  §  7  is 
satisfied, viz.  that  Lt  t(p-p')=o  T2 (p  —  p) uU (p  —  p)  =  0.  Now  put- 
ting P  above  =  VQ  we  see  that  the  first  expression,  viz.  VffjiiVQd% 
can  be  neglected  and  the  second  gives 

V2q  =  ffVudZQ  =  4ttQ, 

where  Q  is  the  mean  value  of  Q  over  the  surface  of  the  sphere  and 
therefore  in  the  limit  =  Q.     Thus  equation  (19)  has  been  proved. 

When  Q  has  a  simple  scalar  value  all  the  above  propositions, 
and  indeed  processes,  become  well-known  ones  in  the  theory  of 
gravitational  potential. 

We  do  not  propose  to  go  further  into  the  theory  of  Potentials 
as  the  work  would  not  have  so  direct  a  bearing  on  what  follows  as 
these  few  considerations. 


SECTION   III. 

*  Elastic  Solids. 

Brief  recapitulation  of  previous  work  in  this  branch. 

11.  As  far  as  I  am  aware  the  only  author  who  has  applied 
Quaternions  to  Elasticity  is  Prof.  Tait.  In  the  chapter  on  Kine- 
matics of  his  treatise  on  Quaternions,  §§  360 — 371,  he  has  considered 
the  mathematics  of  strain  with  some  elaboration  and  again  in  the 
chapter  on  Physical  Applications,  §§  487 — 491,  he  has  done  the 
same  with  reference  to  stress  and  also  its  expression  in  terms  of 
the  displacement  at  every  point  of  an  elastic  body. 

In  the  former  he  has  very  successfully  considered  various 
useful  decompositions  of  strain  into  pure  and  rotational  parts  and 
so  far  as  strain  alone  is  considered,  i.e.  without  reference  to  what 
stress  brings  it  about  he  has  left  little  or  nothing  to  be  done.  In 
the  latter  he  has  worked  out  the  expressions  for  stress  by  means 
of  certain  vector  functions  at  each  point,  which  express  the  elastic 
properties  of  the  body  at  that  point. 

But  as  far  as  I  can  see  his  method  will  not  easily  adapt  itself 
to  the  solution  of  problems  which  have  already  been  considered 
by  other  methods,  or  prepare  the  way  for  the  solution  of  fresh 
problems.  To  put  Quaternions  in  this  position  is  our  present 
object.  I  limit  myself  to  the  statical  aspect  of  Elasticity,  but  I 
believe  that  Quaternions  can  be  as  readily,  or  nearly  so,  applied  to 
the  Kinetics  of  the  subject. 

For  the  sake  of  completeness  I  shall  repeat  in  my  own  no- 
tation a  small  part  of  the  work  that  Tait  has  given. 

*  [Note  added,  1892.  It  would  be  better  to  head  this  section  "  Elastic  bodies  " 
since  except  when  the  strains  are  assumed  small  the  equations  are  equally  true  of 
solids  and  fluids.  I  may  say  here  that  I  have  proved  in  the  Proc.  li.  S.  E.  1890 
— 91,  pp.  106. et  seq.,  most  of  the  general  propositions  of  this  section  somewhat 
more  neatly  though  the  processes  are  essentially  the  same  as  here.] 
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Tait  shews  (§  370  of  his  Quaternions)  that  in  any  small  portion 
of  a  strained  medium  the  strain  is  homogeneous  and  (§  360)  that 
a  homogeneous  strain  function  is  a  linear  vector  one.  He  also 
shews  (§  487)  that  the  stress  function  is  a  linear  vector  one  and  he 
obtains  expressions  (§§  487 — 8)  for  the  force  per  unit  volume  due 
to  the  stress,  in  terms  of  the  space- variation  of  the  stress. 


Strain,  Stress-force,  Stress-couple. 

12.  This  last  however  I  give  in  my  own  notation.  His 
expression  in  §  370  for  the  strain  function  I  shall  throughout 
denote  bv  y.     Thus 

XU>=  W  —  So)W  .  7/ (1), 

where  7]  is  the  displacement  that  gives  rise  to  the  strain. 

Let  %  consist  of  a  pure  strain  ^  followed  by  a  rotation  q  (  )q~v 
as  explained  in  Tait's  Quaternions,  §  365  where  he  obtains  both  q 
and  ty  in  terms  of  %.     Thus 

Xo>  =  qfcoq-1 (2). 

When  the  strain  is  small  i/r  takes  the  convenient  form  %  where 
X  stands  for  the  pure  part  of  x  so  that 

XO)  =  co  -  hScoV  .  7)  -  %V1San]1 (3), 

by  equation  (3)  §  2  above.  Similarly  q  (  )  q'1  becomes  V0  {  )  where 
20  =  FV77.  The  truth  of  these  statements  is  seen  by  putting  in 
equation  (2)  for  q,  1  +  6/2  and  therefore  for  q'1,  1  —  6/2  for  ■ty,  % 
and  then  neglecting  all  small  quantities  of  an  order  higher  than 
the  first. 

13.  Next  let  us  find  the  force  and  couple  per  unit  volume  due 

to   a   stress  which  varies    from  point  to    point.     Let    the    stress 

function  be  cf>.     Then  the  force  on  any  part  of  the  body,  due  to 

stress,  is 

jfydS  =///</>  Ads, 

by  equation  (9)  §  6.  Thus  the  force  per  unit  volume  =  <f>A,  for 
the  volume  considered  in  the  equation  may  be  taken  as  the 
element  d%. 

Again  the  moment  round  any  arbitrary  origin  is 
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by  equation  (9)  §  6.  The  second  term  on  the  right  is  that  due  to 
the  force  </>A  just  considered,  and  the  first  shews  that  in  addition 
to  this  there  is  a  couple  per  unit  volume  =  V£<f)£  =  Uvice  the 
"  rotation  "  vector  of  0.  Let  then  <f>  be  the  pure  part  and  Ve  (  ) 
the  rotatory  part  of  </>.     Thus 

Force  per  unit  volume 

=  0A  =  0A-  FVe    (4-), 

Couple  per  unit  volume 

=  V&Z=2e (5). 

These  results  are  of  course  equivalent  to  those  obtained  by 
Tait,  Quaternions,  §§  487 — 8. 

The  meanings  just  given  to  ?;,  %,  ^,  ty,  q,  (£,  <£  and  e  will  be 
retained  throughout  this  Section.  In  air  cases  of  small  strain  as 
we  have  seen  we  may  use  ty  or  ^  indifferently  and  whenever-  we 
wish  to  indicate  that  we  are  considering  the  physical  phenomenon 
of  pure  strain  we  shall  use  yjr,  %  being  regarded  merely  as  a 
function  of  %.  We  shall  soon  introduce  a  function  zr  which  will 
stand  towards  <f>  somewhat  as  ty  towards  ^  and  such  that  when 
the  strain  is  small  vr  =  <£. 

It  is  to  be  observed  that  </>&)  is  the  force  exerted  on  a  vector 
area,  which  when  strained  is  &>,  not  the  stress  on  an  area  which 
before  strain  is  &>.  Similarly  in  equations  (4)  and  (5)  the  in- 
dependent variable  of  differentiation  is  p  +  rq  so  that  strictly 
speaking  in  (4)  we  should  put  (/>P+7,A  —  Vp+1)Ve.  In  the  case  of 
small  strain  these  distinctions  need  not  be  made. 


Stress  in  terms  of  strain. 

14.  To  express  stress  in  terms  of  strain  we  assume  any 
displacement  and  consequent  strain  at  every  point  of  the  body 
and  then  give  to  every  point  a  small  additional  displacement  8i] 
and  find  in  terms  of  yjr  and  $  the  increment  JjJ8wd%  in  the 
potential  energy  of  the  body,  iudS0  being  the  potential  energy  of 
any  element  of  the  body  whose  volume  before  strain  was  rfg0.    Thus 

ffj8wd%0  =  (work  done  by  stresses  on  surface  of  portion  considered), 

—  (work    done    by   stresses    throughout    volume    of  same 
portion). 
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Thus,  observing  that  by  §  12  the  rotation  due  to  the  small  dis- 
placement Brj  is  Vp+vVSr]/2,  we  have 

The  first  of  the  terms  on  the  right  is  the  work  done  on  the  surface 
of  the  portion  of  the  body  considered ;  the  second  is  —  (work  done 
by  stress-forces  <£p+r,A) ;  and  the  third  is  -  (work  done  by  stress- 
couples  2e).  Thus  converting  the  surface  integral  into  a  volume 
integral  by  equation  (9)  §  6  above 

Limiting  the  portion  of  the  body  considered  to  the  element  d%  we 
get 

8w  =  -mS8Vl(j>p+vV1 (6) 

where  m  is  put  for  d%/d%0  and  therefore  may  be  put  by  §§  3a  and  3 
above  in  the  various  forms 

Gm-aZt&SxtxZxS, (6a) 

-^Wfe^WJi   (66) 

=^r2r8^Kf^r8  c&o 

=  ^W7, 8  (P  +  v\  (p  +  vX  (p  +  v)3 (6d). 

It  is  to  be  observed  that  since  the  rotation-vector  e  of  cp  does  not 
occur  in  equation  (6),  e  and  therefore  the  stress-couple  are  quite  arbi- 
trary so  far  as  the  strain  and  potential  energy  are  concerned.  I  do 
not  know  whether  this  has  been  pointed  out  before.  Of  course 
other  data  in  the  problem  give  the  stress-couple.  In  fact  it  can  be 
easily  shewn  that  in  all  cases  whether  there  be  equilibrium  or  not 
the  external  couple  per  unit  volume  balances  the  stress-couple. 
[Otherwise  the  angular  acceleration  of  the  element  would  be 
infinite.]  Thus  if  M  be  the  external  couple  per  unit  volume  of 
the  unstrained  solid  we  have  always 

M+2me  =  0 (7). 

In  the  particular  case  of  equilibrium  F  being  the  external  force 
per  unit  volume  of  the  unstrained  solid  we  have 

F  +  m(0p+,A-Fp+,Ve)  =  O (8). 

The  mathematical  problem  is  then  the  same  as  if  for  F/m  we 
substituted  F/m  +  Fp+J)V  (M/2m)  and  for  M,  zero.     In  the  case  of 
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small  strains  m  may  be  put  =  1.  In  this  case  then  the  mathe- 
matical problem  is  the  same  as  if  for  F  we  substituted  P  +  FVM/2 
and  for  M  zero. 

Returning  to  equation  (6)  observe  that 

Sxo  =  -  Sxo)p+11V  .  8r] 

[which  is  established  just  as  is  the  equation  xw  =  —  ScoV  .  (p  +  77)], 
so  that  changing  co  which  is  any  vector  into  %-1eo 

Therefore  by  equation  (6)  of  this  section  and  equation  (4)  §  3  above 
we  have 

8w  =  —  mSSxx'1^^ 

We  must  express  the  differential  on  the  right  of  this  equation  in 
terms  of  S-v/r  and  8q,  the  latter  however  disappearing  as  we  should 
expect.  Now  ^eo  =  qtyoocf1  [equation  (2)  §  12]  so  that  remember- 
ing that  8  .  q~l  —  —  q'^Sqq'1 

S^&)  =  8qijrcoq~1  —  qyjrcoq~l  8qq~*  +  qSyjrcoq'1 
=  2VV8qq~1 .  %&)  +  q8yfrcoq~\ 

Sw/m  =  -28.  VSqq'1 .  £$£  -  Sq8fx~^q~lK 

or  since  V£(b%=  0,  cp  being  self-conjugate 

8w  =  —  mS8ylrx~1%q~1<p%q. 

This  can  be  put  into  a  more  convenient  shape  for  our  present 
purpose.  First  put  for  x~l  ^s  value  "fr^q-1  ( )  q  and  then  apply 
equation  (6a)  §  3  above,  putting  for  the  <p  of  that  equation  q~*  (  )  q 
and  therefore  for  the  <p',  q  (  )  q'1.     Thus 

8w  =  —  mS8y\r*tyrX%'&%   (9), 

where  -srco  =  q~*cp  (qcoq-1)  q  (10). 

The  physical  meaning  of  this  last  equation  can  easily  be  shewn. 
Suppose  when  there  is  no  rotation  that  <p  =  ot'.  Then  it  is 
natural*  to  assume — in  fact  it  seems  almost  axiomatic — that  the 
superimposed  rotation  q()  q'1  should  merely  so  to  speak  rotate 
the  stress  along  with  it.  Thus  if  co  is  some  vector  area  before  the 
rotation  which  becomes  co'  by  means  of  the  rotation 

qtx'co  .  q~l  =  cpco'. 

*  Observe  that  we  do  not  make  this  assumption.     We  really  shew  that  it  is 
true. 
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But  &)'  =  qcocf1,  so  that 

tzr'ft)  =  rf1^)  (qwcf1)  q  =  -sTft). 
Thus  Ave  see  that  ct  is  as  it  were  <£  with  the  rotation  undone. 

Before  proceeding  further  with  the  calculation  let  us  see  what 
we  have  assumed  and  what  equation  (9)  teaches  us.  The  one  thing 
we  have  assumed  is  that  the  potential  energy  of  the  body  can  be 
taken  as  the  sum  of  the  potential  energies  of  its  elements,  in  other 
words  that  no  part  of  the  potential  energy  depends  conjointly  on 
the  strains  at  P  and  Q  where  P  and  Q  are  points  separated  by  a 
finite  distance.  This  we  must  take  as  an  axiom.  By  it  we  are 
led  to  the  expression  for  Sw  in  equation  (9).  This  only  involves 
the  variation  of  the  pure  strain  -v|r  but  not  the  space  differential 
coefficients  of  -fy.  This  is  not  an  obvious  result  as  far  as  I  can 
see  but  it  is  I  believe  always  assumed  without  proof. 

We  may  now  regard  w  as  a  function  of  yjr  only.  Therefore  by 
equation  (7)  §  5  above 

therefore  by  equation  (9)  of  this  section 

In  equation  (6a)  §  3  above  putting  <£  —  yjr'1  the  right-hand  member  of 
this  equation  becomes  mSS-^r^yjr'1^.  Now  putting  in  equation  {6a) 
§  3  <f>  =  -stt/t-1  this  member  becomes  mS^Byfr^y1^^  or  mSS-ty-ty'1??^ 
Thus  we  have 

ZSSyfrZ+QwS  =  mSByJr^myjr-1  +  ^~V)  £ 
Now  ot-^""1  4-  ^t_1ot  is  self-conjugate.     Hence  by  §  4  above 

m  (ot^_1  +  i/r"1^)  =  2l/,QW (11). 

This  equation  can  be  looked  upon  as  giving  w  in  terms  of  the 
strain.  We  can  obtain  -cr  however  explicitly  for  -^"Sb-  is  the  conju- 
gate of  -oji/r-1.  Hence  from  equation  (11)  (because  {^ta  +  •cj'\^~1)/2 
is  the  pure  part  of  ota/t-1) 

mTz\\r~lw  —  $(jwa>  +  V0a) 

where  6  is  a  vector  to  be  found.     Changing  &>  into  yjrco 

W5T&)  =  ^(JWyjrco  +  VOyJrco. 

*  By  putting  d\f/  =  uS(  )w'  +  u>'S(  )w  in  this  equation,  equation  (11)  can  be 
deduced  but  as  this  method  has  already  been  applied  in  §  4  I  give  the  one  in 
the  text  to  shew  the  variety  of  Quaternion  methods.  [Note  added,  1892.  If  we 
use  the  theorem  in  the  foot-note  of  §  4,  equation  (11)  follows  at  once.] 
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Now  ct  being  self-conjugate  V^^  —  0.     Hence 

=  0StyZ  +  yjrd, 


whence 


or 

by  equation  (6a)  §  3  above.     Thus  finally 

murto  =  ^(jwyfrco  +  Vdyjrto 
or  mtato  =  ^(jwtyto  —  Vyjrto  (yjr  +  S^^^^V^jr^div^ 


.(13). 


This  completely  solves  the  problem  of  expressing  stress  in 
terms  of  strain  in  the  most  general  case. 

15.  ztco+  Ve'to,  where,  as  we  saw  in  last  section,  e'  is  perfectly 
arbitrary  so  far  as  the  strain  is  concerned,  is  the  force  on  the 
strained  area  to  due  to  the  pure  strain  ty.     And  again 

tj)to  +  Ve'to  or  q  (ixq^toq)  q'1  +  Ve'to 

is  that  due  to  the  strain  q-^r  (  )  q'1  or  y. 

To  find  the  force  on  an  area  which  before  strain  was  to0  let  its 
strained  value  after  -\|r  has  taken  place  be  to.  Then  by  equation 
(4),  §  145  of  Tait's  Quaternions,  mto0  =  tyco.     Hence 

Required  force  =  tsto  +  Veto  =  ^,(JwtoQ  +  V6to0  +  mVe'yjr^to^ 

by  equation  (13)  of  last  section.  If  the  rotation  now  take  place 
this  force  rotates  with  it  so  that  the  force  on  the  area  which  was 
originally  to0  is  after  the  strain  %  or  qyjr  ( )  q'1 

™„  =  q*QWto0q-1  +  qVdtoQ .  q~l  +  V(eq^~1to0 .  a-1) (14) 

where  e  =  mqe'q'1  and 

e  =  ^  the  stress-couple  per  mi.  vol.  of  unstrained  body. 

This  force  then  is  a  linear  vector  function  of  to0,  but  in  general 
even  when  e  =  0  it  is  not  self-conjugate.  When  both  e  and  the 
rotation  are  zero  we  see  that  the  rotation  vector  of  t  is  6  given  by 
equation  (12)  of  last  section. 

The  stress-force  can  be  shewn  as  in  §  13  to  be  tA  per  unit 
volume  of  the  unstrained  body.  Thus  since  the  corresponding 
stress-couple  is  2e  the  moment  exerted  by  the  stresses  on  any 
portion  of  the  body  round  an  arbitrary  origin  is 

jJSV(p+V)ryidS0  +  2jJJed%0. 
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But  this  moment  may  also  be  put  in  the  form 

JJV(p+V)rdX0, 

/FW  +  ff^  +  ,)  r,V A  +  J//F¥V A,   • 
by  equation  (9)  §  6  above.     Comparing  these  results 

2e=VP>rV1  J <14a>' 

in    the    notation    of   next   section.     This    equation    may    also    be 
deduced  from  equation  (15  I)  below  but  not  so  naturally  as  above. 


The  equations  of  equilibrium. 

16.  We  require  equation  (9)  §  14  above  to  prove  the  state- 
ment that  no  space-variations  of  -ty  or  q  are  involved  in  w.  It  is 
also  required  to  shew  that  the  fact  that  q  also  is  not  involved 
in  w  is  a  mathematical  sequence  of  the  assumption  that  the 
potential  energy  of  a  solid  is  the  sum  of  the  potential  energies  of 
its  elements.  Assuming  these  facts  however  we  can  arrive  at  the 
equation  of  stress  (11)  §  14  in  a  different  way  from  the  above. 
We  shall  also  obtain  quite  different  expressions  for  ot,  t,  &c. 
and  most  important  of  all  we  shall  obtain  the  equations  of 
equilibrium  by  obtaining  t  explicitly  in  terms  of  the  displacement 
and  its  space  derivatives.     From  §  12  above  we  have 

X°>  =  ^w?_1>  Xw  =  ^ (fla>q)> 

x>x  =  +*  =  y (15), 

(say)  as  in  Tail's  Quaternions,  §  365.  Thus  Way  =  VficoVfip^pJ 
where  p  is  put  as  it  will  be  throughout  this  section  for  p  +  rj,  the 
vector  coordinate  after  displacement  of  the  point  p.  From  this  as 
we  have  seen  in  the  Introduction  we  deduce  that  the  coordinates 
of  ^  are  the  A,  B,  C,  a,  b,  c  of  Thomson  and  Tait's  Nat.  Phil. 
App.  C. 

We   may  as  do  those  authors  regard  w  as  a  function  of  M^. 
Thus  :— 

Stu  =  —  SSxV£*(jw£, 

=  -  s$x'xZ*<iwS-  sxBx&aw£ 

By  equation  (4a)  §  3  above,  each  of  the  terms  in  this  last 
expression 

=  -  £Hx*awVi» 
8iv  —  —  2S8rj1^(ju<Vl    (15a), 
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Comparing  this  with  equation  (6)  §  14  and  putting  in  both 
8i]1  =  co'Scop'  where  &>',  co  are  arbitrary  constant  vectors,  we  get 

mSco'cpco  =  ^Sco'^qiu^co. 

Hence,  since  co'  is  quite  arbitrary 

m<j>  =  2x*awX  (156). 

From  equation  (10)  §  14  above  which  defines  vr  we  see  that 

m-as  =  2i/r4,QWiJr  (15c), 

which  we  jshould  expect  since  we  have  already  seen  that  ct  is  the 
value  of  <f>  when  there  is  no  rotation  and  therefore  %  =  ■%'  =  ty. 
Now  since 

we  deduce  by  any  one  of  the  processes  already  exemplified  that 

where  of  course  the  differentiations  of  *q  must  not  refer  to  yjr. 
We  see  then  from  equation  (15c)  that 

which  is  equation  (11)  §  14. 

Our  present  purpose  however  is  to  find  the  equations  of 
equilibrium.  Let  co0  be  the  vector  area  which  by  the  strain  % 
becomes  co.     Thus*  as  in  last  section 

mto0=  x'o  (15e). 

Further  let  2e  be  the  stress-couple  per  unit  volume  of  the 
unstrained  solid  so  that  2e/m  is  the  same  of  the  strained  solid. 
As  we  know,  e  is  quite  independent  of  the  strain.  By  §  13 
we  see  that  the  force  rco0  on  the  area  which  before  strain  was  co0  is 
<bco  +  Veco/m.     Therefore 

T0\>  =  2X*dw&>o  +  V*X~1q)o  (15/). 

We  saw  in  last  section  that  the  force  per  unit  volume  of  the  un- 
strained solid  is  tA  and  the  couple  2e.     Hence 

F  +  tA  =  0 (log), 

M  +  2e=0 (lo/i), 

*  See  §  83  below. 
M.  3 
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are  the  equations  of  equilibrium,  where  F  and  M  are  the  external 
force  and  couple  per  unit  volume  of  the  unstrained  solid.  All  that 
remains  to  be  done  then  is  to  express  t  in  terms  of  e  and  the 
displacement.     Putting 

P  +  V  =  P      (15*')> 

as  already  mentioned,  we  have 

X<D  =  -p1'Sa>V1   (15  j), 

so  that  by  equation  (6//)  §  8a  above  we  have 

x'-'o,  =  -  3  vPl'p;s»vxvjss7ytvj3PMP; (is*). 

Therefore  b}^  equation  (15/) 

to,  =  -  2P;SVMiua>-3VeVP;p2'S(oVlVJSV1V2VsSp1'p:p:;...{15l). 

It  is  unnecessary  to  write  down  what  equation  (15#)  becomes 
when  we  substitute  for  t,  changing  e  into  —  M/2  by  equation 
(loh).  In  the  important  case  however  when  M  =  0,  the  equation 
is  quite  simple,  viz. 

T=2Pl'SV1*a_wA   (15m). 

Addition  to  §  10,  Dec,  1887  (sent  in  with  the  Essay).  [The 
following  considerations  occurred  just  before  I  was  obliged  to  send 
the  essay  in,  so  that  though  I  thought  them  worth  giving  I  had 
not  time  to  incorporate  them  in  the  text. 

It  is  interesting  to  consider  the  case  of  an  isotropic  body. 
Here  w  is  a  function  of  the  three  principal  elongations  only  and 
therefore  we  may  in  accordance  with  §§  14  and  15  suppose  it  a 
function  of  a,  b,  c  or  in  accordance  with  §  16  of  A,  B,  G  where 

a  =  -StyZ)    b  =  -Sty*y2,    C  =  -£ftr3t/3 (A). 

A  =  -££¥£  B  =  -SWZI%  G  =  -S^^/S (B). 

Let  us  use  x,  y,  z,  X,  Y,  Z  for  the  differential  coefficients  of  w 
with  respect  to  a,  b,  c,  A,  B,  G  respectively.     Thus 

dw  =  xda  ■+  ydb  4-  zdc  —  —  Sd^  (x  4-  yyjr  +  zyj/2)  £ 

as  can  be  easily  proved  by  means  of  equation  (6a)  §  3  above.    But 

dtu  =  —  SdyJr^Qiu^. 
^<JW  =  x  +  yyjr  +  z^r1  (G), 

by  §  4  above.     Similarly  we  have 

*qw  =  x+  m+zv* (D). 
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(Notice  in  passing  that  to  pass  to  small  strains  is  quite  easy 
for  ^qw  is  linear  and  homogeneous  in  -i/r  —  1  so  that  ^qw  =  x  +  yyjr 
where  y  is  a  constant  and  x  +  y  a  multiple  of  a  —  3.)* 

From  equation  (G)  we  see  that  6  in  equation  (12)  §  14  is  zero 
and  therefore  that  from  equation  (13) 

mts  =  ^(JWi/r  =  X^r  +  y\jr2  +  Z\jr3 (E). 

Again  from  equation  (14)  §  15 

too  =  xqcorf1  +  yXa>  +  zy^rw  +  VeX~lw    (F). 

Similarly  from  equations  (156)  and  (15/)  §  16  we  may  prove  by 
equation  (D)  that 

^4  =  xxx'+Y(xxy+z(xxy (G), 

and  r  =  2XX  +  2  Yxx'x  +  2ZXXXX'X  +  VeX'  (  ) (H). 

If  we  wish  to  neglect  all  small  quantities  above  a  certain  order 
the  present  equations  pave  the  way  for  suitably  treating  the 
subject.  I  do  not  however  propose  to  consider  the  problem  here 
as  I  have  not  considered  it  sufficiently  to  do  it  justice.] 

Variation  of  Temperature. 

16a.  The  iu  which  appears  in  the  above  sections  is  the  same  as 
the  w  which  occurs  in  Tait's  Thermo-dynamics,  §  209,  and  therefore 
all  the  above  work  is  true  whether  the  solid  experience  change  of 
temperature  or  not.  w  will  be  a  function  then  of  the  temperature 
as  well  as  i/r.  To  express  the  complete  mathematical  problem  of 
the  physical  behaviour  of  a  solid  we  ought  of  course  instead  of  the 
above  equations  of  equilibrium,  to  have  corresponding  equations  of 
motion,  viz.  equations  (15A),  (15£)  and  (instead  of  (I5g)) 

Dp'=F  +  rA    (15n), 

where  D  is  the  original  density  of  the  solid  at  the  point  con- 
sidered. Further  we  ought  to  put  down  the  equations  of  con- 
duction of  heat  and  lastly  equations  (16/)  and  (16fZ)  below. 

We  do  not  propose  to  consider  the  conduction  of  heat,  but 
it  will  be  well  to  shew  how  the  thermo-dynamics  of  the  present 
question  are  treated  by  Quaternions. 

*  [Note  added,  1892.  In  the  original  essay  there  was  a  slip  here  which  I  have 
corrected.  It  was  caused  by  assuming  that  \p  instead  of  \p  -  1  was  small  for  small 
strains.     In  the  original  I  said  "  where  y  is  a  constant  and  x  is  a  multiple  of  a."] 

3—2 
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Let  t  be  the  temperature  of  the  element  which  was  originally 
d%0  and  Ed%0  its  intrinsic  energy.     Let 

where  M,  JV  are  a  linear  self-conjugate  function  and  a  scalar 
respectively,  both  functions  of  ^r  and  t.  Here  BHd%0  is  the  heat 
required  to  be  put  into  the  element  to  raise  its  temperature  by  Bt 
and  its  pure  strain  by  Byjr.  Now  when  t  is  constant  BH  must  be  a 
perfect  differential  so  that  we  may  put 

where /is  some  function  of -^  and  t.     Thus 

m=-tSty^<jfZ+N8t    (16). 

Now  we  have  seen  in  §  14  that  the  work  done  on  the  element 
during  the  increment  By{r,  divided  by  d%0 

=  —  mSh'ty-ty''1^^. 

Thus  by  the  first  law  of  Thermo-dynamics 

BE  =  J8H  -  mStyf-'frZ 

=  JNBt  -  JtSSyfr&afZ  -  mSBf^1  fa£, 

where  J  is  Joule's  mechanical  equivalent.     Thus 

JN=dE/dt (16a), 

and  |m(wf1  +  f1w)=^w; (166), 

where  w  =  E  -  Jtf (16c). 

To  apply  the  second  law  we  go  through  exactly  the  same  cycle 
as  does  Tait  in  his  Thermo-dynamics,  §  209,  viz. 

(yjr,  t)  (yjr  +  Byfr,  t)  (f  +  8f,t  +  St)  (f,  t  +  Bt)  ty,  t). 

We  thus  get* 

or  Jf=-dw\dt (16d), 

the  arbitrary  function  of  t  being  neglected  as  not  affecting  any 
physical  phenomenon.     Substituting  for  w  from  equation  (16c), 

Jtdf/dt  =  dEjdt  =  JN (We) 

by  equation  (16a).     Thus  from  equation  (16) 

BH=tBf (16/), 

[By  assuming  from  the  second  law  that  the  work  done  by  the  element  in  the 
cycle,  i.e.  the  sum  of  the  works  done  by  it  during  the  first  and  third  steps  is  J5t/t 
multiplied  by  the  heat  absorbed  by  the  element  in  the  third  step.    Note  added,  1893.] 
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so  that  in  elastic  solids  as  in  gases  we  have  a  convenient  function 
which  is  called  the  "entropy".  Thus  the  intrinsic  energy  Ed%, 
the  entropy  fd%3  and  the  stress  w  have  all  been  determined  in 
terms  of  one  function  w  of  \fr  and  t  which  function  is  therefore  in 
this  general  mathematical  theory  supposed  to  be  known. 

If  instead  of  regarding  w  (which  with  the  generalised  meaning 
it  now  bears  may  still  conveniently  be  called  the  "potential  energy" 
per  unit  volume)  as  the  fundamental  function  of  the  substance  we 
regard  the  intrinsic  energy  or  the  entropy  as  such  it  will  be  seen 
that  one  other  function  of  -v/r  must  also  be  known.  For  suppose  / 
the  entropy  be  regarded  as  known.     Then  since  dwldt  =  —  Jf 

w=  W-Jjfdt (16*7), 

where  the  integral  is  any  particular  one  and  W  is  a  function  of  yjr 
only,  supposed  known.     Again 

E  =  w  +  Jtf 

or  E=W+J(tf-ffdt).. (16/0- 

Thus  all  the  functions  are  given  in  terms  of/'  and  W.     Similarly 
if  E  be  taken  as  the  fundamental  function 

w  =  t(W'-fEdt/f) (160, 

Jf=E/t  +  fEdt/t*-W (16/), 

where  as  before   the  integral  is  some  particular  one  and  W  is  a 
function  of  i/r  only. 

Small  strains. 

17.  We  now  make  the  usual  assumption  that  the  strains  are 
so  small  that  their  coordinates  can  be  neglected  in  comparison 
with  ordinary  quantities  such  as  the  coefficients  of  the  linear 
vector  function  ^(Jio.  We  can  deduce  this  case  from  the  above 
more  general  results. 

To  the  order  considered  q  =  1  so  that  by  equation  (10)  §  14  above 

(f)  =  ts. 

We  shall   use -the  symbol  ot  rather  than  cf>  for  the  same  reasons 
explained  in  §  13  above  as  induce  us  to  use  ^  rather  than  %. 

Remembering  that  (o0  and  &>  are  now  the  same  and  that  ty  may 
be  put  =  1  so  that  V^^(j_w^  =  V^(jw^  =  0  and  therefore  6  of 
equation  (12)  =  0,  we  have  from  equation  (13), 

■UTO)  =  tyQIUOJ (17). 


o 
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Of  course  we  do  not  require  to  go  through  the  somewhat  com- 
plicated process  of  §§  14,  ^  5  to  arrive  at  this  result.  In  fact  in 
equation  (6)  §  14,  we  may  put  m  =  1  and  p+r)V  =V  so  that 

hw  —  —  S8f]1(j}^1  =  —  SBrjjVrV^ 

and  therefore  by  §  3  above 

8iu  =  -  S8xfaS=  -  SSffaZ. 

But  by  §  5  above 

Bw  =  —  StyZfQwZ, 

and  therefore  by  §  4  we  get  equation  (17). 

18.  It  is  convenient  here  to  slightly  change  the  notation. 
For  yjr,  x  we  shall  now  substitute  ^  + 1,  %+  1  respectively.  This 
leads  to  no  confusion  as  will  be  seen. 

With  this  notation  the  strain  being  small  the  stress  is  linear 
in  i/r  i.e.  ^qw  is  linear  and  therefore  w  quadratic.  Now  for  any 
such  quadratic  function 

iu  =  -S^„l(jw^2 (18), 

for  we  have  by  §  5, 

dw  =  —  SdyjrZ  ^(jw£. 

Put  now  \|r  =  n-^r'.     Then  because  ^qw  is  linear  in  \|r 

^qiu  =  n  ^(Jto, 

where  ^(j'w  is  put  for  the  value  of  ^(jw  when  the  coordinates  of 

•x/r'  are  substituted  for  those  of  i/r.     Thus  keeping  ty'  constant  and 

varying  n, 

dw  =  —  ndnS-^r'^d'w^, 

whence  integrating  from  n  =  0  to  w  =  1  and  changing  -vj/  into  ty 
we  get  equation  (18). 

From  equation  (18)  we  see  that  for  small  strains  we  have 

W  =  -Syjr^/2 (19). 

Now  w  is  quadratic  in  \Jr  and  therefore  also  quadratic  in  ■or,  so 
that  regarding  w  as  a  function  of  ct  we  have  as  in  equation  (18) 

w  =  -  Svr£  m(jwy2, 

so  that  by  equation  (19)  and  §  4  above 

*  =  WQW (20). 

All  these  results  for  small  strains  are  well-known  in  their  Car- 
tesian  form,  but  it   cannot  be  bias  that  makes  these  quaternion 
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proofs  appear  so  much   more  natural  and  therefore  more  simple 
and  beautiful  than  the  ordinary  ones. 

19.     Let  us  now  consider  (as   in  §  1G  is  really  done)  w  as  a 

function  of  the  displacement.     Now  w  is  quadratic  in  ty,  and  ^  is 

linear  and  symmetrical  in  Vx  and  rjt.     In  fact  from  equation  (3) 

§  12  above,  remembering  that  the  %  of  that  equation  is  our  present 

t|t  +  1,  we  have 

2-v/ro)  =  —  ?;1»S'&)V1  —  V1S(i)7)1. 

Therefore  we  may  put 

iv  =  iv(Vl,  Vt,  t]2,  V2) (21), 

where  w  (a,  ft,  y,  8)  is  linear  in  each  of  its  constituents,  is  sym- 
metrical in  a  and  ft,  and  again  in  y  and  8,  and  is  also  such  that 
the  pair  a,  ft  and  the  pair  y,  8  can  be  interchanged.  [This  last 
statement  can  be  made  true  if  not  so  at  first,  by  substituting  for 
w  (a,  ft,  y,  8),  w  (a,  ft,  y,  8)/2  +  w  {y,  8,  a,  ft)/2,  as  this  does  not  affect 
equation  (21).]  Such  a  function  can  be  proved  to  involve  21  inde- 
pendent scalars,  which  is  the  number  also  required  to  determine 
an  arbitrary  quadratic  function  of  -yjr,  since  -v/r  involves  six  scalars. 

Thus  we  have  the  two  following  expressions  for  8iu,  which,  we 

equate 

-S8f^=w{8VvV1>v.2,V2)  +  iu{Vl,Vi,8V2,V2) 

=  2w(8ri1,V1,Vs,Vs), 

or*  by  §  3  above, 

-S8Vl^l  =  2w(8i]uV1,V2,V2). 

Now  let  us  put  8r)  =  co'Scop  where  co'  and  co  are  arbitrary  constant 
vectors.     We  thus  get 

Sco'vTCO  =  —  2lU  (O)',  CO,  7]v  Vt)  =  2Sco'%W  (£,  CO,  7]^  Vj). 

Whence  since  co'  is  quite  arbitrary, 

™  =  2^(£&>,7?1,V1) (22). 

The  statical  problem  can  now  be  easily  expressed.     As  we  saw 
in  §  14.  equations  (7)  and  (8),  it  is  simply 

F  +  FVM/2  +  otA  =  0 (22a), 

*  [Note  added,  1892.    Better  thus:— by  §  3  above, 

therefore  by  §  4,  z?«  =  2fio  (j",  w,  yv  Vi) 

for  i'w  (£,  w,  t)v  Vi)  regarded  as  a  function  of  w  is  clearly  self  conjugate. 
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throughout  the  mass;  and  at  the  surface 

Fs-  VUdt7HL/2-<srUd2=0 (226), 

where  F,  M  are  the  given  external  force  and  couple  per  unit 
volume  and  F,,  is  the  given  external  surface  traction  per  unit 
surface.     Substituting  for  w  from  equation  (22) 

F  +  FVM/2  +  2£w  (£,  A,  Vl ,  Vt)  =  0 
F,  -  VUdZNLfi  -  2£w (£  UdX,  Vl,  ^)  =  0[ 


.(23). 


Isotropic  Bodies. 

20.  The  simplest  way  to  treat  these  bodies  is  to  consider  the 
(linear)  relations  between  zr  and  i/r. 

In  the  first  place  notice  that  i/r  can  always  be  decomposed  into 
three  real  elongations  (contractions  being  of  course  considered  as 
negative  elongations).  Thus  i  being  the  unit  vector  in  the  direc- 
tion of  such  an  elongation, 

"tyco  =  —  %eiSioi. 

The  elongation  —  eiSicu  will  cause  a  stress  symmetrical  about  the 
vector  i,  i.e.  a  tension  Ae  in  the  direction  of  i  and  a  pressure  Be  in 
all  directions  at  right  angles  ;  A  and  B  being  constants  (on  account 
of  the  linear  relation  between  ot  and  yjr)  independent  of  the  direc- 
tion of  i  (on  account  of  the  isotropy  of  the  solid).  This  stress 
may  otherwise  be  described  as  a  tension  ( A  +  B)  e  in  the  direction 
of  i  and  a  hydrostatic  pressure  Be.     Thus 

OT&)  —  —  (A  +  B)  XeiSiro  —  BcoXe 
=  (A+B)fco  +  BcoStyt;. 

To  obtain  the  values  of  A  and  B  in  terms  of  Thomson  and 
Tait's  coefficients  k  and  n  of  cubical  expansion  and  rigidity  respec- 
tively ;  first  put 

yjrco  —  eo)  and  •era)  =  Skeco, 
and  then  put 

yfrco  =  V\cofi  and  ot&)  =  2nV\u>ix, 

X  and  /x  being  any  two  vectors  perpendicular  to  each  other.     We 
thus  get 

A  +  B  =  2n,  B=-(k-  2n/3), 

whence 

wo)  =  2nfco  -  (k  -  2n/3)  w»S'£tyf   (24). 
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From  this  we  have 

yfra)  =  *ra>/2ib  +  coStyt  (k  -  2nfd)/2n, 
but  from  the  same  equation 

■'■   ^v'+iik-tk)*8™ (25>- 

Equation  (24)  gives  stress  in  terms  of  strain  and  (25)  the  converse. 


21.  We  can  now  give  the  various  useful  forms  of  w  for 
isotropic  bodies  for  from  equation  (19)  §  18, 

W  =  -  SyJr^/2. 

Therefore  from  equations  (24)  and  (25)  respectively 

w  =  -  n  mf  +  $  (k  -  2n/8)  flft^Sftf  &    {26), 

w=~i  w  -  ${l  -  ik)  *.<*«6 w)- 

Therefore  again  from  §  3  above  and  from  equation  (26), 

w  =  -  nSVl^Vl  +%(k-  2n/3)  (SVV)2 
or 

2w (Vl,  Vlt  v,,  V2)  =  nSVinJWt%  +  nSWJSw, 

+  (k  -  2n/3)  SVlVlSVa%  (28). 

Hence  from  equation  (22) 

-osw  =  —  nScoV .  ?;  —  nVfirj^  —  (m  —  ?i)  w$  v7^  ( 29), 

where  m  is  put  for  k  +  n/S.     This  last  could  have  been  deduced  at 
once  from  equation  (24)  by  substituting  for  -v/ro). 

Thus  the  equations  (23)  for  the  statical  problem  are 

F  +  V VM/2  =  nV2,;  +  mVS  Vy  (30), 

-  Fs  +  VUdXM/2  =  iiSUdXV  .  v  +  nVfirj^dl 

+  (m-ri)UdXSVr) (31). 

We  now  proceed  to  apply  these  results  for  small  strains  in 
isotropic  bodies  to  particular  cases.  These  particular  cases  have 
all  been  worked  out  by  the  aid  of  Cartesian  Geometry  and  they 
are  given  to  illustrate  the  truth  of  the  assertion  made  in  the 
Introduction  that  the  consideration  of  general  problems  is  made 
simpler  by  the  use  of  Quaternions  instead  of  the  ordinary  methods. 
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Particular  integral  of  equation  (30)*. 

22.     Since  from  equation  (30)  (F  being  put  for  simplicity  in- 
stead of  F  +  FVM/2)  we  have 

nV2v  =  F-  mVSVij 

we  obtain  as  a  particular  case  by  equations  (18)  and  (19)  §  10, 

47T/M7  =fffu  (F  -  mrVSVij)  d%, 

where  u  has  the  meaning  explained  in  §  9,  and  the  volume  inte- 
gral extends  over  any  portion  (say  the  whole)  of  the  body  we  may 
choose  to  consider.  To  express  ///ttV S  v7 i]d%  as  a  function  of  F  put 
in  this  term  u  =  —  \Up^l  where  p  is  taken  for  the  p  —  p  of  §  9, 
and  apply  equation  (9)  §  6.     Thus 

*7mv  =  fffuFd%  +  ±mfffUPlVyS?vd$ 

=  fffuFd%  -  %mfffUpV*SVVd%  +  a  surf.  int. 

=  JJJuFd$  -  .y— ^ — -JJJUpSVFdZ  +  the  surf.  int. 

for  by  equation  (30)  #vF  =  (m  +  n) V2SVt].  Now  (in  order  to 
get  rid  of  any  infinite  terms  due  to  any  discontinuity  in  F)  apply 
equation  (9)  §  6  to  the  second  volume  integral.     Thus 

47^  =  JJJtiFd%  +  -j-    -?- JffSF \7.Upd%  +  a  surf.  int. 

The  surface  integral  may  be  neglected  as  we  may  thus  verify. 
Call  the  volume  integral  <kiriirj '.     Thus 

4™vy  =  4ttF  +  2  ^L—fffSFV .  VUpto 

4,TrmVSVV'  =  -  JffSFV  .  Vadg  +  =-^-  -.  fffSFV  .  VSV  Upd%, 
n  Zn{m  +  it)  r 

so  that  putting  V  Up  =  —  2u  we  get 

/tVy  +  mV8Vv'  =  F, 
whence  we  have  as  a  particular  solution  of  equation  (30)  rj  =  r/  or 

'  =  SSJD4**  +  s™^)///'^  •  Dp*     (32). 

*  [Note  added,  1892.     For  a  neater  quaternion  treatment  of  this  problem  see 
Phil.  Mug.  June,  1892,  p.  493.] 
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This  is  generally  regarded  as  a  solution  of  the  statical  problem 
for  an  infinite  isotropic  body.  In  this  case  some  law  of  con- 
vergence must  apply  to  F  to  make  these  integrals  convergent. 
Thomson  and  Tait  {Nat.  Phil.  §  730)  say  that  this  law  is  that  Fr 
converges  to  zero  at  infinity.  This  I  think  can  be  disproved  by 
an  example.  Put  *  F?*  =  r~a\  where  X  is  a  constant  vector  and  a  a 
positive  constant  less  than  unity.  Equation  (32)  then  gives  for  the 
displacement  at  the  origin  due  to  the  part  of  the  integral  extend- 
ing throughout  a  sphere  whose  centre  is  the  origin  and  radius  R 

m  +  Sn     Rl~a  „ 
3?i  (m  +  n)  1  -  a 

Putting  R=  oo  ,  7]  also  becomes  x  .  The  real  law  of  convergence 
does  not  seem  to  me  to  be  worth  seeking  as  the  practical  utility 
of  equation  (32)  is  owing  to  the  fact  that  it  is  a  particular  integral. 

The  present  solution  of  the  problem  has  only  to  be  compared 
with  the  one  in  Thomson  and  Tait's  Nat.  Phil.  §§  730-1  to  see  the 
immense  advantage  to  be  derived  from  Quaternions. 

It  is  easy  to  put  our  result  in  the  form  given  by  them.  We 
have  merely  to  express  SFV .  Up  in  terms  of  F  and  r2$FV .Vu 
where  r  is  put  for  the  reciprocal  of  u.     Noting  that 

Vm  =  - u3p,     Up  =  up,     SFV  .p  =  -T, 

we  have  at  once 

SFV.Up  =  -uF-u3pSFp 

SFV  .  Vu  =  w3F  +  Su'pSFp 
therefore  eliminating  pSFp, 

SFV  .Up  =  -  uF2J'3  -  r2SFV  .  Vt«/3, 

.-.    v  =  {247rn  (m  +n)}*7//cfe  {2(2m+Sn)  wF-mr^FV.Vw}...(33), 
which  is  the  required  form. 

23.     Calling  the  particular  solution  77'  as  before  and  putting 

rj  =  7]  +7] 
the  statical  problem  is  reduced  to  finding  r/"  to  satisfy 

nVV  +  mV>Sr\7V'  =  0 

*  [Note  added,  1892.  This  is  not  legitimate  since  it  makes  F=  oo  for  r  =  0.  The 
reasoning  is  rectified  in  the  Phil.  Mag.  paper  just  referred  to  by  putting  F  — 0  from 
r  —  0  to  r  —  b  and  Fr  =  r~"  X  from  r  —  b  to  ;•  —  s.  .  ] 
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and  the  surface  equation  either 

r(  +  7]"  =  given  value, 
i.  e.  rj"  =  given  value, 

or  zjlldlt  =  given  surface  traction, 

i.e.  by  equation  (29)  §  21  above, 
nSUdXV .  v"  +  nVfiri"  UdX  +  (m  -  n)  Ud&SVv"  =  known  value. 
This  general  problem  for  the  spherical  shell,  the  only  case 
hitherto  solved,  I  do  not  propose  to  work  out  by  Quaternions,  as 
the  methods  adopted  are  the  same  as  those  used  by  Thomson  and 
Tait  in  the  same  problem.  But  though  each  step  of  the  Cartesian 
proof  would  be  represented  in  the  Quaternion,  the  saving  in 
mental  labour  which  is  effected  by  using  the  peculiarly  happy 
notation  of  Quaternions  can  only  be  appreciated  by  him  who  has 
worked  the  whole  problem  in  both  notations.  The  only  remark 
necessary  to  make  is  that  we  may  just  as  easily  use  vector,  surface 
or  solid,  harmonics  or  indeed  quaternion  harmonics  as  ordinary 
scalar  harmonics. 

Orthogonal   coordinates. 

24.  It  is  usual  to  find  what  equations  (22a)  of  §  19  and  (3)  of 
§  12  become  when  expressed  in  terms  of  any  orthogonal  coordi- 
nates. This  can  be  done  much  more  easily  by  Quaternions  than 
Cartesian  Geometry.  Compare  the  following  investigation  with 
the  corresponding  one  in  Ibbetson's  Math.  Theory  of  Elasticity, 
Chap.  V. 

Let  x,  y,  z  be  any  orthogonal  coordinates,  i.e.  let  x  =  const, 
y  =  const.,  z  —  const.,  represent  three  families  of  surfaces  cutting 
everywhere  at  right  angles.  Particular  cases  are  of  course  the 
ordinary  Cartesian  coordinates,  the  spherical  coordinates  r,  6,  <f> 
and  the  cylindrical  coordinates  r,  (f>,  z.  Let  Dx,  Dv,  Dz  stand  for 
differentiations  per  unit  length  perpendicular  to  the  three  coordi- 
nate surfaces  and  let  X,  fit  v  be  the  unit  vectors  in  the  correspond- 
ing directions.     Thus 

V  =  \DX  +  fxDv  +  vD,. 
Thus,  using  the  same  system  of   suffixes  for  the  D's  as  was  ex- 
plained in  connection  with  V  in  §  1, 

<£A  =  D^X  +  Dn4>ji  +  AA* (34), 


or 


4>A  =  Dx  (</>\)  +  Du  (M  +  A  (<M  -  0  (Dx\  +  Dulx  +  Dzv). .  .(35). 
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25.  Now  to  put  equation  (22a)  §  19  into  the  present  coordi- 
nates all  that  is  required  is  to  express  otA  in  terms  of  those  co- 
ordinates. Let  the  coordinates  of  ot  be  PQRSTU.  Thus  from 
equation  (35)  we  have 

^A  =  DX{P\  +  Ufi  +  Tv)  +  Dy(U\  +  Qfi+  Sv)  +  DZ(T\  +  S/x  +  Rv) 

-  Z7  (Dx\  +  DyfJL  +  Dzl>). 

The  first  thing  then  is  to  find  Dx\,  Dxjjl,  Dy\  &c. 

Let  p2,  p3  be  the  principal  curvatures  normal  to  x  =  const.,  i.e. 
(by  a  well-known  property  of  orthogonal  surfaces)  the  curvatures 
along  the  lines  of  intersection  of  x  =  const.,  with  z  =  const.,  and 
y  =  const.  p2,  p3  will  be  considered  positive  when  *the  positive 
value  of  dx  is  on  the  convex  side  of  the  corresponding  curvatures. 
Similarly  for  q3qj\i\.     Thus  for  the  coordinates  r,  6,  <f> ; 

V*  =Pz  =  l!r>  ?a  =  cot  elr>  ?!  =  ri  =  r2  =  °- 
Again  for  r,  </>,  z ;  p.2  =  1/r  and  the  rest  are  each  zero. 

With  these  definitions  we  see  geometrically  that 

DXX  =  —  /j.ql  —  vi\,  Dx/m  =  \qt ,  Dxv  =  \i\ (*>6). 

Similarly  for  Dy\  D!ffi,  Dz\  &c.     Thus 

v&=\(DxP  +  DyU+DzT)  \-fj,()  +  v() 

+  \(-Qp2-RP3  +  rr1+l\)  +  v()  +  v() 

+  *>{(P>+Ps)P  +  (qB  +  ql)U  +  (r1  +  rJT}  +  M,{}  +  v{}, 
or 

m A  =  \  [DXP  +  DyU  +  DzT+P(p2  +  p3) -  Qp2 - Rps 

+  T(2r1  +  rz)+U(q3  +  2qt)}+H,{}  +  v{} (37). 

26.  The  other  chiefly  useful  thiug  in  transformation  of  coor- 
dinates in  the  present  subject  is  the  expression  for  the  strain 
function  yjr  in  terms  of  the  coordinates  of  displacement.  Let 
u,  v,  w  be  these  coordinates.  Now  by  equation  (3)  §  12,  remem- 
bering (§  18)  that  -»/r  =  ^  —  1  we  have 

whence         —  2-v/r\  =  \S\Dxr)  +  /j,S\.Dy7j  +  vS\Dzr}  —  Dxrj. 
But     Dxt)  =  \Dxii  +  fiDxv  +  vDxw  +  uDx\  -f-  vDxfi  4-  wDxv 

=  X  (Dxii  +  vq1  +  wrt)  +  fM  ( Dxv  -  uqx)  +  v  (Dxio  -  urt). 

*  [Note  added,  1892.     This  is  contrary  to  the  usual  convention.] 


y (39). 
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Similarly 

Dyt)  =  \  (Dyu  -  vp2)  +  ii  (Dyv  ■+  wr2  +  ups)  +  v  (Dyw  -  vr2) 
Dzr)  =  \  {Dzu  -  wp3)  +  jx  (Dzv  -  wq3)  +  v  (Dzw  +  ups  +  vq3) 

.  • .     2t/tA,  =  2\  (Dxv  +  vqt  +  wr\)  +  fi  (Dyu  +  Dxv  -  uq1  -  vp2) 

+  v  (Dzii  +  Dxw  -  w\  -wp3). .  .(38). 

But  with  Thomson  and  Tait's  notation  for  pure  small  strain 

2\jr\  =  2\e  +  fie  +  vb, 
*e  =  Dxu  +vq1    +wr1  \ 

f  =  Dyv  +  wr2    +  up, 
g  =  Dziu  -f  up3    +  vq3 
a  =  Dyw  +  Dzv  —  wq3  —  vr2 
b  =  Dzu  +  Dxw  —  ui\  —  wp3 
c  =  Dxv  +  ByU  -  vp2  -  uql 

Thus  we  have  efgabc  in  terms  of  the  displacement  and  we 
have  already  in  equation  (24)  §  20,  which  expresses  ot  in  terms  of  i/r, 
found  the  values  of  PQRSTU  in  terms  of  e,  &c.  Finally  the 
expression  in  equation  (37)  for  otA  gives  us  the  equations  of  equi- 
librium in  terms  of  P,  &c.  Thus  we  have  all  the  materials  for 
considering  any  problem  with  the  coordinates  we  have  chosen. 

All  these  results  can  be  at  once  applied  to  spherical  and  cylin- 
drical coordinates,  but  as  this  has  nothing  to  do  with  our  present 
purpose — the  exemplification  of  Quaternion  methods — we  leave 
the  matter  here. 

Let  us  as  an  example  of  particular  coordinates  to  which  this 
section  forms  a  suitable  introduction  consider  St  Venant's  Torsion 
Problem  by  means  of  cylindrical  coordinates. 


Saint-  Venant's  Torsion  Problem. 

27.  In  this  problem  we  consider  the  equilibrium  of  a  cylinder 
with  any  given  cross-section,  subjected  to  end-couples,  but  to  no 
bodily  forces  and  no  stress  on  the  curved  surface. 

*  [Note  added,  1892.  In  the  Phil.  Mag.  June,  1892,  p.  488,  there  is  a  mistake  in 
the  equation  just  preceding  equation  (31)  and  there  are  two  mistakes  in  equation  (31). 
In  the  first  of  these  I  (2Dtu  -  tx/v  -  /zzr  w)  should  he  2I(T)tv  -  vs,v  -  jus  ?r).  In 
equation  (31)  all  the  2's  should  be  dropped.] 
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We  shall  take  r,  <f),  z  as  our  coordinates,  the  axis  of  z  being 
parallel  to  the  generating  lines  of  the  cylinder.  Let  X,  fi,  v  be 
the  unit  vectors  in  the  directions  of  dr,  d(p,  dz  respectively  and  let 

7)  —  UX  +  VfJL  +  wv 

as  before. 

We  shall  follow  Thomson  and  Tait's  lines  of  proof — i.e.  we 
shall  first  find  the  effect  of  a  simple  torsion  and  then  add  another 
displacement  and  so  try  to  get  rid  of  stress  on  the  curved  surface. 

Holding  the  section  z  =  0  fixed  let  us  give  the  cylinder  a  small 
torsion  of  magnitude  t,  i.e.  let  us  put 

7]  =  TZl'fX (40), 

for  all  points  for  which  rz  is  small. 

The  practical  manipulation  of  such  expressions  as  this  is  almost 
always  facilitated  by  considering  the  general  value  of  Q(y&1,  tjJ 
where  Q  is  any  function  linear  in  each  of  its  constituents.  Thus 
in  the  present  case 

Q  (Vv  Vl)  =  r  {zQ  (X,  fj,)  -  zQ  (ft,  X)  +  rQ  (v,  fi)}. 

[If  Q  is  symmetrical  in  its  constituents,  e.g.  in  the  case  of  stress 
below  this  reduces  to  the  simple  form  Q(V1(  rj^  =  rrQ  (v,  /j,).] 
From  this  we  at  once  see  that  rj  satisfies  the  equation  of  internal 
equilibrium 

nV'rj  +  rnVSVy  =  0, 

for  putting  Q  (a,  /3)  —  a/3 

VV=.T  (2zv  -  rX)  =  tV  O2  -  r2/2) 

so  that  both  SVr]  and  V2??  =  0. 

Again  the  value  for  Q  at  once  gives  us  the  stress  for 

otg>  =  —  nScoV  .  7)  —  nVlSo}7)1  —  (m  —  n)  gjSVt}, 

or  ctg)  =  —  mr  (fjuScov  +  vScofi)  (41), 

which  is  a  shearing  stress  mr  on  the  interfaces  perpendicular  to 
fi  and  v. 

Putting  (o  =  the  unit  normal  of  the  curved  surface  we  have  for 

the  surface  traction 

vru)  =  —  nrrvSoi/u,. 

In  the  figure  let  the  plane  of  the  paper  be  z  =  0,  0  the  origin, 
P  a  point  on  the  curved  surface  and  OM  the  perpendicular  from 
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0  on  the  tangent  at  P.     Thus  -  rSatfi  =  OP  cos  OPM  =  PM,  PM 
being  reckoned   positive  or  negative  according   as    it    is  in   the 


positive  or  negative  direction  of  rotation  round  Oz.     Thus  we  see 
that  the  surface  traction  is  parallel  to  Oz  and  =  nrPM. 

Hence  in  the  case  of  a  circular  cylinder  a  torsion  round  the 
axis  satisfies  all  the  conditions  of  our  original  problem,  but  this  is 
true  in  no  other  case. 

The  surface  traction  at  any  point  on  the  plane  ends  necessary 
to  produce  this  strain  is  vjv=nrr/M  by  equation  (41)  so  that  its 
moment  round  the  origin  is  nrjjr'dA,  where  dA  is  an  element  of 
area  and  the  integral  extends  over  the  whole  cross-section. 

28.     Let  us  now  assume  a  further  displacement 

7]  =  wv (42), 

where  wisa  function  of  r,  </>  only,  and  let  us  try  to  determine  w 
so  that  there  is  still  internal  equilibrium  and  so  that  the  stress  on 
the  curved  surface  due  to  w  shall  neutralise  the  surface  traction 
already  considered. 

In  the  present  case 

Q(V1,i?1)  =  Q(V«;,  v). 

Thus  SVi)  =  0  (since  v  is  perpendicular  to  Vw)  and  therefore  the 
equation  of  internal  equilibrium  gives 

V*w  =  0. 

Again  zsco  =  —  nSwV  .  ??  —  nVl8corjv 

or  ct&)  =  —  n  (vSaVw  +  VivScov)     (43), 

a   shear  =  ?iTVw  on    the    interfaces  perpendicular  to  Vw  and   v. 
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Thus  putting  co  for  the  unit  normal  to  the  curved  surface,  the 
present  surface  traction  will  neutralize  the  former  if 

ScoViv  =  —  ttSw/j,  =  —  tScov  (\r) 
=  |  tSvcoV  (r2), 
i-e.  dw/dn  =  d(Tr2/2)/ds, 

where  d/dn  represents  differentiation  along  the  normal  outwards 
and  d/ds  differentiation  along  the  positive  direction  of  the  bound- 
ing curve. 

We  leave  the  problem  here  to  the  theory  of  complex  variables 
and  Fourier's  Theorem. 

Observe  however  that  the  surface  traction  at  any  point  on  the 
plane  end  =  ctv  =  nViv  by  equation  (43),  and  therefore  that  the 
total  couple  =  -  nffrSfiVwdA  =  nff(dw/d<f>)  dA.  This  leads  to  the 
usual  expression  for  torsional  rigidity. 


Wires. 

29.  In  the  following  general  treatment  of  Wires  some  of  the 
processes  are  merely  Thomson  and  Tait's  translated  into  their 
shorter  Quaternion  forms ;  others  are  quite  different.  The  two 
will  be  easily  distinguished  by  such  as  are  acquainted  with  Thom- 
son and  Tait's  Nat.  Phil. 

The  one  thing  to  be  specially  careful  about  is  the  notation  and 
its  exact  meaning.     This  meaning  we  give  at  the  outset. 

The  wires  we  consider  are  not  necessarily  naturally  straight 
but  we  assume  some  definite  straight  condition  of  the  wire  as  the 
"  geometrically  normal  "  condition. 

The  variable  in  terms  of  which  we  wish  to  express  everything 
is  s  the  distance  along  the  wire  from  some  definite  point  on  it. 

Any  element  of  the  wire,  since  it  is  only  slightly  strained,  may 
be  assumed  to  have  turned  as  a  rigid  body  from  its  geometrically 
normal  position.  This  rotation  is  expressed  as  usual  (Tait's 
Quaternions,  §  354)  by  the  quaternion  q ;  the  axis  of  q  being  the 
axis  of  rotation,  and  the  angle  of  q,  half  the  angle  turned  through. 

co  is  taken  so  that  the  rate  of  this  turning  per  unit  length  of 
the  wire  is  qtoq'1  so  that  co  is  the  rate  of  turning  per  unit  of  length 
when  the  whole  wire  is  moved  as  a  rigid  body  so  as  to  bring  the 
M.  4 


50  ELASTIC    SOLIDS.  [§  29. 

element  under  consideration  back  to  its  geometrically  normal 
position.  Of  course  a>  is  a  function  of  q  and  its  derivative  with 
reference  to  s.  This  function  we  shall  investigate  later.  The 
resolved  part  of  qeocf1  parallel  to  the  wire  is  the  vector  twist  and 
the  resolved  part  perpendicular  to  the  wire  is  in  the  direction  of 
the  binormal  and  equal  to  the  curvature.  In  fact  to  is  the  vector 
whose  coordinates  are  the  k,  \,  t  of  Thomson  and  Tait's  Nat.  Phil. 
§  593.  When  we  are  given  g  or  w  for  every  point  we  know  the 
strain  of  the  wire  completely.  wQ  is  defined  as  the  naturally 
normal  value  of  &>,  i.e.  the  value  of  w  when  the  wire  is  unstressed. 
As  usual  we  take  p  as  the  coordinate  vector  of  any  point  of  the 
wire,  p  like  the  rest  of  the  functions  being  considered  as  a  func- 
tion of  s.  We  shall  denote  (after  Tait,  Quaternions,  Chap,  ix.) 
differentiations  with  regard  to  s  by  dashes. 

We  now  come  to  the  dynamical  symbols.  P  and  M  are  the 
force  and  couple  respectively  exerted  across  any  normal  section  of 
the  wire  on  the  part  of  the  wire  which  is  on  the  negative  side  of 
the  section  by  the  part  on  the  positive  side. 

Finally  let  X,  L  be  the  external  force  and  couple  per  unit 
length  exerted  upon  the  wire. 

30.  When  the  wire  is  strained  in  any  way  let  us  impose  a 
small  additional  strain  represented  by  an  increment  Bco  in  on  and 
an  increment  Be  in  the  elongation  at  any  point.  Then  the  work 
done  on  the  element  ds  by  the  stress-force  =  —  BeSFp'ds  and  that 
done  by  the  stress-couple  =  -  SqBcoq^Mds.  If  (as  we  assume, 
though  the  assumption  is  not  justified  in  some  useful  applications 
of  the  general  theory  of  wires)  F  and  M  to  be  of  the  same  order 
of  magnitude  the  former  of  these  expressions  can  be  neglected  in 
comparison  with  the  latter  for  Be  is  a  quantity  small  compared 
with  S&>.  Now  the  work  done  on  the  element  by  the  stress  =  the 
increment  of  the  element's  potential  energy  =  Bwds  where  w  is 
some  function  of  the  strain.     Hence 

Biv  =  —  SBcoq^Mq. 

Thus  w  is  a  function  of  to  ouly  and 

Btu  —  —  SBcoJVtu, 

whence  we  see  that 

q~1Mq  =  0>Vw  (46). 

Notice  that  M  is  Thomson  and  Tait's  £  rj,  £  and  q^Mq  their 
KLM  {Nat.  Phil.  %  594,  614). 
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31.  Now  since  the  strain  is  small,  q~1-M.q  is  linear  in  terms  of 
the  strain  and  therefore  in  terms  of  &).  Hence  we  see  that  w 
is  quadratic  in  terms  of  co.     Let  us  then  put 

W  =  W2  (ft)  -  ft)0,  ft)  -  ft)0)  +  w1  (ft)  -  ft)0) 

+  a  function  of  temperature  only, 

w2  and  w1  being  linear  and  homogeneous  in  each  of  their  con- 
stituents.    This    is    the    most    general    quadratic    function    of  &>. 

Now 

*>  =  Kw%  (£  »  -  «0)  +  £wa  (ft)  -  ft)0,  ^)  +  £»£ 

Putting  then  ft)  =  &)0  and  wVw  =  0  we  get  ^w^=0.  Operating 
on  the  last  by  Scr  (  )  where  a  is  any  vector  we  see  that  w1  =  0. 
Thus  putting 

where  cf>  has  the  value  given  by  the  last  equation  and  is  therefore 
self-conjugate  we  get  the  two  following  equations 

w  =  -  S  (&)  -  ft)0)  <j>  (g)  -  ft)0)/2  +  iv0 (47), 

[as  can  be  seen  by  a  comparison  of  the  last  three  equations 
w0  being  the  function  of  the  temperature]  and 

q-'Mq  =  <f>  (&)  -  tw0) (48). 

When  the  natural  shape  of  the  wire  is  straight  these  become 

iv  =  -So)(f)(o/2+w0 (49), 

q~1Mq  =  (f>(D  (50), 

and  when  further  the  wire  is  truly  uniform  <f)  and  w0  are  constant 
along  the  wire. 

32.  Assuming  the  truth  of  these  restrictions  let  us  conceive  a 
rigid  body  moving  about  a  fixed  point  which,  when  placed  in  a 
certain  position  which  we  shall  call  the  normal  position,  has,  if 
then  rotating  with  any  vector  angular  velocity  &),  a  moment 
of  momentum  =  <j>co  where  (f>  has  the  meaning  just  given.  If  the 
rigid  body  be  made  to  take  a  finite  rotation  q()  q'1  and  then 
to  move  with  angular  velocity  qcoq'1  its  moment  of  momentum 
will  be  qfywf1.  Now  let  a  point  move  along  the  wire  with  unit 
velocity  and  let  the  rigid  body  so  move  in  unison  with  it  that 
when  the  moving  point  reaches  the  point  s  the  rigid  body  shall 
have  made  the  rotation  represented  by  q  (§  29  above).  Thus  by 
the  definition  of  co  and  q  its  angular  velocity  at  any  instant 
is  qcoq'1  and  its  moment  of  momentum  therefore  qcfxoq"1  or  M. 

4—2 
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Now  consider  the  equilibrium  of  the  wire  when  no  external 
force  or  couple  acts  except  at  its  ends.  In  this  case  F  is  constant 
throughout  and  it  is  easy  to  see  (what  indeed  is  a  particular  case 
of  equation  (53)  below)  that 

M'+Vp'F  =  0 (51). 

Interpreting  this  equation  for  our  rigid  body  we  get  as  the  law 
which  governs  its  motion 

d  (vect.  mom.  of  mom.)/dt  =  —  Vp'T. 

Thus  the  rigid  body  will  move  as  if  acted  upon  by  a  constant  force 
F  at  the  end  of  the  unit  vector  p  or — since  this  vector  is  fixed  in 
the  body — as  if  acted  upon  by  a  constant  force  acting  through  a 
point  fixed  in  the  body.  From  this  kinetic  analogue  of  Kirchhoff's 
the  mathematical  problem  of  the  shape  of  such  a  wire  as  we 
are  now  considering,  under  the  given  circumstances,  is  shewn  to 
be  identical  with  the  general  problem  of  the  pendulum  of  which 
the  top  is  a  variety. 

33.  We  will  now  give  the  general  equations  for  any  wire 
under  any  external  actions.  The  comparison  of  the  Quaternion 
treatment  of  this  with  the  Cartesian  as  given  in  Thomson  and 
Tait's  Nat.  Phil.  §  614  seems  to  me  to  be  all  in  favour  of  the 
former. 

The  equations  of  equilibrium  of  an  element  ds  are  with  the 
notation  explained  in  §  29  above 

dF  +  Xds  =  0, 

VdpT  +  dm  +  I*ds  =  0, 

or  dividing  by  ds 

F'  +  X  =  0    (52), 

Vp'T  +  M'  +  L  =  0 (53). 

Operating  on  the  last  equation  by  Vp'(  )  noting  that  p'2  =  —  1  and 

putting  Sp'F  =  —  T  we  get 

F  =  pT+  Ty(M'+L) (54), 

whence  by  equations  (52)  and  (53)  respectively 

X.+  d[p'T+Vp'CM.'  +  I*)}/ds  =  0 (55), 

Sp  (M'  +  L)  =  0  (56). 
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Now  by  equation  (48)  above 

q-1Mq  =  <p(a>-a>0) , (48). 

Also  by  the  definition  of  q 

p'  =  qM~l (57), 

where  \  is  some  given  constant  unit  vector.     Finally  as  we  are 
about  to  prove 

co=2Vq-1q     (58). 


• 


It  is  usual  in  the  Cartesian  treatment  to  leave  the  problem  in 
the  form  of  equations  equivalent  to  the  above  (55)  to  (58),  13 
scalar  equations  for  the  13  unknown  scalars  of  p,  T,  M,  q  and  &>. 
We  can  however  as  we  shall  directly  .in  the  Quaternion  treatment 
quite  easily  reduce  the  general  problem  to  one  vector  and  one 
scalar  equation  involving  the  four  unknown  scalars  of  T  and  q 
in  terms  of  which  all  the  other  unknowns  are  explicitly  given. 

To  prove  equation  (58)*  observe  that 

(q  +  dq)  o-(q  +  dq)'1  =  q(<r  +  Vaxrds)  q~l, 

where  er  is  any  vector.  The  truth  of  this  is  seen  by  noticing  that 
(q  +  dq)  (  )  (q  +  dq)  _1  is  the  operator  that  rotates  any  vector  of 
the  element  s  +  ds  from  its  geometrically  normal  position  to  its 
strained  position.  But  we  can  also  get  to  this  final  position 
by  first  in  the  geometrically  normal  wire  making  the  small  strain 
aids  at  the  given  element  and  then  performing  the  strain  of  the 
wire  up  to  the  point  s.  The  first  process  is  represented  on  the 
left  of  the  last  equation  and  the  second  on  the  right.  Thus 
we  get 

q^dqa  +  ad  (g-1) .  q  =  Vcoads, 

or  v  d  (q'1)  =  —  q~ldqq~l, 

— 1     /  — 1     '  T7" 

q   qcr  —  aq   q  =  Vaxr, 
whence  co  =  2  Vq^q. 

Returning  to  equations  (55)  to  (58)  observe  that  equations 
(57)  and  (58)  give  p  and  to  as  explicit  functions  of  q.  Hence 
by  equation  (48) 

M  =  #(2F-<Z-Y-*>0)(f    (59), 

*  This  could  be  deduced  from  Tait's  Quaternions,  §  356,  equation  (2).     His  € 
is  our  quq~x  and  his  dots  our  dashes. 
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which  gives  M   also  explicitly.     Substituting    for    M  and  p    in 
equations  (55)  and  (56)  we  have 

d 
df> 


X  +  ~  iqXq-'T  +  Vq\q'x  (l*  + 


SqXq'1  (L  + 


1 

ds 


qcp  (2  PgY  -  »0)  q'1 

0 (60), 

=  0 (61), 


q(f)(2Vq-1q  -a}0)q 


-i 


which  are  sufficient  equations  to  determine  T  and  q,  whereupon 
M  is  given  by  equation  (59),  w  by  equation  (58)  and  p  and 
therefore  also  p  by  equation  (57). 

Spiral  springs  can  be  treated  very  simply  by  means  of  the 
above  equations,  but  we  have  already  devoted  sufficient  space  to 
this  subject. 


SECTION    IV. 
Electricity  and  Magnetism. 

Electrostatics — General  Problem. 

34.  Merely  observing  that  all  the  theorems  in  integration 
given  in  the  Preliminary  and  ivth  chapters  of  Maxwell's  treatise 
on  Electricity  and  Magnetism,  Part  I.,  are  easy  particular  cases  of 
equations  (8)  and  (9)  §  6  above,  we  pass  on  to  the  one  application 
of  Quaternions  that  we  propose  to  make  in  Electrostatics. 

This  is  to  find  the  most  general  mechanical  results  arising  from 
Maxwell's  theory  of  Electrostatics,  and  to  see  if  they  can  be 
explained  by  stress  in  the  dielectric.  This  problem  as  far  as  I  am 
aware  has  not  hitherto  in  all  its  generality  been  attacked  though 
the  most  important  practical  cases  have  been,  as  we  shall  see, 
considered  by  Maxwell,  Helmholtz,  Korteweg,  Lorberg  and 
Kirchhoff. 

It  is  necessary  first  of  all  to  indicate  as  clearly  as  possible  what 
I  take  to  be  Maxwell's  theory  of  Electricity. 

He  assumes*  all  space  to  be  uniformly  filled  with  a  certain 
substance  called  Electricity.     Whatever  electrical  actions  take  place 

*  Prof.  J.  J.  Thomson  in  his  paper  on  Electrical  Theories,  BA.  Reports,  1885, 
p.  125,  does  not  credit  Maxwell  with  such  a  definite  and  circumscribed  theory  as 
that  described  in  the  text,  and  he  is  thereby  led  to  find  fault  with  Maxwell's  term 
"Displacement"  and  points  out  that  there  is  an  assumption  made  with  reference 
to  the  connection  between  the  true  current  and  this  polarisation  (displacement). 
He  says  moreover,  "  It  is  rather  difficult  to  see  what  is  meant  in  Maxwell's  Theory 
by  the  phrase  '  Quantity  of  Electricity.'  "  None  of  these  remarks  are  called  for  if 
the  view  I  take  of  Maxwell's  theory  be  correct,  and  these  grounds  alone  I  consider 
sufficient  for  taking  that  view.  The  paper  of  Thomson's  here  mentioned  I  shall 
frequently  have  to  refer  to.  [Note  added,  1892.  In  the  text  I  have  given  much  too 
rigid  a  form  to  Maxwell's  theory.  What  I  have  called  his  theory  I  ought  rather 
to  have  called  his  analogy.  Still  I  think  the  present  foot-note  is  in  the  main  just. 
In  my  opinion  it  is  no  more  and  no  less  difficult  to  see  what  is  meant  in  Maxwell's 
Theory  by  "  Quantity  of  Electricity  "  than  by  "  displacement"  since  the  two  are 
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depend  on  the  continued  or  past  motion  of  this  substance  as  an 
incompressible  fluid.  If  electricity  is  brought  from  a  distance  by 
any  means  and  placed  in  a  given  space  there  must  be  a  displace- 
ment of  the  original  electricity  outwards  from  that  space  and  the 
quantity  of  foreign  electricity  is  conveniently  measured  by  the 
surface  integral  of  that  displacement. 

Dielectrics  are  substances  in  which  this  displacement  tends  to 
undo  itself,  so  to  speak,  i.e.  the  original  electricity  tends  to  go  back 
to  its  primitive  position.  In  conductors,  on  the  other  hand,  there 
is  no  property  distinguishing  any  imported  electricity  from  the 
original  electricity. 

The  rate  of  variation  of  displacement,  whether  in  dielectrics  or 
conductors,  of  course  constitutes  an  electric  current  as  it  is  conve- 
niently called. 

We  have  next  to  consider  a  vector  at  each  point  of  space 
called  the  electro-motive  force,  which  depends  in  some  way  at 
present  undefined  on  the  distribution  of  the  displacement  in  the 
dielectrics,  the  distribution  of  currents  whether  in  dielectrics  or 
conductors,  and  on  extra-electrical  or  semi-electrical  action,  e.g. 
chemical  or  mechanical. 

If  at  any  point  the  electro-motive  force  be  multiplied  by  a 
scalar  the  medium  at  the  point  remaining  (except  electrically) 
unchanged,  the  current  in  the  case  of  conductors  and  the  displace- 
ment in  the  case  of  dielectrics  is  altered  in  the  same  ratio.  In 
other  words  the  current  or  the  displacement,  as  the  case  may  be, 
is  a  linear  vector  function  of  the  electro-motive  force,  and  the 
coordinates  of  the  linear  vector  function*  at  any  point  depend 

connected  by  perfectly  definite  equations.  Of  course  it  is  wrong  to  define  "  dis- 
placement" as  "  displacement  of  quantity  of  electricity"  and  then  to  define  "quan- 
tity of  electricity"  in  terms  of  "displacement,"  but  Maxwell  does  not  seem  to  me 
even  tacitly  to  do  this.  Rather  he  says — the  dielectric  is  polarised  ;  this  polarisa- 
tion can  be  represented  by  a  vector  D  ;  electrical  quantity  can  be  expressed  in 
terms  of  D ;  the  mathematical  connections  between  electrical  quantity  and  D 
are  the  same  as  those  between  quantity  of  matter  in  a  space  and  the  displacement 
out  of  that  space  made  by  other  matter  to  make  room  for  the  given  matter ;  we 
will  impress  this  useful  analogy  firmly  on  our  minds  by  calling  D  the  displacement. 
But  I  have  expressed  my  present  views  on  the  meaning  of  Maxwell's  theory  much 
more  fully  in  Phil.  Trans.  1892,  p.  685.] 

*  This  frequently  recurring  cumbrous  mode  of  description  must  be  tolerated 
unless  a  single  word  can  be  invented  for  "a  linear  vector  function  of  a  vector." 
Might  I  suggest  the  term  "  Hamiltonian  ? "  Thus  we  should  say  that  the  dis- 
placement is  a  Hamiltonian  of  the  electro-motive  force,  the  Hamiltonian  at  any 
point  being  a  function  of  the  state  of  the  medium. 
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solely  on  the  state  of  the  medium  (whether  fluid,  solid,  &c.,  or 
again  strained  or  not)  at  that  point. 

To  complete  the  theory  we  have  to  explain  how  the  part  of  the 
electro-motive  force  which  is  a  function  of  the  distribution  of  dis- 
placement and  current  depends  on  this  distribution.  This  expla- 
nation is  obtained  by  making  the  assumption  that  the  electro- 
motive force  bears  to  electricity  defined  as  above  exactly  the  same 
energy  relation  as  ordinary  force  does  to  matter,  i.e. — 

(wk.  dn.  on  electricity  moved  during  any  displacement) 

=  (total  displacement  of  elect.)  x  (resolved  part  of  e.m.f.  in  the 

direction  of  displacement). 

In  the  ivth  part  of  Maxwell's  treatise  he  gives  complete  inves- 
tigations of  the  mechanical  results  flowing  from  this  theory  so  far 
as  it  refers  to  currents,  but  he  has  not  given  the  general  results  in 
the  case  of  Electrostatics.  Nor  has  he  shewn  satisfactorily,  it 
seems  to  me,  that  the  ordinary  laws  of  Electrostatics  flow  from  his 
theory.     It  is  these  investigations  we  now  propose  to  make. 

35.  Our  notation  will  be  as  far  as  possible  the  same  as 
Maxwell's.  Thus  for  the  displacement  at  any  point  we  use  D,  and 
for  the  E.m.f.  E.  From  the  connection  explained  in  last  section 
between  D  and  E  we  have 

D  =  ZE/4tt  (1), 

where  K  at  any  point  is  some  linear  vector  function  depending  on 
the  state  of  the  medium  at  the  point.  If  the  medium  change  in 
any  manner  not  electrical,  e.g.  by  means  of  ordinary  strain  K  will 
in  general  also  suffer  change. 

Let  w  be  the  potential  energy  per  unit  volume  due  to  the 
electrical  configuration.  Thus  if  a  small  increment  BD  be  given 
to  D  at  all  points,  the  increment  JJJ8wd%  in  jjjwd%,  the  potential 
energy  of  the  electrical  configuration  in  any  space,  =  work  done  on 
the  electricity  in  producing  the  change, 

.-.     Jffitod%  =  -  JfJSESDd*, 

by  the  relation  stated  in  §  34  existing  between  E  and  D.  Thus 
limiting  the  space  to  the  element  d$ 

Bw  =  -SESB (2). 
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Now  suppose  D  =  nD'  so  that  by  equation  (1)  E  =  hE'  where  E', 
D'  are  corresponding  E.M.F.  and  displacement  respectively.     Thus 

8w  =  -  nSnSE'B'. 

Integrating  from  n  —  0  to  n=  1,  and  finally  changing  D',  E'  into 
D,  E  we  get 

w  =  -SDE/2    (3). 

From  this  we  get 

Bw  =■  -  SE8T>/2  -  SD8E/2, 

so  that  by  equation  (2)  #ESD  =  $DSE  or  by  equation  (1) 

SEK8E  =  SSEKE. 

Hence  (because  E  and  SE  are  quite  arbitrary)  K  is  self- conjugate 
and  therefore  involves  only  six  instead  of  nine  coordinates*. 

In  electrostatics  the  line  integral  of  E  round  any  closed  curve 
must  be  zero,  for  otherwise  making  a  small  conductor  coincide 
with  the  curve  we  shall  be  able  to  maintain  a  current  by  §  34,  and 
so  (by  the  same  section)  constantly  do  work  on  it  (i.e.  as  a  matter 
of  fact  create  heat)  without  altering  the  statical  configuration. 
Hence  E  must  have  a  potential,  say  v.     Thus 

E  =  -  Vv   (4). 

Since  in  an  electrostatic  field  there  is  no  current  in  a  conductor, 
E  =  0  throughout  any  such  conductor  and  therefore  v  =  const. 

36.  The  charge  in  any  portion  of  space  is  defined  as  the 
amount  of  foreign  electricity  within  that  space.  Thus  the  charge 
in  any  space  is  the  surface  integral  of  the  displacement  outwards. 
Thus  if  there  be  a  charge  on  the  element  rfS  of  a  surface  in  the 
dielectric  this  charge  =  Sd%cJDa  +  #d26D;,  where  a,  b  denote  the 
two  faces  of  the  element  (so  that  dXa  =  —  dXi)  and  in  accordance 
with  |  1  above  dSa  points  away  from  the  region  in  which  the  dis- 
placement is  Da.     Thus  a  being  the  surface  density 

a  =  [SBUdt]a+b (5), 

where  the  notation  []a+b  is  used  for  []a  +  [  ]b.  Similarly  if  there 
be  finite  volume  density  of  foreign  electricity,  i.e.  finite  volume 
density  of  charge  in  any  space,  the  charge  =  —JJSDdX  =  —  j$jSVDd%, 
so  that  if  D  be  the  volume  density 

D  =  -SVD (6). 

*  We  see  from  this  that  H  =  j,yw  or  E=jjv"'«>  according  as  le  is  looked  upon  as 
a  function  of  E  or  D. 
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[The  reason  for  having  +  Sd^aT)a  before  and  -  SdU'D  here  is  that 
in  the  former  case  we  were  considering  a  charge  outside  the  region 
where  Dft  is  considered — between  the  regions  a  and  b  in  fact — 
whereas  in  the  latter  case  we  are  considering;  the  charge  inside  the 
region  where  D  is  considered.  The  same  explanation  applies  to 
the  sign  of  —  $D  TJd%  for  the  surface  of  a  conductor  given  below.] 

In  conductors,  as  we  saw  in  §  34,  the  displacement  has  virtually 
no  meaning  (except  when  it  is  changing  and  so  the  phenomenon 
of  a  current  takes  place)  for  the  foreign  electricity  and  the  original 
electricity  are  not  to  be  distinguished.  Not  so  however  with  the 
surface  of  the  dielectric  in  contact  with  the  conductor.  We  may 
therefore  regard  the  electricity  within  the  body  of  the  conductor 
as  the  original  electricity  so  that  the  charge  is  entirely  at  the 
surface.  Thus  the  surface  density  will  be  —  SDUdX  where  dX 
points  away  from  the  conductor  and  D  is  the  displacement  in  the 
dielectric.  This  may  be  regarded  as  a  particular  case  of  equation 
(5)  D  being  in  accordance  with  what  we  have  just  said  considered 
as  zero  in  the  conductor. 

37.  All  the  volume  integrals  with  which  we  now  have  to  deal 
may  be  considered  either  to  refer  to  the  whole  of  space  or  only  to 
the  dielectrics,  a,s  the  conductors  (except  at  their  surfaces)  in  all 
cases  contribute  nothing.  The  boundary  of  space  will  be  con- 
sidered as  a  surface  at  infinity  and  all  surfaces  where  either  D  or  E 
is  discontinuous. 

Putting  W  =  JJJtudS  we  have  already  found  one  expression  for 
W,  viz. 

2W=-JJJSDEd%. 

We  now  give  another.     By  equation  (4) 

2W  =  fffSl>Vvds 

=  JJvSd^-JJjvSVDds, 
by  equation  (9)  §  (J  above.     Thus  by  equations  (5)  and  (6)  §  36, 

2W  =  JJv<rds+JJJvDd% (7), 

where  ds  is  put,  as  it  frequently  will  be,  for  an  element  of  surface, 
i.e.  Tdl,.  The  value  of  W  which  we  shall  use*  is  obtained  by 
combining  these  two,  viz. 

W  =  fjvads  +  JJJvDdS  +  \  JjJSBEd% (8). 

*  This  is  for  the  general  case  following  the  example  of  Helmholtz  in  the  par- 
ticular case  when  K  reduces  to  a  single  scalar.  See  Wiss.  Abh.  vol.  I.  equation  (2d), 
p.  80.j.     The  method  adopted  in  the  following  investigation  is  also  similar  to  his. 
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So  far  we  have  merely  been  shewing  that  all  the  above  results 
of  Maxwell's  flow  from  what  in  §  34  has  been  described  as  his 
theory.  We  now  proceed  to  the  actual  problem  in  hand  which  is 
proved  from  these  results  however  they  may  be  obtained.  I  may 
remark  that  some  such  investigation  as  the  above  seems  to  me 
necessary  to  make  the  logic  of  Maxwell's  treatise  complete. 

38.  Suppose  now  that  W  is  the  potential  energy  of  some 
dynamical  system  extending  throughout  space.  Let  us  give  to 
every  point  of  space  a  small  displacement  Br]  vanishing  at  infinity 
and  find  the  consequent  increment  BW  in  W.  If  this  can  be  put  in 
the  form 

BW=-SJJSBVl<l>Vld% (9), 

we  shall  have  the  following  expression  for  F  the  force  per  unit 
volume  due  to  the  system 

F  =  0A (10), 

and  the  following  expression  fur  F&.  the  force  per  unit  surface  at 
any  surface  of  discontinuity  in  cf> 

Ts  =  -[cpUdX]a+b (11), 

the  notation  being  the  same  as  in  equation  (5)  §  36. 

Moreover  if  </>  be  self-conjugate  the  forces  both  throughout  the 
volume  and  at  surfaces  of  discontinuity  are  producible  by  the 
stress  <f>  as  can  be  seen  by  §  13  above.  [Compare  all  these  state- 
ments with  §  14  above.] 

For  proof,  we  have  by  equation  (9)  §  6 

B  W  =  -  PM^i*  =  -  /JB&tfcG  +  fffS8v<j>1V1dz, 

where  of  course  the  element  dX  is  taken  twice,  i.e.  once  for  each 
face.     But 

BW  =  —  (work  done  by  the  system  F,  Fg  of  forces) 

=  Jj8BvFsds  +  JJJSBr]Fd%> 

where  the  element  ds  is  taken  only  once.  Equating  the  coefficients 
of  the  arbitrary  vector  Br)  for  each  point  of  space  we  get  the 
required  equations  (10)  and  (11). 

39.  We  must  then  put  B  W  where  W  is  given  by  equation  (8) 
in  the  form  given  in  equation  (9). 

We  must  first  define  B  when  applied  to  a  function  of  the  posi- 
tion of  a  point.     Suppose  by  means  of  the  small  displacement  By 
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any  point  P  moves  to  P'.  Then  Q  being  the  value  at  P,  before 
the  displacement  8t],  of  a  function  of  the  position  of  a  point,  Q  +  8Q 
is  denned  as  the  value  of  the  function  at  P'  after  the  displace- 
ment. Thus  even  in  the  neighbourhood  of  a  surface  of  discon- 
tinuity 8Q  is  a  small  quantity  of  the  same  order  as  Srj. 

Now  the  charge  within  any  space,  that  is  the  quantity  of 
foreign  electricity  within  tha*t  space  will  not  be  altered  by  the 
strain. 

.-.     8(Dd%)=0,      8(<rds)  =  0 (12). 

To  find  SV  we  have 

S(dp+8dp)(V  +  8\7).=SdpV. 

or,  since  8dp  =  —  SdpV  .  8tj, 

SdpBV  =  SdpVfi&vy, 

whence  Sv7=V1^1V  (13). 

The  part  of  8W  depending  on  the  first  two  terms  of  equation 
(8)  is  by  equation  (12) 

fffSvDds  +  JJSvads 
=  _  JffSvSVDdS  +  JJSvSd^jy      [by  equations  (5)  and  (6)  §  36] 
=  JJJSl)V8vd%  [by  equation  (9)  §  6]. 

Noticing  that  8d%  =  -  d%SV8r)  and  that  4ttSD  =  K8E  +  8KE  we 
see  that  the  last  term  in  equation  (8)  contributes 

///SDSEf/s  -$ftfSDHSV$vd%  +  (Stt)-1  jffSEBKEd*. 

Combining  the  last  result  with  the  first  term  of  this  we  get 

fffSD  (VSv  -  8  (Vv))  d$  =  -  tffSD8V.vd% 
=  -fffSDV1S$rltVvd% 

Thus  we  have 

SW=-fffS8Vl{±V1SDE-ESV1l))d%  +  (Stt)"1  JffSESKEdi...(U). 

40.  Now  the  increment  8K  in  K  is  caused  by  two  things 
viz.  the  mere  rotation  of  the  body  and  the  change  of  shape  of  the 
body.     Let  us  call  these  parts  8Kr  and  8KS  respectively. 

First  consider  8Kr.  Suppose  the  rotation  is  e  so  that  any 
vector  which  was  &>  becomes  thereby  co  +  Few.     Thus  the  result 
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of  operating  onw+  Veco  by  K  +  hKr  is  the  same  as  first  operating 
on  a  by  K  and  then  rotating.     In  symbols 

( K  +  BKr)  (ft)  +  Veto)  =  Zo)  +  FeZa>, 

whence  hl\,w  =  Fe/vtw  —  K  (Vew). 

Thus  SESAME  =  SEeZE  -  SEZ  ( FeE ) 

=  2£e7v'EE  =  8tt  SeDE, 

whence  giving  e  its  value  WSrj/2, 

(Stt)"1  SESKrE  =  SVSyVDE/Z. 

Substituting  in  equation  (14) 

BW=  -  JfJS8Vl  {FV1DB/2  -  ES^D]  d%  +  (8tt)-x  JJJSEBK8Ed%, 
or 

STf  =  fcJP&fcFDVjSds  +  (Stt)"1  ftfSEBK&d* (15). 

It  only  remains  to  consider  BKS.  K  is  a  function  of  the  pure 
strain  of  the  medium  and  BK^  is  the  increment  in  K  due  to  the 
increment  in  pure  strain  owing  to  8r).  Calling  this  increment  of 
pure  strain  $yjr  so  that  by  equation  (3)  §  12  above 

Byfra  =  -  8Vl8(oV  J2  -  V  fiaBvjZ, 
we  have 

BK8  =  -8B^a^-^  (l6)> 

by  equation  (7)  §  5  above.     This  gives  by  equation  (5)  §  3 

BK8  =  -8BVli,a^1.^  (17), 

SESh\E  =  -  A'^^^V^A^E (18). 

Now  by  equations  (1)  and  (3)  §  35, 

8ttw  =  -  #E7vE, 

so  that  w  is  a  function  of  the  independent  variables  E,  yjr  (because 
K  is  a  function  of  i/r).     Therefore 

(87r)_15ESZ'aE  =  S&^flwV,  (19). 

This  equation  might  have  been  deduced  at  once  thus 

(H7r)-\SE8KhE  =  S8Tlr^aw$=SB>nl+DwVv 

but  equation  (17)  is  itself  of  importance  so  the  above  proof   is 
preferable. 

Thus  finally  from  equation  (15)  wc  get 

SF  =  ///^&?1(FDVJE/2  +  ^(iwV1)ds (20). 
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We  therefore  have  for  (f>  in  equations  (10)  and  (11)  §  38, 

<j)co  =  -  PTDooB/2  -  ^QWw (21). 

This  is  a  self-conjugate  function  so  that  as  we  saw  in  §  38  it  is 
a  stress  which  serves  to  explain  forces  both  throughout  the  volume 
of  the  dielectric  and  over  any  surfaces  of  discontinuity  in  D  or  E*. 

The  force  in  particular  cases. 

41.  Let  us  first  consider  that  part  —  FDAE/2  of  the  force 
(equations  (10)  §  38  and  (21)  §  40)  which  does  not  depend  on  the 
variation  of  K  with  the  shape  of  the  body. 

Suppose  our  dielectric  is  homogeneous  and  electrically  isotropic 
so  that  K  is  a  simple  constant  scalar.     In  this  case 

4ttD  =  -  KVv  (22), 

by  equations  (1)  and  (4)  §  35.  Therefore  by  equations  (5)  and 
(6)  §  36, 

4>irD  =  KV2v    (23), 

47TO-  =-K[SUdtVv]a+b  (24). 

From  these  we  at  once  get  by  the  theory  of  potential  that 

Kv  =  JffuDds  +  JJuads (25). 

From  this  we  know  by  the  theory  of  potential  that  at  the  surface 
where  the  charge  a  resides  Vv  is  discontinuous  only  with  regard 
to  its  normal  component  and   at  all  other  points   is  continuous. 

Thus 

V  va  =  Vi>6  +  x  UdSa 

and  by  equation  (24)  x  =  ^ira-jK  so  that 

(W^K  ( Vva  -  Vvb)  =  a  Udta , 

whence  (^y'KUdSa  (Vva  -  Vvb)  =  -<r  (26). 

Now  the  force  F  per  unit  volume  is 

-  FDAE/2  =  -  KWvAVv/S-tt 
=  -  KV2vVvJ4nr, 

or  F  =  -DVv    (27), 


* 


As  far  as  I  am  aware  nobody  has  hitherto  attempted  to  find  the  electrical 
forces  much  less  the  stress  except  in  the  case  when  D  is  parallel  to  E  i.e.  the 
dielectric  is  electrically  isotropic  when  unstrained.  The  particular  results  contained 
in  §  45  below  have  been  obtained  by  Korteweg,  Lorberg  and  Kirchhoff  as  is  stated 
in  Prof.  J.  J.  Thomson's  paper  (p.  155)  referred  to  in  §  34. 
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and  the  force  per  unit  surface  Fs  is  by  equation  (11)  §  38, 
i  [VJ)UdfE]a+b  =  (Hirr'KlVVvUd-ZVvl^ 
=  (8TrrKV(VvaUdtaVva  -  VvbUdtaVvb) 
=  (SttPK  V(Vvn  +  Vvb)  UdSa  (Vva  -  Vvb), 

whence  by  equation  (26) 

Fs  =  -i<r[Vv]a+b  (28). 

Thus  we  see  that  Maxwell's  theory  as  given  in  §  34  above 
reduces  to  the  ordinary  theory  when  K  is  a  single  scalar.  In  fact 
two  particles  containing  charges  e  e'  apparently  repel  one  another 
with  a  force  ee  /Kr2  where  r  is  the  distance  between  them,  for  by 
equations  (25),  (27)  and  (28)  the  force  in  any  charged  body  is  that 
due  to  a  field  of  potential  v  given  by 

Kv  =  2ue  (29). 

42.  If  the  medium  when  strained  remain  electrically  isotropic 
^(jK  as  well  as  K  must  be  a  simple  scalar.  Thus  with  Thomson 
and  Tait's  notation  for  strain,  which  makes  the  coordinates  of 
yfr,  e,f,  g,  a/2,  6/2,  c/2  we  have 

^K  =  dK/de  =  dK/df=  dKjdg 
dKjda  =  dK/db  =  dK/dc  =  0. 

Therefore  K  is  &  function  of  e  +f+g  only,  i.e.  of  the  density  (m) 
of  the  medium.     Thus  because 

de  +  df+  dg  =  —  dm/m  =  —  d  log  m 

we  get  ^(jK  =  —  dK/d  log  m  =  —  k  suppose.     Hence 

^qw  =  (SttY'SE  (dK/d  log  m)  E  =  A:E2/87r. 

Thus  the  force  —  ^qwA  [equations  (21)  §  40  and  (10)  §  38]  result- 
ing from  the  change  of  K  with  pure  strain  is  in  the  case  we  are 
now  considering 

-&VE8/87r (30), 

and  is*  therefore,  according  as  k  is  positive  or  negative,  in  the 
direction  of  or  that  opposite  to  that  of  the  most  rapid  increase  of 
the  square  of  the  electromotive  force.  Thus  even  in  the  case  of  a 
fluid  dielectric  which  has  no  internal  charge  but  which  forms  part 
of  a  non-uniform  field  of  (electromotive)  force  the  surfaces  of  equal 
pressure  and  therefore  the  free  surface  will  if  originally  plane  no 
longer  remain  so. 

*  This  is  the  same  result  as  Helmhnltz's  on  the  same  assumption  Wiss.  Abh. 
i.  p.  798. 
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Nature  of  the  Stress. 

43.  We  have  seen  that  the  stress  which  serves  to  explain  the 
electrostatic  forces  is  that  given  by  equation  (21)  §  40,  viz. 

<f)a)  =  -  VDcoE/2  -  ^qwco (21). 

Let  us  first  consider  the  part  —  FDwE/2  which  does  not 
depend  on  the  variation  of  K.  Putting  co  first  =  Z7D  and  then 
=  UE  we  get 

<f>TTD  =  TBTE .  Z7E/2, 

<f>UE  =  Tl>TE.Ul>/2. 

Therefore  putting  <u  first  =  any  multiple  of  Z7D  +  UE  and  then 
=  any  multiple  of  t/"D  —  UE  we  get 

(jxo  =  TBTEco/2, 

<f>a  =  -  T-DTEco/2. 

Lastly,  since 

-  VB(oE  =  (oSBE  -  BScoE  -  E&uD, 

we  see  that  if  we  put  w  =  any  multiple  of  FDE 

4>(jo  =  w8DEj2  =  -  cow. 

Thus  we  see  that  the  stress  now  considered  is  a  tension  aloncj 
one  of  the  bisectors  of  D  and  E  (the  bisector  of  the  positive  direc- 
tions or  the  negative  directions  of  both)  =  TDTE/2,  an  equal 
pressure  along  the  other  bisector  and  a  pressure  =  w  perpendicular 
to  both  these  directions.  When  D  is  parallel  to  E  this  at  once 
reduces  to  Maxwell's  case,  viz.  a  tension  in  the  direction  of  E  and 
a  pressure  in  all  directions  at  right  angles  each  =  w. 

44.  We  have  now  to  consider  the  other  part  of  the  stress,  viz. 

<f)CO  =  —  ^Qivw, 
or  <f>G)  =  ^(oSEK^/Stt     (31). 

If  we  assume  that  K  is   a  function  of   the   density  (m)  of  the 
medium  only  we  shall  have 

dKJde  =  dK/df=  dK/dg  =  -  dK/d  log  m  =  -  k, 
say,  and 

dK/da  =  dK/db  =  dK/dc  =  0, 

as  in  §  42.     Here  however  k  is  not  in  general  a  mere  scalar  but  a 
self- conjugate  linear  vector  function.     We  have  then  in  this  case 

(/>&)  =  -  (oSEkE/S7r, 
M.  5 
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which  is  a  hydrostatic  pressure  or  an  equal  tension  in  all  direc- 
tions according  as  $E&E  is  positive  or  negative.  In  this  case  the 
36  coordinates  of  ^(J^w  .  Kx  reduce  to  the  6  of  k  for  each  point  of 
space. 

A  more  general  assumption  is  that  8KS  (§  40)  depends  only  on 
the  elongations  in  the  directions  of  the  principal  axes  of  K. 
Taking  i,j,  k  as  unit  vectors  in  these  directions  we  again  have 

dK/da  =  dK/db  =  dK/dc  =  0, 

and  thus  <f>i  =  iSE  (dK/de)  E/8tt, 

and  similarly  for  j  and  k,  so  that  the  principal  axes  of  the  stress 
now  considered  are  the  principal  axes  of  K. 

45.*  The  most  natural  simple  assumption  for  solid  dielectrics 
seems  to  me  to  be  that  the  medium  is  electrically  isotropic  before 
strain,  and  also  isotropic  with  regard  to  the  strain  in  the  sense 
that  if  the  strain  be,  so  to  speak,  merely  rotated,  8KS  will  suffer 
exactly  the  same  rotation.  We  may  treat  this  problem  exactly  as 
we  did  (§  20)  that  of  stress  in  terms  of  strain  for  an  isotropic  solid. 
Thus  splitting  up  Si|r  into  its  principal  elongations,  i.e.  putting 

8-v^w  =  —  'SSeiSia), 
we  shall  get,  as  in  §  20, 

8Ks(o  =  —  (a  —  ft)  'ZSeiSico  +  ftcoXSe 
=  (a-ft)Byjr(o-ftcoS^^. 

But  SKs(o  —  —  SSyfr^(j^ .  Kja  by  equation  (16)  §  40,  so  that  from 
equation  (31) 

-  87nS&KK  =  («-£)  SES^E  -  /3E2£^£ 

whence  we  see  by  §  4  above  that 

<fxo  =  {(a  -  ft)  ESEco  +  /3E2<»}/87r (32), 

which  consists  of  a  pressure  in  the  directions  of  the  lines  of  force 
=  —  aE2/87r  and  another  pressure  in  all  directions  at  right  angles 

=  -  /3E2/87T. 

*  For  references  to  former  proofs  of  this  see  foot-note  to  §  40  above. 
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Magnetism. 

Magnetic  potential,  force,  induction. 

46.  We  now  go  on  to  the  ordinary  theory  of  magnetism  ;  and 
here  we  shall  merely  follow  Maxwell  in  his  General  Theory,  so  as 
to  give  an  opportunity  of  comparing  Quaternion  proofs  with  Car- 
tesian, as  we  have  already  done  in  Elasticity. 

We  shall  not  consider  in  detail  the  effect  of  one  small  magnet 
upon  another,  as  this  has  already  been  done  by  Tait.  In  connec- 
tion with  this  I  am  content  to  remark  that  I  think  the  treatment 
of  this  problem  can  be  made  somewhat  simpler  than  Tait's  by 
means  of  potential. 

Suppose  we  have  a  pole  -mat  0  and  a  pole  +  m  at  0'  where 
00'  is  small.  Calling  the  vector  from  0  to  0'  00',  let  us  call  the 
vector  mOO'  fju,  so  that  fi  is  the  vector  magnetic  moment  of  the 
magnet.  The  potential  of  —  m  at  any  point  P  is  —  mu,  where  u  as 
usual  =  PO'1.  Similarly  the  potential  of  +  m  is  mu',  where 
u'  =  P0'~l.     Therefore  the  potential  of  the  magnet 

=  m  (V  —  u)  =  mSO  O'Vu, 

where  of  course  P  is  the  variable  point  implied  by  V.  Thus  the 
potential  of  a  small  magnet  /jl  at  any  j)oint  =  S/mVu. 

Hence  the  potential  of  any  magnet  whose  magnetic  moment 
per  unit  volume  at  an}7  point  is  I  is 

O  =  SVfffuIdZ  =  -  JffSlVudz (33), 

according  to  the  convention  of  §  9  above.  By  equation  (9)  §  6  this 
may  be  put 

n  =  -  JfvSldX  +  fffuSVId*  (34), 

which  shews  that  we  may  consider  it  due  to  a  volume  density 
$VI  and  a  surface  density  —  SlUdX*  of  magnetic  matter,  the  sur- 
face density  occurring  wherever  there  is  discontinuity  in  I. 

By  again  considering  the  poles  m  and  —  m  of  the  small  magnet 
IX  we  see  that  its  potential  energy  when  placed  in  a  field  of  mag- 
netic potential  12  is  —  S/xVD,,  whence  just  as  we  obtained  equation 

*  [Note  added,  1892.     More  generally  and  better  -  [SlUd1,]a+h.] 

5—2 
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(33)  we  now  see  that  the  potential  energy  (W)  of  any  magnet  in 
such  a  field  is 

w  =  -jfjsivnd%  (35), 

or  W=-  fjClSldZ  +  fjfnSVlds (36), 

by  equation  (9)  §  6  above,  so  that  the  potential  energy  is  just  the 
same  as  it  would  be  for  the  imaginary  distribution  of  magnetic 
matter. 

47.  The  force  (H)  on  a  unit  magnetic  pole  at  any  point 
external  to  the  magnet  is  given  by 

H  =  -  VH  =  -  VSVfffuIdS  =  -  V*fffuld%  +  V  VVJJJuIds 

=  -  4ttI  +  VjjJVV7ud%  =  -  4ttI  +  VA, 
where 

A=fffVlVud% (37). 

Thus  we  see  that  for  all  external  points  H  =  VA,  so  that  A  is 
called  the  vector  magnetic  potential.  [It  is  to  be  observed  that 
since  VA  =  H  +  47rl  =  a  vector,  #VA  =  0.]  VA  is  called  the  mag- 
netic induction  and  for  it  we  use  the  single  symbol  B  so  that 

B  =  VA   (38). 

Thus  £VB  =  SV2A  =  0    (39), 

and  also  by  the  equation  for  H  just  given 

B  =  H  +  4ttI    (40). 

This  is  not  the  way  in  which  Maxwell  defines  the  magnetic 
force  and  induction,  but  lie  shews  quite  simply  (Elect,  and  Mag. 
§§  398 — 9)  that  his  definition  and  the  present  one  are  identical. 
This  can  be  shewn  as  easily  without  analysis  at  all. 

48.  Where  I  is  discontinuous  both  H  and  B  are  also  discon- 
tinuous. From  the  surface  density  view  we  gave  in  equation  (34) 
we  see  that,  just  as  we  have  the  expression,  given  in  §  41  for 
Vva—Vvb,  so  now 

Ub-Ha  =  -4>7rUdla[SlUd^}a+b (41), 

so  that  the  discontinuity  in  H  is  entirely  normal  to  the  surface  of 
discontinuity.     Further  from  this  equation  we  have 

S(Ht  -  H„)  dZa  =  4tt  [SldX]a+h, 

i.e.  [£(H  +  47rI)eZ2]a+&  =  0, 

or  [SBdZ\a+b  -  0 (42), 

so  that  the  discontinuity  in  B  is  entirely  tangential. 
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From  this  equation  we  see  that  for  any  closed  surface  whatever 
whether  it  include  surfaces  of  discontinuity  in  I  or  not 

JfSdSB  =  0. 

For  adding  these  surfaces  to  the  boundary  of  the  inclosed  space, 
in  accordance  with  §  7  above,  we  see  by  equation  (42)  that  they 
contribute  zero  to  the  surface  integral ;  but  the  total  surface  inte- 
gral is  by  equation  (9)  §  6  fjfSVBd%  =  0  by  equation  (39). 


Magnetic  Solenoids  and  Shells. 

49.  A  magnet  is  said  to  be  solenoidal  if  the  imaginary  mag- 
netic matter  of  equation  (34)  is  entirely  on  the  surface.  Thus 
for  a  solenoidal  distribution 

£VI  =  0   (43). 

In  this  case  the  potential  is  by  equation  (34) 

a  =  -jjuSIdt (44). 

50.  A  simple  magnetic  shell  is  defined  as  a  sheet  magnetised 
everywhere  normally  to  itself  and  such  that,  at  any  point,  the 
magnetic  moment  per  unit  surface  is  a  constant  called  the  strength 
of  the  shell. 

Calling  the  strength  cf>  we  have  for  the  potential  energy  at  any 
point  by  equation  (33) 

Q,  =  -<f>fJSd$Vu (45). 

Now  —  SdXVu  is  the  solid  angle  subtended  by  the  element  dS  at 
the  point  considered,  so  that 

12  =  </>  x  (solid  angle  subtended  by  shell  at  point) (46). 

Thus  if  P  be  a  point  on  the  positive  side  of  the  shell  and  P  a 
point  infinitely  near  P  but  on  the  negative  side 

Potential  at  P  —  Potential  at  P'  =  4nr(p, 

or  what  comes  to  the  same  thing 

-  )     SUdp  =  4tt(£. 
J  p 

This  integral  may  be  taken  along  any  path,  e.g.  along  a  path 
which  nowhere  cuts  the  shell.  The  same  integral  is  true  if  H  be 
the  magnetic  force  due  to  a  whole  field  of  which  the  shell  is  only 
one  of  several  causes,  for  the  part  contributed  by  the  rest  of  the 
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field  is  zero  on  account  of  the  infinite  proximity  of  P  and  P'. 
For  future  use  in  electro-magnetism  observe  that  this  statement 
cannot  be  made  if  for  H  in  the  integral  be  substituted  B. 

51.  The  condition  that  any  magnet  can  be  divided  up  into 
such  shells  is  at  once  seen  to  be  that  I  can  be  put  in  the  form 

I  =  V<£  (47), 

where  <p  is  some  scalar. 

In  this  case  the  potential  is  by  equation  (33) 

O  =  -JffSlVud%  =  -JffSVcpVudz, 

or  by  equation  (9)  §  6 

a  =  -  JftSdZVu  +  JJf4>V2ud%, 

i.e.  Q,  =  -  ff<f>SdZVu  +  4ir<l> (48). 

Remarking  that  the  solid  angle  again  occurs  here  it  is  needless  to 
interpret  the  equation  further.  By  equation  (37)  we  have  for  the 
vector  potential 

A  =JJJVlVud%  =  JJJVV(pVud%, 

or  by  equation  (9)  §  6 

A=ff<f>Vd%Vu (49). 

52.  The  potential  energy  of  a  magnetic  shell  of  strength  <p 
placed  in  a  field  of  potential  O  is  of  importance.     We  see  by 

■  equation  (35)  that  it  is 

W=-<f>fJSd'2Va 

If  then  the  magnets  which  cause  12  do  not  cut  the  shell  anywhere, 
so  that  —  V12  =  VA,  we  shall  have 

W  =  cf>fjSd$VA  =  <f>jjSd%B (50), 

or  W  =  <f>JSdPA (51), 

by  equation  (8)  §  6. 

Suppose  now  that  A  is  caused  by  another  shell  of  strength  </>'. 
Then  by  equation  (49) 

A  =  f  JJVdZVu  =  4>'judp, 

by  equation  (8)  §  6.  Thus  finally  the  potential  energy  M  of  these 
two  shells  is  given  by 

M  =  <pcp'  JJuSdpdp' (52). 

53.  The  general  theory  of  induced  magnetism  when  once  the 
proposition  (given  in  equation  (42)  §  48)  that  [$c2XB]a+6  is  zero  is 
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established,  is  much  the  same  whether  treated  by  Quaternion  or 
Cartesian  notation.  We  shall  therefore  not  enter  into  this  part  of 
the  subject. 


Electro-magnetism. 

General  Theory. 

54.  We  now  propose  to  prove  the  geometrical  theorems  con- 
nected with  Maxwell's  general  theory  of  Electro-magnetism  by 
means  of  Quaternions. 

We  assume  the  dynamical  results  of  Chaps.  V.,  VI.  and  VII., 
and  the  first  six  paragraphs  of  Chap.  VIII.  of  the  fourth  part  of 
his  treatise. 

These  assumptions  amount  to  the  following.  Connected  with 
any  closed  curve  in  an  electro-magnetic  field  there  is  a  function 

p  =  -JSAdP  (53), 

where  A  is  some  vector  function  at  every  point  of  the  field.  The 
function  p  has  the  following  properties.  If  any  circuit  be  made 
to  coincide  with  the  curve  the  generalised  force  acting  upon  the 
electricity  in  the  circuit  is 

E=-p (54). 

Again,  if  there  be  a  current  of  electricity  7  flowing  round  this 
circuit,  the  generalised  force  X,  corresponding  to  any  coordinate  x 
of  the  position  of  the  circuit  due  to  the  field  acting  upon  the  con- 
ductor, is 

X  =  ydp/dx    (55). 

bb.  The  first  thing  to  be  noticed  is  that  p  can  be  transformed 
into  a  surface  integral  by  equation  (8)  §  G  above. 

Thus  p  =  -ffSBd2 (56), 

where  B  =  FVA    (57), 

so  that  *S'VB  =  0. 

Next  we  see  by  the  fundamental  connection  (§  34  above) 
between  the  E.M. f.  E  and  electricity,  that  E  must  equal  the  line 
integral  of  E  round  the  circuit,  or 

E=-jSEdp (58). 
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Wc  arc  now  in  a  position  to  find  El  in  terms  of  A  and  B,  i.e. 
of  A. 

56.  The  rate  of  variation  of  p  is  due  to  two  causes,  viz.  the 
variation  of  the  field  (A)  and  the  motion  of  the  circuit  (p).  In 
the  time  St  then  there  will  be  an  increment  8A  in  A  and  an  incre- 
ment SdX  in  d'Z  to  be  considered.     Thus 

pSt  =  -  JSSAdp  -  Jj'SBSdX. 

[This  amounts  to  assuming  that  JSSAdp  =  JfSSlBd^,  which  of 
course  is  true  by  equation  (8)  §  6.]  Now  when  the  circuit  changes 
slightly  we  may  suppose  the  surface  over  which  the  new  integral 
extends  to  coincide  with  the  original  surface  and  a  small  strip  at 
the  boundary  traced  out  by  the  motion  (pSt)  of  the  boundary. 
Thus  SdX  is  zero  everywhere  except  at  the  boundary  and  there  it 

=  V(pSt)  dp, 

so  that  pSt  =  —  JSSAdp  +  St  JSp"Bdp, 

whence  dividing  by  St 

p=JS(-A+  VpB)dp (59), 

but  by  equations  (54)  and  (58)  p  =JSEdp.     Thus 

E  =  -A+F>B~V^ (60), 

where  yjr  is  a  scalar  and  —  V-v/r  is  put  as  the  most  general  vector 
whose  line  integral  round  any  closed  curve  is  zero. 

57.  We  now  come  to  the  mechanical  forces  exerted  on  an 
element  through  which  a  current  C  per  unit  volume  flows. 

We  see  by  equation  (55)  that  the  work  done  by  the  mechanical 
forces  on  any  circuit  through  which  a  current  of  magnitude  7  flows 
in  an}7  small  displacement  of  the  circuit  equals  7  x  the  increment  in 
p  caused  by  the  displacement.  Give  then  to  each  element  dp  of 
the  circuit  an  arbitrary  small  displacement  8p  and  let  F'  be  the 
mechanical  force  exerted  by  the  field  upon  the  element.  Thus  as 
in  last  section 

-JSF'Sp  =  ySp  =  -  yff&BSdZ  =  -  yfSSpdpB. 

Thus  the  force  F'  on  the  element  dp  is  yVdpB.  But  we  may 
suppose  this  element  to  be  an  element  d$  of  volume  through  which 
the  current  C  flows.  Thus  if  for  ydp  we  write  Gd%,  and  for  F', 
Frfs,  where  F  is  the  force  per  unit  volume  exerted  by  the  field,  we 


get 


F=  VCB    (61). 
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58.  So  far  we  have  been  able  to  go  by  considering  the  electric 
field  as  a  mechanical  system,  but  to  go  further  (as  Maxwell  points 
out)  and  find  how  B  or  A  depends  on  the  distribution  of  current 
and  displacement  in  the  field  we  must  appeal  to  experiment.  It 
has  been  shewn  by  experiment  that  a  small  circuit  produces 
exactly  the  same  mechanical  effects  on  magnets  as  would  a  small 
magnet,  at  the  same  point  as  the  circuit,  placed  with  its  positive 
pole  pointing  in  the  direction  of  the  positive  normal  to  the  plane 
of  the  circuit  when  the  positive  direction  round  the  circuit  is 
taken  as  that  of  the  current*.  Moreover  the  magnetic  moment 
of  the  magnet  which  must  be  placed  there  is  proportional  to  the 
strength  of  the  current  x  the  area  of  the  circuit.  Further,  the 
effect  of  this  circuit  upon  other  such  small  circuits  is  the  same  as 
the  mutual  effects  of  corresponding  magnets.  We  have  now  only 
to  consider  a  finite  circuit  as  split  up  in  the  usual  way  into  a 
number  of  elementary  circuits  to  see  that  a  finite  circuit  will  act 
upon  magnets  or  upon  other  circuits  exactly  like  a  magnetic  shell 
of  strength  proportional  to  the  strength  of  the  current  and 
boundary  coinciding  with  that  of  the  circuit.  The  unit  current 
in  the  electro-magnetic  system  is  so  taken  as  to  make  this  propor- 
tionality an  equality. 

The  one  difference  between  the  circuit  and  the  magnetic  shell 
is  that  there  is  no  discontinuity  in  the  magnetic  potential  in  going 
round  the  circuit,  so  that  by  §  50  above  the  line  integral  of  H 
round  the  circuit  will  be  4vr  x  the  strength  of  the  current.  In 
symbols 

fSRdp  =  4nrffSCd2 

for  any  curve,  so  that  by  equation  (8)  §  6 

4ttC  =  7VH   (02), 

whence  SVC  =  0 (63), 

which  of  course  is  a  direct  result  of  our  original  assumption  that 
electricity  moves  like  an  incompressible  fluid.  Maxwell  tacitly 
assumes  this  by  making  the  assumption  that  only  one  coordinate 
is  required  to  express  the  motion  of  electricity  in  a  circuit. 

*  See  §  1  above  for  the  convention  with  respect  to  the  relation  between  the 
positive  side  of  a  surface  and  the  positive  direction  round  its  boundary.  Hitherto 
we  have  had  no  reason  for  choosing  either  the  right-handed  or  the  left  handed 
screw  as  the  type  of  positive  and  negative  rotation.  But  to  make  the  statement  in 
the  text  correct  we  must  take  the  former. 
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59.  We  are  now  in  a  position  to  identify  the  B  we  are  now 
using  with  the  magnetic  induction  for  which  we  have  already  used 
the  same  symbol. 

We  see  by  equation  (50)  §  52  that  the  mechanical  force  on  the 
shell  corresponding  to  any  coordinate  x  is 

-QdffSdZB'/dx, 

where  B'  is  the  magnetic  induction ;  and  by  equations  (53)  and 
(55)  that  the  force  on  the  corresponding  electric  circuit  is 

-  <f>dfJSdtB/dx, 

therefore  B  =  B'  wherever  there  is  no  magnetism.  And  where 
there  is  magnetism  B  is  not  =  H  for  #VB  =  0,  as  we  have  seen. 
Thus  B  =  B'  at  all  points.  In  other  words  the  two  vectors  are 
identical  and  we  are  justified  in  using  the  same  symbol  for  the 
two. 

This  practically  ends  the  general  theory  of  electro-magnetism. 
We  content  ourselves  with  one  more  application  of  Quaternions 
in  this  subject.  We  give  it  because  it  exhibits  in  a  striking 
manner  the  advantages  of  Quaternion  methods. 


Electro-magnetic  phenomena  explained  by  Stress. 

60.  In  §  46  equation  (35)  we  have  seen  that  the  potential 
energy  of  a  magnetic  element  =  Ids  in  a  magnetic  field  is  SlRdS 
when  H  has  a  potential.  Maxwell  assumes  that  the  same  ex- 
pression is  true  whether  H  have  a  potential  or  not.  Assuming 
this  point*  with  him  we  can  find  the  force  and  couple  acting  on 
the  medium  and  a  stress  which  will  produce  that  force  and 
couple.    The  force  and  couple  due  to  the  magnetism  of  an  element 

*  I  do  not  defend  the  legitimacy  of  this  assumption.  It  seems  to  me  bold  to 
assume  that  a  magnet  possesses  any  such  thing  as  potential  energy  in  a  field  which 
has  no  potential.  If  we  assume  H  and  its  derivatives  to  be  continuous  throughout 
our  typical  element  cZg  of  volume  containing  a  great  number  of  molecules  (both 
material  and  magnetic)  the  force  on  a  magnetic  molecule  /j.  consisting  of  two  poles 
is  -  Syiiy-  H  and  the  force  per  unit  volume  -  Sly.  H,  which  is  only  identical  with  the 
expression  -  Vi^111!  obtained  below  when  H  has  a  potential.  With  this  expression 
a  stress  cannot  be  found  that  produces  the  force.  If,  however,  H  and  its  derivatives 
be  not  assumed  continuous  in  this  manner  the  force  on  the  magnet  /x  is  quite  inde- 
terminate whether  the  magnetic  pole  or  the  molecular  current  view  of  magnetism 
be  taken,  unless  it  be  specified  in  what  way  the  poles  and  currents  are  distributed  in 
the  element  of  volume. 
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is  obtained  by  giving  the  element  an  arbitrary  translation  and 
rotation  and  assuming  that  the  work  done  by  this  force  and 
couple  =  the  decrement  in  the  potential  energy  of  the  element. 
Thus  the  force  per  unit  volume  is  -  V^IHj  for  the  decrement  in 
the  potential  energy  due  to  a  small  translation  Sp  is  SSpV^IH^ 
Similarly  the  couple  M  is  given  by 

M=  FIH (64) 

for  the  decrement  —  $MS&>  in  the  potential  energy  due  to  a  small 
rotation  So>  is  —  SVSwI  .  H  =  —  SScoVlH.  The  total  force  F  per 
unit  volume  is  the  sum  of  that  just  given  and  that  given  by 
equation  (61),  so  that 

F  =  FCB-V^IH,. 

Therefore  by  equation  (62) 

4ttF  =  VVVH. .  B  -  4ttV1^IH1 

=  -  H^B  +  V^H,  -  4ttV1^IH1 
=  -H1^V1B  +  V1>SHH1. 

Now  £VB  =  0  so  that  H^B  =  H#AB,  and  again 

V15fHH1  =  V(H2)/2, 

so  that  F  =  <£A  (65), 

where  Sircfxo  =  -  2H&bB  +  coU2    (66). 

From  this  we  get 

87r7#>£=  2FBH  =  SttFIH, 
so  that  M=  F£</>£ (67). 

From  these  two  results  (equations  (65)  and  (67))  we  see  by 
§13  above  that  the  stress  (p>  will  produce  all  the  mechanical  effects 
of  the  field. 

This  stress,  as  can  be  seen  by  giving  w  the  required  values  in 
equation  (66),  is  one  of  pressure  —  H2/8tt  in  all  directions  at  right 
angles  to  B  and  of  tension  —  $H  (2B  —  H)/87r  in  the  direction  of 
H.  When  there  is  no  magnetism  H  =  B  so  that  this  pressure  and 
tension  become  equal  and  their  directions  at  right  angles  to  and 
along  B  respectively.     In  fact  we  then  have 

8tt<£g>  =  -  BwB     (68). 


SECTION  V. 
Hydrodynamics. 

61.  In  the  applications  I  am  about  to  make  in  this  I  have 
practically  nothing  new  to  shew  except  the  utility  of  Quaternion 
methods  in  the  general  theory  of  Hydrodynamics  in  all  its  parts. 

I  therefore  take  a  treatise  (Greenhill's  article  in  the  Encyc. 
Brit,  on  this  subject)  and  work  out  the  general  theory  on  the  lines 
of  the  treatise.  This  is  more  necessary  than  at  first  sight  it  would 
seem,  for  I  believe  mathematicians  who  have  studied  Quaternions 
are  under  the  impression  that  the  method  does  not  lend  itself 
conveniently  to  the  establishment  and  treatment  of  such  equations 
as  the  Lagrangian  and  those  of  Cauchy's  integrals.  With  our 
meaning  of  V,  however,  the  Quaternion  treatment  of  these  equa- 
tions is  as  much  simpler  than  the  Cartesian  as  in  the  case  of  the 
Eulerian  equations. 

Down  to  equation  (13)  below  the  subject  has  already  been 
handled  by  Prof.  Hicks*  in  his  Quaternion  treatment  of  Strains 
and  Fluid  motion  (Quart.  Journ.  Math.  xiv.  [1877]  p.  271).  I  do 
not  hesitate  to  go  over  the  ground  again  as  my  methods  are 
different  from  his. 


Notation. 

62.  For  the  vector  velocity  at  any  point  we  shall  use  a,  for 
the  density  D,  for  the  force  per  unit  mass  P,  for  its  potential,  when 
it  has  one,  v,  and  for  the  pressure  p.  For  time-flux  which  follows 
a  particle  we  shall  use  d/dt  or  the  Newtonian  dot,  and  for  that 
which  refers  to  a  fixed  point  of  space  d/dtf. 

*  [Note  added,  1892.     Prof.  Tait's  name  ought  to  be  added  to  Prof.  Hicks's.] 
+  [Note  added,  1892.     I  am  aware  that  this  is  contrary  to  the  usual  English 
custom,  but  that  custom — of  interchanging  the  meanings  of  d/dt  and  bjdt  as  given 
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Thus 

d/dt  =  d/dt-Sa-V (1). 

This  equation  is  given  on  p.  446  of  Greenhill's  article  already 
mentioned.  In  future  we  shall  refer  to  this  article  simply  as 
"  Greenhill's  article." 

Eider's  Equations. 

63.  To  find  the  equation  of  continuity,  with  Greenhill,  we 
merely  express  symbolically  that  the  rate  of  increase  of  the  mass  cf 
the  fluid  in  any  space  equals  the  rate  at  which  it  is  flowing 
through  the  boundary.     Thus  M  being  the  mass  in  any  space, 

dM/dt=ffDScrd2 (2). 

This  is  Greenhill's  equation  (1)  p.  445.  By  it  and  equation  (9) 
§  6  above  we  have 

dM/dt=fffSV(D<r)d%, 

whence  reducing  the  volume  to  the  element  d%, 

dD/dt  =  SV(Da) (3). 

This  is  Greenhill's  equation  (2).  Thus  by  equation  (1)  of  last 
section, 

D/D  =  SVa-    (4). 

64.  To  obtain  Euler's  equations  of  motion  we  express  that 
the  vector  sum  of  the  impressed  forces  for  any  volume  equals  the 
vector  sum  of  the  bodily  forces  plus  the  vector  sum  of  the  pressures 
on  the  surface.     Thus 

fffDad%=fffDFdS-ffpdZ 
=  fff(DF-Vp)d%, 

by  equation  (9)  §  6  above.  Applying  this  to  the  element  d%  and 
dividing  by  Dd%  we  have 

&  =  T-VpjD  (5), 

or  by  equation  (1)  §  62 

da/dt-SaV.<T  =  F-Vp/D  (6). 

This  is  Greenhill's  equations  (4)  (5)  (6). 

in  the  text — seems  to  me  out  of  harmony  with  the  meaning  attached  to  r  in  other 
hranches  of  Mathematics.  At  any  rate  I  have  respectable  fellow-sinners,  e.g. 
Kirchhoff  in  his  Mechanik,  zweite  Vorlesung,  et  seq.] 
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If  F(p,  t)  =  0  (F  a  scalar)  be  the  equation  of  a  surface  always 
containing  the  same  particles  F=0,  or  by  equation  (1) 

dF/dt-SaVF=0 (7). 

This  is  Greenhill's  equation  (7). 

65.  Let  us  now  put 

fdp/D  =  P (8). 

This  of  course  assumes  that  D  is  a  function  of  p  only,  which  is  not 
always  the  case,  for  instance  in  a  gas  where  diffusion  of  heat  is 
taking  place.  If  F  have  a  potential  v,  F  =  —  Vv.  Thus  equation 
(5)  becomes 

&  =  -V(v+P) (9), 

and  equation  (6) 

d*/dt-S<rV.v  =  -V(v  +  P)    (10). 

Now 

SaV1 .  at  =  VaVV1a1  +  V^o^  =  VaVVa  +  V  (<r2)/2. 
Thus  equation  (10)  becomes 

d<r/dt  +  2Ve*  +  VR  =  0 (11), 

where  the  scalar  R  is  put  for  P  +  v  —  o-a/2  and  the  vector  e  for 
VVa/2.     This  is  Greenhill's  equations  (8)  (9)  (10). 

If  e  =  0,  i.e.  Wa  =  0,  we  have  <r  =  V<£,  whence  our  equation 
becomes 

V  (dcf>/dt  +  R)  =  0, 
so  that  d(f>jdt  +  R  =  H    (12), 

where  if  is  a  function  of  t  only.  In  next  section  we  shall  obtain 
in  the  case  of  an  infinite  fluid  a  generalisation  of  this  which  I 
believe  has  not  hitherto  been  obtained.  Here  we  have  made  the 
assumption  that  if  Wa  =  0  at  one  epoch  it  will  be  so  always. 
This  we  shall  prove  later. 

66.  Greenhill  next  considers  the  case  of  steady  motion.     In 
this  case  dajdt  =  0,  so  that  equation  (11)  becomes 

2Vea  +  VR  =  0  (13), 

and  therefore  the  surface  R  =  const,  contains  both  vortices  and 
stream-lines  and  the  relation  dR/dn  =  2qco  sin  6  given  on  p.  446  of 
Greenhill's  article  is  the  natural  interpretation  of  our  equation. 

So  far  we  have  been  going  over  much  the    same  ground   as 
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Hicks,  but  now  we  enter  upon  applications  of  Quaternions  that  I 
think  have  not  been  made  before. 

G7.  Greenhill  next  considers  rotating  axes  and  finds  the  form 
of  the  equations  of  motion  when  referred  to  these.  Let  a  be  the 
velocity  referred  to  them  ;  so  that  if  the  axes  have  at  any  time 
made  the  rotation  qOq'1  the  real  velocity  will  be  qo-q'1.  Thus*, 
as  always  with  rotating  axes,  if  a.  be  any  vector  function  of  a 
particle  the  rate  of  increase  of  a  in  space  is 

a.  +  Vcoa, 

where  w  is  the  angular  velocity  referred  to  the  same  system  of 
rotating  axes.  Thus  we  see  that  the  acceleration  of  a  particle 
=  &  +  Vwa.  The  velocity  a  —  p  +  Va>p ',  where  p  is  the  vector 
coordinate  referred  to  our  present  axes  (i.e.  in  the  notation  of  the 
footnote  p  =  qpcf1).     Thus  our  equation  of  motion  is  now 

d-  +  Va>a  =  P  -  P'Vp/D, 
but  now 

d/dt  =  d/dt  -  Sp'pV  .  =  d/dt  -S(a-  cop')  P-V, 

whence  changing  p  into  p  we  have 

d<r/dt-S((T-a,p)V  .a+  V(oa  =  T-Vp/D  (14), 

which  is  the  equation  Greenhill  gives  on  p.  446. 

*  This  is  quoted  as  a  known  result  because  it  occurs  generally  in  the  subject  of 
Rigid  Dynamics.  No  Quaternion  proof,  however  as  far  as  I  am  aware,  has  been 
given.  We  therefore  give  one  here.  What  is  meant  by  rotating  axes  may  be  thus 
explained.— Instead  of  choosing  as  our  coordinates  the  vectors  a,  /3' ...  which  occur 
in  any  problem,  we  take  others  a,  /3...  such  that  a'  =  qaq~\  §'  =  q^q~1 ....  q(  )q~* 
may  be  called  the  integral  rotation  of  the  axes.  Thus  if  we  say  that  the  vector 
angular  velocity  of  these  axes  themselves  is  a  we  mean  that  the  real  angular 
velocity  is  quq-1,  so  that,  as  can  be  seen  by  putting  in  equation  (a)  below  a = const., 
or  as  in  equation  (58)  §  33  above,  a=2Vq~1q.  (This  maybe  established  also  by 
Tait's  Quaternions,  §  356,  equation  (2),  from  which  gwgr1=2F£g-1  or  a=2Vq~1q). 
Again,  when  we  say  that  the  rate  of  increase  of  a  in  space  is  r,  we  mean  that 
d'  =  (/T</-1  or  r  =  q~1  .  d  (qaq-^/dt .  q,  or 

T  =  d+2VVq-lq.a  («), 

or  r  =  d+  Vwa (^)- 

This  could  have  been  proved  with  fewer  symbols  and  more  explanation,  but  the 
above  seems  to  me  the  most  characteristic  Quaternion  proof.  We  might  have 
started  with  not  quite  so  general  an  explanation  of  reference  to  rotating  axes  and 
so  refrained  from  introducing  the  integral  rotation,  and  therefore  also  q. 
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The  Lagrangian  Equations 


* 


68.     We  now  consider  the  history  of  a  single  particle  and  for 
this  we  require  different  notation. 

We  consider  the  vector  coordinate  (p)  of  a  particle  as  a  func- 
tion of  some  other  vector  a  (say  the  initial  value  of  p)  and  of  t. 

We  first  require  the  connection  between  aV  and  pV.     We  shall 
drop  the  affix  a  and  retain  p,  so  that  now  not  Q  (Vx,  p^  but  Q  (V1?  ax) 

Now  SdaV  =  8dPpS7, 

but  dp  =  —  pfidaV^ 

.-.     SdaV  =  -SdofflSPlpV, 

or  since  da  is  perfectly  arbitrary 

V  =  -VX^1PV    (15). 

Thus  operating  upon  equation  (9)  §  65  by  V^pj  (  )  and  remember- 
ing that  in  that  equation  V  must  be  changed  to  PV  we  get 

V18Pl&  =  V(v  +  P)    (16), 

and  this  is  our  new  equation  of  motion  (Greenhill's  (1)  (2)  (3) 
p.  448). 

In  equation  (1)  §  12  above  let  us  put  for  p  +  7],  p  ;  and  for  p,  a. 

Thus  Xc0  =  -Sa>V.p    (17). 

Hence  by  equation  (Qd)  §  14  above, 

strained  vol.  of  el.  /  original  vol.  =  SV 1V 2V '^p^jjJG, 
whence  we  see  that 

^,Wpm  =  6  A (18)> 

where  D0  is  a  constant  which  when  a  is  taken  as  the  initial  value 
of  p  is  the  original  density.     This  is  the  equation  of  continuity. 


* 


[Note  added,  1892.  At  the  time  of  writing  the  essay  I  did  not  notice  that 
these  equations  are  a  particular  case  of  the  general  equation  for  an  elastic  body 
already  established  (see  equations  (15  m)  §  16  and  (15  n)  §  16  a) 
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Gauchy's  integrals  of  these  equations. 

69.  To  obtain  these  integrals  we  require  to  express  pV  in 
terms  of  V.  Now  equation  (15)  of  last  section  expresses  V  as  a 
linear  function  of  pV.  The  converse  we  have  already  seen  how  to 
get.     In  fact  from  equation  (6A)  §  8a  above 

2JpV  =  -VPlp2SVy2V (19), 

where,  with  Greenhill,  for  brevity  J  is  put  for  SVy^^p^pJQ. 
The  only  function  we  wish  to  apply  this  to  is  VpVa.     We  have 

2JVpV<r=V<r3Vp1pj3V1V2Va 

or  interchanging  the  suffixes  1  and  2  in  the  last  term 

JVpVa  =  p2SPla3SV1V2V3 (20), 

which  gives  the  spin  at  any  instant  in  terms  of  our  present  inde- 
pendent variables  a  and  t. 

70.  In  order  to  obtain  Cauchy's  integral  of  equation  (16) 
operate  on  it  by  W  (  ).     Thus 

0  =  V(V1  +  V2)  VtSPla2  =  SPl*2 .  VV2VV 
Now 

d  (SPla2 .  VV^J/dt  =  Strs,  ■  VV2V1  +  SPl&2 .  FV2V, . 

The  first  term  of  this  last  expression  is  zero  since  the  sign  is 
changed  by  interchanging  the  suffixes.  From  the  last  expression 
then 

8pl<r2 .  W^t  =  constant 

=  initial  value  of  £fo-a7V8f  =  -  VVa0, 

where  a0  is  the  initial  value  of  a.  Changing  the  suffixes  and  sub- 
stituting in  equation  (20)  we  have 

JVpVa  =  -PlSV<T0Vlf 

or  giving  Jits  value  DJD  from  equation  (18)  and  putting  V^cr  =  e, 
Wa0  =  e0,  we  have 

e/D  =  -S(eJD0)V.p    (21). 

This  is  Greenhill's  equations  (4)  (5)  (6)  p.  448. 

The    physical    interpretation    of   this  equation   is  quite   easy. 
Consider  a  small  vector  da  drawn  in  the  fluid  initially.     At  the 
M.  0 
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time  t  this  will  have  become  dp  =  —  SdaV  .  p.  Thus  we  see  that  if 
a  small  vector  ceJD0  be  drawn  in  the  fluid  initially  it  will  at  the 
time  t  be  cefD,  from  which  we  infer  that  an  element  of  a  vortex 
filament  will  always  remain  an  element  of  a  vortex  filament ; 
or,  a  vortex  filament  or  tube  always  remains  a  vortex  filament 
or  tube.  Again  we  see  that  Te/D  at  any  time  varies  directly 
with  the  elongation  in  the  direction  of  the  vortex  filament  so 
that  Te  varies  as  that  elongation  x  the  density,  i.e.  inversely  as 
the  cross-section  of  a  small  vortex  tube  at  the  point.  This  is  not 
the  easiest  way  of  arriving  at  these  results,  but  it  is  well  to  show 
in  passing  how  easy  of  interpretation  are  our  results. 

We  see  from  equation  (21)  that  if  e0  =  0,  e  =  0.  In  other  words, 
if  the  motion  have  a  velocity  potential  at  one  instant  it  will  have 
one  always. 


Flow,  circulation,  vortex-motion. 

71.  We  are  about  to  consider  vortex-motion  from  another 
point  of  view,  viz.  that  of  circulation. 

In  §  12  we  saw  that  a  strain  due  to  a  small  displacement  n 
could  be  decomposed  into  a  pure  strain  followed  by  a  rotation,  the 
vector  rotation  being  Wij/2.  If  now  for  n  we  put  the  small 
vector  aSt  we  see  the  propriety  of  calling  FV<x/2  the  spin.  This 
therefore  is  taken  as  a  definition  of  spin.  Greenhill  does  not  take 
this  (usual)  course  but  uses  the  property  we  shall  immediately 
prove  concerning  circulation  and  spin  to  lead  to  his  definition. 

It  is  not  necessary  here  to  define  flow  and  circulation.  Putting 
a  =  V(£  we  see  that  for  irrotational  motion  the  flow  =  —  JSadp  =  Jdcp 
from  one  point  to  another  is  the  increment  in  the  velocity  potential. 
Thus  for  mutually  reconcilable  paths  it  is  always  the  same. 

Taking  the  circulation  round  a  closed  curve  in  the  general  case, 
the  curve  not  inclosing  any  singular  region  of  the  fluid,  we  may 
transform  the  line  integral  into  a  surface  integral  by  equation  (8) 
§  6  above.     Thus 

-JS<rdp  =  -ffSdZVa    (22), 

so  that  the  circulation  round  the  curve  equals  twice  the  surface 
integral  of  the  spin.     Hence  Greenhili's  definition  of  the  spin. 
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72.     From  equation  (9)  §  65  we  see  that 


01 


J  A 


d(v+P), 


v  +  P+^/2 


.(23). 


-dt{jA8trdP)  = 

Therefore  for  a  closed  curve 

d  (JSadp)/dt  =  0, 

so  that  the  circulation  round  the  curve,  and  therefore  the  cor- 
responding surface  integral  of  the  spin  remains  always  constant. 
Taking  the  curve  round  a  small  vortex  tube  we  once  more  arrive 
at  the  propositions  enunciated  in  §  70  about  vortices. 

Thus  at  any  point  of  a  vortex  tube  the  strength  which  is 
defined  as  the  product  of  the  spin  into  the  cross-section  is  constant 
throughout  all  time.  Also  it  is  the  same  for  all  points  of  the  tube, 
for  by  equation  (9)  §  6  we  have 

fJSd2V<r=JfJSV*od%  =  0 

for  any  portion  of  the  tube.  But  the  only  part  contributing  to 
the  surface  integral  is  the  ends  of  the  part  of  the  tube  considered, 
so  that  the  strength  at  these  two  ends  is  the  same,  and  is  there- 
fore constant  for  the  whole  tube. 

73.  These  propositions  are  proved  in  the  paper  by  Hicks 
already  referred  to.  There  is  yet  a  third  method  of  proof  which, 
like  Hicks's,  is  derived  directly  from  the  equation  (9)  of  motion. 

We  have  from  equation  (1)  §  62 


dt  dt 


therefore 


dVVa/dt  -  VV&  =  FCV^o-^.  <r2). 


.(24)* 


*  [Note  added,  1892.  When  first  giving  this  in  the  Mess,  of  Math.  1884  and 
when  again  putting  it  in  the  present  essay,  I  was  unaware  of  Prof.  Tait's  paper  in 
Proc.  R.  S.  E.  1809—70,  p.  143,  where  in  the  present  notation  he  has  for  an  incom- 
pressible fluid,  (1)  djdt  =  dldt-  Say  [given  in  Cartesian  form],  (2)  <r=  -yv-yp/D, 
(3)  Sv<r=0,  (4)  yff-d(v(r)ldt=:-SV(rv.<r,  (5)  Vyar=0,  (6)  d^trfdt=  -  SyV  •  «r- 
Perhaps  I  have  not  interpreted  Prof.  Tait's  notation  which  is  but  briefly  described 
correctly,  but  (4)  should  apparently  be  V\& -  d  (\'a)ldt  =  S^<r\'  .a,  and  so  agree 
with  (5)  and  (6).  It  will  be  seen  that  the  whole  of  §§  512 — 13  of  Tait's  Quaternions, 
3rd  ed.,  is   contained  in  this  essay.    I  cannot  at  present  recall  whether  this  is 

6—2 
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Now  by  equation  (9)  §  65,  W&  =  0,  so  that 

dVValdt=VVl(T2SV2<T1 (25), 

.-.   d  (p  +  ccW<r)/dt  =<r+  xWcr  +  ccVV l(r2SV .a^ 

Now   ^SV^V^Vct^SV^  +  FV,V2,SV.,+  W^JSV^, 
but  W1V2S(T1a2  =  0,  for  an  interchange  of  suffixes  leads  to  a  change 
of  sign.     Therefore 

VV^SV^  =  VVaSVa  -  SVo-V  .  a. 

Put  now  cc  =  c/D  where  c  is  some  small  constant.     Thus 

w/c  =  -D/D*  =  -SV<rlD, 

and  we  get      d (p  +  cWa ./D)dt  =  a-  cSV<tV  .a.(D (26). 

This  shews  that  p'  =  p  +  cVS/a/D  is  the  variable  vector  of  a 
particle,  for  the  equation  asserts  that  the  time-flux  of  p  =  the 
velocity  at  p.  Thus  again  we  get  the  laws  of  vortices  given  in 
§70. 

Irrotational  Motion. 

74.  We  use  this  heading  merely  to  connect  what  follows  with 
what  Greenhill  has  under  the  same  on  p.  449.  It  is  not  very 
appropriate  here.  . 

For  irrotational  motion  we  have  seen  that  we  may  put 

o-  =  V0    (27). 

If  the  fluid  be  incompressible  we  have  further 

0  =  £V<r=  V2<£ (28). 

Let  T  be  the  kinetic  energy  of  this  liquid.     Thus 
-  2T=  DffJ(V<f>yd%  =  Dff4>Sd2V<l>  -  Dfff<f>V2<f>d% 
by  equation  (9)  §  6  above.     Thus  by  equation  (28) 

T=-±Dff(f>SdlV(p (29). 

Hence  we  see  that  if  we  have  a  vector  a  which  satisfies  the 
equations  Wa  =  0  and  SVa  =  0  throughout  any  singly-connected 
space,  or  more  generally  satisfies  the  second  of  these  equations 

owing  wholly  or  in  part  to  my  being  indebted  to  other  old  papers  of  Prof.  Tait,  or 
whether  in  writing  his  3rd  ed.  he  has  arrived  independently  at  the  same  treatment, 
but  am  inclined  to  the  latter  belief.] 
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throughout  any  space,  and  also  has  zero  circulation  round  any 
closed  curve  in  that  space*;  and  further  satisfies  the  equation 
SdSa  =  0  at  the  boundary  ;  then  the  function 

fffo>d%  =  0, 

or  o-=  0  at  all  points  of  the  space,  for  each  element  of  the  integral 
being  essentially  negative  must  be  zero. 

75.  The  following  important  theorem  now  follows  :  If  there 
are  given  ;  at  every  point  of  any  region  the  convergence  $Ver  and 
the  spin  W<r,  at  every  point  of  the  boundary  the  normal  velocity 
(and  therefore  Sd^a),  and  the  circulation  round  every  cavity  which 
increases  the  cyclomatic  number  of  the  region  ;  then  the  motion 
(as  given  by  a  at  every  point)  if  possible,  is  unique.  For  let  a  be 
one  possible  velocity,  and  if  possible  <r+  r  another.  Thus  r 
satisfies  all  the  conditions  that  a  does  at  the  end  of  last  section, 
and  therefore  t  =  0  at  every  point  or  a  is  unique. 


Motion  of  a  solid  through  a  liquid. 

76.  We  assume  that  there  is  no  circulation  of  the  liquid  for 
any  cycle ;  in  other  words,  that  if  all  the  solids  be  brought  to  rest, 
so  will  be  also  the  liquid. 

We  shall  take  axes  of  reference  fixed  in  the  (one)  moving  solid. 
In  the  foot-note  to  §  67  is  explained  how  the  rotation  of  these  axes 
is  taken  account  of.  The  effect  of  the  translation  of  the  origin 
will  be  easily  found. 

Let  a,  a)  be  the  linear  and  angular  velocities  respectively  of 
the  moving  solid.     Let  us  put 

<P  =  -S<rf-Sa>x  (30), 

where  <f>  is  the  required  velocity  potential  of  the  liquid,  and  yjr  and 
^  are  two  vector  functions  of  the  position  of  a  point  independent 
of  a  and  a).  Let  us  see  whether  yjr  and  %  can  be  found  so  as  to 
satisfy  these  conditions,  a  and  co  are  of  course  assumed  as  quite 
arbitrary. 

The  conditions  are  first  the  equation  (28)  of  continuity 

V2<£  =  0  (28), 

*  To  insure  that  <p  is  single  valued  and  so,  that  the  part  of  the  surface  integral 
of  equation  (29)  due  to  "  barriers  "  is  zero. 
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which  gives 

V^  =  0,  V^  =  0  (31), 

and  second,  the  equality  of  the  normal  velocity  —  SUdXV<f>  of  the 
liquid  at  any  point  of  the  boundary  with  that  of  the  boundary  at 
the  same  point.     Thus  for  a  fixed  boundary 

So-  (SdXV  .f)  +  Sa>  (SdZV  .  x)  =  0, 

whence  on  account  of  the  arbitrariness  of  a-  and  &> 

SdZV.iJr  =  Sd2,V.x  =  0    (32). 

For  a  moving  boundary  again  we  have 

Sd2V<j>  =  Sd%(<r+Va>p), 

or  So-  (dX  +  SdXV  .f)  +  Sco(  VpdX  +  SdXV  .  %)  =  0, 

which  gives 

dZ+Sd&V  .^  =  0 (33), 

Vpd£+Sd2tf.x  =  0   (34). 

Now  it  is  well  known  that  yfr  and  %  can  be  determined  to  satisfy 
all  these  conditions.  In  fact  x  being  a  coordinate  of  either  tJt  or  % 
these  conditions  amount  to  : — V2*-  =  0  throughout  the  space,  and 
SdX^x  =  given  value  at  the  boundary. 

77.  We  do  not  propose  to  find  i|r  and  %  in  any  particular 
case.  We  leave  p.  454  of  Greenhill's  article  and  go  on  to  p.  455, 
i.e.  we  proceed  to  find  the  equations  of  motion  of  the  solid.  For 
this  purpose  we  require  the  kinetic  energy  of  the  system.  Calling 
this  T  we  shall  have 

T=  -  S<rXa/2-  Sg>®<t  -  SaClm/V (35), 

where  %,  CI  and  <X>  are  linear  vector  functions  and  the  first  two  are 
self-conjugate.  This  is  the  most  general*  quadratic  function  of 
w  and  a.  It  involves  21  independent  constants,  six  in  X,  six  in  CI 
and  nine  in  <I>.  When  -vjr  and  %  are  known  X,  <E>  and  CI  are  all 
known.  We  will  obtain  the  expressions  for  them,  and  for  simplicity 
will  assume  that  the  origin  is  at  the  centre  of  gravity  of  the 
moving  solid.  Let  M  be  the  mass  of  this  last  and  /j,a>  (where,  as  is 
well-known,  ja  is  a  self-conjugate  linear  vector  function)  the 
moment  of  momentum.     Thus  by  equation  (29) 

2T  =  -  DfJ<pSdXV<p  -  Ma1  -  SufM. 

*  For  let  A  (c,  a) +13  (w,  <r)+  G  (w,  w)  be  this  general  function  A,  B,  and  C  being 
scalar  functions  each  linear  in  each  of  its  constituents,  and  let 
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Putting  SdSV<f>  =  SdS(<r+V<op)  at  the  moving  boundary  and 
zero  at  the  fixed, 

2T  =  DJKSaf  +  Scox)  (S<rdX  +  ScopdX)  -  Ma2  -  Sofia 

where  the  surface  integral  must  be  taken  only  over  the  moving 
boundary.     Thus  noting  that  2  and  12  are  self-conjugate 

XX  =  MX  -  ±DJJ  (yjrSXdX  +  dtSXf) (36), 

OX  =  pX -  iDff(xS\pdX  +  VpdZSXtf (37), 

-  SX®X'  =  ±DJf(S\xSX'dX  +  SXpdXSX'yfr)    (38). 

These  surface  integrals  can  be  simplified,  for  in  each  of  these 
equations  the  first  surface  integral  equals  the  second.  In  equation 
(36)  in  order  to  prove  this  we  have  merely  to  put  for  dX,  —  $dSV.-^r 
(equation  (33))  when  we  shall  find  by  equation  (9)  §  6  above 
(since  we  may  now  suppose  the  integrations  to  extend  over  the 
ivhole  boundary),  that 

JffSxdX  =  -  JfJf.SVy^Sx^ds  =JJdXSXyfr. 
Similarly  for  equation  (37).     Again  in  equation  (38), 

JjSXxSX'dX  =  -JJjS^1ViSXXlSX,yfrid%  =  JjSXpdXSX'f. 
Thus  finally  for  X,  £1  and  3> 

XX  =  MX  -  DJjfSXdX  =  MX  -  DjjdZSXyJr  (39), 

£IX  =  pX-  DjJxSXpdX  =pX-  DJJVpdXSXx    . .  .(40), 
_  S\®X'  =  DJjSxxSX'dX      =  DfJSXpdXS\'^ (41 ). 

This  last  may  be  put  in  the  following  four  forms 

cE>'\  =  -  DffdXSXx  =  -  DJM  SXpdZ (42), 

4>X  =  -  DJJxSXdX  =  -  DjjVpdXSX^ (43). 

Thus  assuming  that  t/t  and  x  are  determined  we  have  found  T 
as  a  quadratic  function  of  a  and  to. 

78.  If  P,  G  be  the  linear  and  angular  impulses  of  the  system 
respectively  our  equations  of  motion  are 

P  +  VcoT  =  F, 
G+Fg>G-M, 

by  the  footnote  to  §  67  above.  Here  F  and  M  are  the  external 
force  and  couple  applied  to  the  body.  Now  at  the  instant  under 
consideration  P  =  ffVT,  G  =  WVT.  But  if  the  origin  had  been  at 
the  point  -  p,  P  would  still  be  ^JT  whereas  G  would  be  WVT+  VpP. 
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Thus  differentiation  with  regard  to  t  does  not  affect  the  form  of  P 
but  does  that  of  G.  In  fact  using  the  last  stated  values  of  P  and 
G  along  with  the  last  two  equations  and  eventually  putting 
p  =  a,  p  =  0  we  get 

d(rVT/dt  +  Vco(rVT=F  (44), 

d<aVT/dt+Va>tyT+  V<t„VT=M (45). 

Now  from  equation  (35)  we  see  that 

ffVT=2<7 +  <£'*> (46), 

a>V77=^0-  +  na>.- (47), 

whence  from  equations  (44)  and  (45) 

%&  +  &6>  +  Vco^a  +  Vco&co  =  F (48), 

®&  +  Ho)  +  VaZa  +  ( V<r$>'co  +  V<o<&cr)  +  FwHw  =  M. .  .(49). 

Thus  we  see  that  with  Quaternion  notation  even  the  general 
equations  of  motion  are  not  too  complicated  to  write  down  con- 
veniently. 

We   now   leave   Greenhill's   article   and    proceed    to    certain 
theorems   nut   contained   therein.     The    Cartesian   treatment   of 
these  subjects  will  be  found  in  Lamb's  treatise  on  The  Motion  of 
Fluids,  Chapters  VI.  and  IX.,  which  are  headed  Vortex  Motion  and 
Viscosity  respectively. 

The  velocity  in  terms  of  the  convergences  and  spins. 

79.  We  have  seen  in  §  75  that  if  the  spin  and  convergence 
are  given  for  each  point  of  a  bounded  fluid  and  the  normal 
velocity  at  each  point  of  the  boundary,  there  is  if  any,  but  one 
possible  motion.  We  shall  see  directly  that  this  unique  value 
always  exists. 

At  present  we  observe  that  we  cannot  tind  a  motion  giving  an 
assigned  spin  and  convergence  for  each  point  and  any  assigned 
velocity  for  each  point  of  the  boundary.  If  however  with  such 
data  a  motion  be  possible  we  can  find  the  velocity  at  any  point 
explicitly.     By  equation  (19)  §  10  above  we  have 

4^(7  =  -  Vjjp7uad% 
but  by  equation  (9)  §  G 

JJjVu<rd%  =fjud%*  -fffuVad* 
...  4<Tr<r=fj[V]Vud2a-Jff[V]VuV<Td% (50), 
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where  the  square  brackets  indicate  that  we  may  or  may  not 
retain  the  V  at  our  convenience.     This  equation  may  be  put 

4tt<7  =//  [V]  (Vudta  -  udXVa)  +fffuV'*d% (51). 

Both  of  these  equations  solve  the  question  now  proposed. 

If  the  fluid  be  considered  infinite  the  surface  integral  of 
equation  (50)  vanishes  if  cr  converges  to  zero  at  infinity  and  also 
that  of  equation  (51)  if  in  addition  the  spin  and  convergence 
converge  to  a  quantity  infinitely  small  compared  with  the  re- 
ciprocal of  the  distance  of  the  surface. 

80.  We  may  now  consider  the  case  when  Vcr  (spin  and  con- 
vergence) is  given  at  all  points  and  SdX<r  (normal  velocity)  at  the 
boundary.  It  must  be  observed  that  the  given  value  of  the  spin 
must  be  distributed  in  a  solenoidal  manner  for 

SVVVcr  =  0. 

Whenever  the  quaternion 

4tt2  =  -JffVuVadZ  =  VJjjuVadS  (52), 

is  a  vector,  we  see  that  all  the  conditions  except  the  boundary 
ones  are  satisfied  by  putting  <r  =  q.  Let  us  then  see  when  q 
reduces  to  a  vector.     By  equation  (9)  §  6  above 

47nSf?  =  -  JjuSdXVa  +  fffuSVod*. 

Thus  q  is  a  vector  if  this  surface  integral  vanish.  The  surface 
integral  vanishes  if  all  the  vortices  form  closed  curves  within  the 
space.  If  they  do  not  we  must  extend  the  space  and  make  them 
form  closed  curves  outside  the  original  space.  Extending  the 
volume  integral  accordingly,  we  may  put 

47ro-/=-jy/Vi£V0-dS (53), 

and  now  Vcr'  =  Vcr.     Suppose  then  that 

cr  =  cr'+cr// (54). 

We  thus  get  Vcr"  =  0  and  may  therefore  put 

cr"  =  V0 (56), 

where  c/>  satisfies  the  equations 

V2<£  =  0 (57), 

and  SdlVcp  =  SdZ  (<r  -  a)  =  known  quantity    (58). 

Now  it  is  well  known  that  <f>  can  be  determined  so  as  to  satisfy 
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these  two  equations.  Therefore  the  problem  under  discussion 
always  admits  of  solution. 

The  above  is  equivalent  to  Lamb's  §§  128 — 131.     His  §  130  is 
the  natural  interpretation  of  the  equation 

4>7ra  =  -JJJWuVa-d%. 

His  §  132  is  seen  at  once  from  equation  (50)  above.  For  in  the 
case  he  considers  we,  in  accordance  with  §  7  above,  take  each  side 
of  the  surface  of  discontinuity  as  a  part  of  the  boundary.  Now 
[Sd%(r]a+b  —  0  so  that  for  this  part  of  the  boundary  we  can  leave 
out  the  part  Sd^cr  and  we  get 

4,77*  =  fjVVuVdZa-fJJVS7aVad% (59), 

so  that  if  we  regard  [VdXcr]a+b  as  —  2  x  an  element  of  a  vortex, 
we  get  the  same  law  for  these  vortex  sheets  as  for  the  vortices  in 
the  rest  of  the  fluid. 

81.  The  velocity  potential  due  to  a  single  vortex  filament  of 
strength  dd  is  at  once  obtained  by  putting  in  equation  (50)  for 
Wads,  IdOdp.  Thus  calling  a  the  part  of  the  velocity  due  to 
the  filament 

27TO-'  =  -  ddfWudp  =  -  dej[u1VVlVdSV1 

by  equation  (8)  §  G  above.     Now 

FV17aSV1  =  -  dtV?  +  Vl8dS,Vl. 

Hence,  since  V2?t  =  0  for  all  points  not  on  the  filament 

2TrcT,  =  V(ddjjSdZVu)  (60). 

Thus  the  velocity  potential  =  (27r)_1  x  the  strength  x  the  solid 
angle  subtended  by  the  filament  at  the  point  considered. 

Kinetic  Energy. 

82.  Let  us  put 

4m-q=JJJuVad» (61), 

so  that  tr  =  Vq  (62). 

We  may  now  find  expressions  for  T,  the  kinetic  energy,  in  one 
or  two  interesting  forms.     We  have 

2T=-JJfDa2d%, 
whence     2T  =  -  jjfDSaVqdS  =  -JJDSadSq  +JJJJ)l8<r1V1qd%. 
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We  assume  that  the  fluid  extends  to  infinity  and  that  the 
vortices  and  convergences  are  all  at  a  finite  distance,  so  that  at 
infinity  a  is  of  order  1/R2  and  q  of  order  1/R.  Thus  the  surface 
integral  vanishes  and  we  have 

ZT-fffB^^qd* (63), 

or,  putting  in  the  value  of  q  from  equation  (01), 

T  =  (Stt)-1  IfIfIIuD1Sa-1V1Vi<r2d%1d%2 (64). 

These  are  Lamb's  equations  (28)  and  (29),  p.  100,  generalised. 
Similarly  his  equation  (30)  generalised  is 

2T  =  fffDJSfHrlVl<r&   (65), 

for  JffD.Spa^^dS  =ffDSpadS<r  -JfJDSfcfcdz, 

by  equation  (9)  §  0  above.     But 

and  the  surface  integral  vanishes  as  before. 


Viscosity. 

83.  To  consider  viscosity  the  assumption  is  made  that  the 
shearing  stress  which  causes  it  is  p  x  the  rate  of  shear  of  the 
moving  fluid,  p  is  assumed  independent  of  the  velocity  and 
experiment  seems  to  shew  that  it  is  also  independent  of  the 
density.  This  last  we  assume  though  the  variation  with  density 
can  be  easily  treated. 

Consider  a  general  strain  %co.     Since 

mSXfiv  =  Sx^XllXv  =  SxX  vXfMXv  *> 

vxiixv  =  mx~ivtxv- 

Putting  p  and  v  for  any  two  vectors  perpendicular  to  w  we  see 

that  the  normal  Uco  to  any  interface  becomes  £7%'_1&)  by  the  strain. 

Hence  the  interface  a>  experiences  a  shear  (strain)  which  equals 

the  resolved  part  of  %£7a>  perpendicular  to  the  vector  %'"'&>  equals 

resolved  part  of  (%  —  ^~1)  Uco  perpendicular  to  %'_1&>.     Now  when 

^  is  the  strain  function  due  to  a  small  displacement  adt,  by  §  12 

above 

^co  =  (o  —  aS'coV  .  adt    (06)) 

*  Kelland  and  Tait's  Quaternions,  chap.  x.  equation  (n).     We  have  already  used 
a  particular  case  of  this  in  §  15,  above. 
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whence  x'^w  =  ro +  VlSa><r1dt     (67), 

and  we  see  that  the  shear  is  the  resolved  part  of 

-(SUcoV.cr  +  VfiUcoaJdt 

perpendicular  to  co,  i.e.  parallel  to  the  interface  w.  From  this  we 
see  by  our  assumption  concerning  viscosity  that  any  element  of 
the  fluid  is  subject  to  a  stress  <£  given  by 

</>&)  =  —  Ro)  —  fi  ($&>V  .  cr  +  Vj/SWj) 

R  being  a  scalar.     Now  we  define  the  pressure  p  by  putting 

3p  =  S&Z=3R  +  2fiSV<r 

.'.     <f>co  =  -p(o  +  IfmSVa  -  fx  (StoV  .  a  +  VfitoaJ. .  .(68). 

The  equation  of  motion  is 

Dfr  =  DF  +  </>A, 
or 

D  (dv/dt  -  SaV  .a)  =  D&  =  DF-Vp-  fiVSVa/S  -  /*vV. .  .(69). 

If  fx  be  not  as  stated  independent  of  D  this  equation  must  be 
made  to  contain  certain  space  fluxes  of  /a  which  are  quite  easy  to 
insert. 

84.  In  §  179  of  Lamb's  treatise  he  considers  the  dissipation 
function  due  to  the  viscosity  of  a  fluid  in  a  way  that  seems  to  me 
misleading  if  not  wrong.  It  appears  as  if  he  should  add  to  the 
first  expression  of  that  section 

udpxx/dx  +  vdpXyJdx  +  wdpxz[dx. 

He  refers  to  Stokes,  Camb.  Trans,  vol.  IX.  p.  58.  Let  us  give  the 
quaternion  treatment  of  Stokes's  method. 

If  T  is  the  kinetic  energy  of  any  portion  of  the  fluid 

2T=-JJJDa2d%. 

Thus  v     d(Dd%)/dt  =  0, 

T=-JJJDSa&d%, 

but  D&  =  DF  +  <f>1V1, 

so  that  t  =  -  JJJDSaFdS  -  JJJSvfaV \d% 

or  f  =  -  JjjDSaFdS  -JJSacpdZ  +///&r10V1ds (70). 

If  now  </>  be  expressed  in  terms  of  p  and  /u,  we  have  T  depend- 
ing on  F,  p  and  p.  The  part  of  T  depending  on  F  and  p  repre- 
sents energy  stored  up  as  potential  energy  of  position  and  potential 
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energy  of  strain  respectively,  but  the  part  depending  on  /x  repre- 
sents a  loss  of  energy  to  the  system  we  are  considering  the  energy 
being  converted  into  heat. 

Thus  putting  cf>=  p  +  rs  we  have  by  equation  (68) 

otw  =  §yu,a>$v7a  —  /jb(Sa)V  .  a  +  V^'coo-,)  I  7 

=  % /jLcoStyZ  +  2}ifa)  J  *•     >' 

where  ty  is  given  by  equation  (75)  below. 

Thus  we  see  that  the  rate  of  loss  of  energy  is 

JJSazrdS  -  JJJSa^Vfo  (72). 

The  surface  integral  is  the  work  done  by  viscosity  against  the 
moving  fluid  at  the  boundary  and  the  volume  integral  is  con- 
sidered due  to  the  work  done  against  the  straining  of  the  fluid. 
Thus  we  put 

F=-S<rt*rV1 (73), 

and  call  F  the  "  dissipation  function." 
By  equation  (6)  §  3  above 

F=-SW*r£    (74), 

where  -ty  is  the  rate  of  pure  strain  of  the  fluid,  i.e. 

2i/r&)  =  -  ScoV  .  a-  -  V1S(oa1 (75). 

Again  by  equation  (18)  §  18  above  because  F  is  quadratic  in  yjr 

F=-Sf^a^m (76)> 

so  that  from  equation  (74)  we  have,  by  §  4  above, 

«  =  *a*72  on 

Substituting  for  ■sr  from  equation  (71)  in  equation  (74) 

F  =  -  %»  (S^y  -  2rtW (78), 

which  gives  F  in  terms  of  yfr. 


SECTION    VI. 
The  Vortex- Atom  Theory. 

85.  If  Quaternions  can  give  valuable  hints  or  indicate  a 
promising  method  of  dealing  with  the  highly  interesting  mathe- 
matical theory  of  Vortex-Atoms,  I  think  this  alone  ought  to  be 
sufficient  defence  of  its  claims  to  be  within  the  range  of  practical 
methods  of  investigation. 

In  what  follows  I  think  I  may  be  said  to  have  indicated  a 
hopeful  path  to  follow  in  order  to  test  to  some  extent  the  sound- 
ness of  this  theory. 

Statement  of  Sir  Wm.  Thomsons  and  Prof.  HicJcs's  theories. 

86.  Sir  Wm.  Thomson's  theory  is  so  well-known  that  it  is  not 
necessary  to  state  it  in  detail.  Matter  is  some  differentiation  of 
space  which  can  vary  its  position  carrying  with  it  so  to  speak 
certain  phenomena,  some  of  which  admit  of  definite  quantitative 
measurement.  Perhaps  the  most  important  of  these  phenomena 
is  what  is  called  mass.  Now,  says  in  effect  Sir  Wm.  Thomson,  if 
we  suppose  all  space  filled  with  an  incompressible  perfect  fluid  the 
vortices  in  it  are  just  such  differentiations  of  space.  They  also 
carry  about  with  them  definite  quantitative  phenomena.  Can  we 
prove  that  these  hypothetical  vortices  would  act  upon  one  another 
as  do  the  atoms  of  matter,  the  laws  of  whose  action  are  contained 
in  the  various  Sciences,  e.g.  Physics,  Chemistry  and  Physiology  ? 
The  problem  in  its  first  stages  at  any  rate  is  a  mathematical  one, 
but  during  the  many  years  it  has  been  before  the  mathematical 
world  very  little  progress  has  been  made  with  it. 

Hicks's  extension  of  this  theory  is  perhaps  not  so  well  known, 
but  it  seems  to  me  quite  as  interesting  and  more  likely  to  tally 
with  the  known  phenomena  of  matter.     He  enunciates  his  theory 
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in  the  Proceedings  of  the  Cambridge  Philosophical  Society,  vol.  III. 
p.  27  G.  It  differs  from  Thomson's  simply  in  assuming  that  the 
fluid  does  not  quite  fill  space — that  there  are  in  it  bubbles*. 
These  bubbles  will  find  their  way  to  where  the  pressure  is  least, 
i.e.  speaking  generally  to  the  centre  of  some  at  least  of  the  vor- 
tices. Thus  we  have  another  source  of  differentiation  of  space 
and  general  considerations  seem  to  point  to  these  differentiations 
being  the  true  atoms,  though  of  course  "  atom "  is  no  longer  a 
descriptive  term. 


General  considerations  concerning  these  theories. 

87.  In  Maxwell's  article  Atom  in  the  Encyc.  Brit.  p.  45,  he 
says — "  One  of  the  first  if  not  the  very  first  desideratum  in  a 
complete   theory  of  matter  is  to   explain  first  mass  and  second 

gravitation In  Thomson's  theory  the  mass  of  bodies  requires 

explanation.  We  have  to  explain  the  inertia  of  what  is  only  a 
mode  of  motion  and  inertia  is  a  property  of  matter,  not  of  modes 
of  motion.  It  is  true  that  a  vortex  ring  at  any  given  instant  has 
a  definite  momentum  and  a  definite  energy,  but  to  shew  that 
bodies  built  up  of  vortex  rings  would  have  such  momentum  and 
energy  as  we  know  them  to  have  is  in  the  present  state  of  the 
theory  a  very  difficult  task. 

"  It  may  seem  hard  to  say  of  an  infant  theory  that  it  is  bound 
to  explain  gravitation." 

Now  as  Hicks  tells  us  what  induced  him  to  give  his  theory 
was  the  promise  it  gave  of  explaining  gravitation.  But  I  believe 
nowhere  does  he  point  out  the  still  more  important  result  that 
probably  on  his  theory  we  can  explain  inertia. 

The  statement  of  the  principal  property  of  inertia  put  scien- 
tifically is  that  the  motion  of  the  centre  of  gravity  of  any  two 
bodies  approximates  more  and  more  nearly  to  uniform  velocity  in 
a  straight  line,  the  more  nearly  they  are  isolated  from  external 
influence.  But  this  property  is  probably  true  of  Hicks's  bubbles. 
The  centre  of  gravity  of  any  portion  of  the  fluid  containing  certain 
bubbles  (1),  (2)  ...  (n)  will  if  approximately  isolated  from  all  the 

*  By  "tmbbles"  of  course  I  do  not  mean  spaces  occupied  by  another  kind  of 
fluid  of  smaller  density  but  actual  vacua.  Thus  a  bubble  may  start  into  existence 
where  none  previously  existed,  or  again  a  bubble  may  completely  disappear. 
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rest  of  the  fluid  move  approximately  in  a  straight  line,  but  this 
amounts  to  saying  that  the  centre  of  volume  of  the  n  bubbles  will 
also  move  uniformly  in  a  straight  line  if  the  centre  of  the  whole 
volume  (bubbles  and  fluid)  considered,  move  similarly.  These 
conditions  are  not  evidently  true,  but  I  think  they  are  probably 
so  when  we  consider  groups  of  large  numbers  of  bubbles. 

Whether  this  theory  explains  gravitation  is  one  of  the  principal 
questions  to  be  considered  in  the  following  first  trial. 


Description  of  the  method  here  adopted. 

88.  Finding  it  practically  impossible  to  consider  the  real 
problem  of  a  number  of  bubbles  in  the  liquid,  I  consider  the  fluid 
continuous  and  not  incompressible.     Let  us  assume  that 

D  =  tanh(p/c)  (1), 

where  c  is  some  very  small  constant  which  in  the  limit  =  0.     For 

ordinary   values   of  p,  D  very  nearly  =  1.     It   is   only   when  p 

becomes  comparable  with  c  that  D  varies.     When  p  =  0,  D  =  0. 

Also 

P=J(dp/D)  =  c\ogsmh(p/c) (2). 

For  ordinary  values  of  p  then,  P  =p  but  when  p  becomes  compar- 
able with  c,  P  varies  and  when  p  =  0,  P  =  -  co  .  If  we  assume 
then  for  the  greater  part  of  the  fluid  that  p  is  large  compared 
with  the  small  quantity  c  we  see  that  the  fluid  will  have  almost 
exactly  the  same  properties  as  a  liquid  containing  bubbles.  We 
may  now  apply  the  equations  of  motion  of  §§  64,  65. 

Now  we  know  the  velocity  at  any  point  of  an  infinite  fluid  in 
terms  of  the  spins  and  convergences  at  all  points.  From  this  we 
may  deduce  the  acceleration  in  terms  of  the  spins,  convergences 
and  time-fluxes.  It  so  happens  that  by  the  equations  in  §  73  we 
can  get  rid  of  the  time-fluxes  of  the  spins.  This  greatly  simplifies 
the  further  discussion  of  the  problem.  The  equation  that  we  thus 
deduce  forms  the  starting  point  of  our  investigation,  for  the  phe- 
nomena of  gravitation,  electro-magnetism,  stress,  &c,  are  exhibited 
and  measured  by  the  acceleration  in  bodies  due  to  their  relative 
positions. 

The  equation  we  obtain  at  once  gives  a  generalisation  of  the 
integral,  equation  (12)  §  65  above. 

At  the  end  of  this  section  I  give  an  equation  which  is  rather 
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complicated    but    which    promises    to   enable    us    to   deal    more 
rigorously  with  our  problem  than  I  profess  to  have  done  below. 

We  proceed  to  the  investigation  just  indicated. 

Acceleration  in  terms  of  the  convergences,  their*  time-fluxes,  and 

the  spins. 

89.     Let  us  put  for  the  convergence  and  twice  the  spin  on  and 
t  respectively,  so  that 

SVa  =  m (3), 

FV<t  =  t  (4). 


We  have  seen  that  in  equation  (50)  §  79  we  may  neglect  the 
surface  integral.     Thus 

4tto-  =  VJJJu  (m  +  r)d% (5), 

.-.   4nrdtr/dt  =  VJJJud  (m  +  r)/dt .  d%. 


o\ 


By  §  73  f  =  r/SfVo-  -  StV  . 

but  by  equation  (1)  §  62 

dr/dt  =  f  +  SaV  .  r, 
.-.      dr/dt  =  SaV  .  t  +  rSVa  -  StV  .  a, 
or  because  SVr  =  0 

dr/dt  =  tSA<t  -  <tS±t  =  VVVar  (0) 

.-.     kirda/dt  =  VJJJu  (  VV  Var  +  dm/dt)dS (7). 

This  equation  is  not  in  a  convenient  form  for  our  purpose,  for 
t  and  m  may  be  and  most  probably  will  be  discontinuous,  so  that 
it  is  advisable  to  get  rid  of  their  derivatives.  Moreover  it  is 
advisable  to  allow  the  time-flux  of  m  to  appear  only  under  the 
form  d  (md%)/dt  because  in  case  of  pulsations  the  time  integral  of 
this  expression  will  be  zero,  whereas  we  can  predicate  no  such 
thing  of  (dm/dt)ds,  or  indeed  of  md%.  Let  us  first  then  consider 
the  first  term  in  equation  (7),  viz.  VJJJu  FV  VardS.  Using  equation 
(9)  §  6,  and  neglecting  the  surface  integral  at  infinity,  this 
becomes 

-  VfffWu  Yard*  =  V  WJJJu  VardS    [§  9  ] 

=  VJJJu  Yards  -  VSVJJJu  VardS 

=  4tt  Var  +  VJJJSarViidS, 
by  equation  (19)  §  10  above. 

m.  7 


98  THE  VORTEX-ATOM  THEORY.  [§  89. 

Remembering  now  (§  65  above),  that 

&=d<rldt-V<TT-V(a*)l2, 
we  see  by  equation  (7)  that 

4tt6-  =  -  2ttV  .  <r2  +  VJJJSarVudS  +  Vfffu  (dm/dt)  d%  . .  .(8), 
whence  by  equation  (9)  §  65  above 

P  +  v-  o-2/2  +  (4aryLJff(S<rrVu  +  udm/dt)d%  =  H  ...(9), 

where  H  is  a  function  of  the  time  only.  This  is  a  generalisation 
of  equation  (12)  §  65,  for  assuming  that  t  =  0  we  know  by  §  79 
that  d<f>/dt  =  (47r)_1///M  (dm/dt)  d%. 

90.  This  integral  equation  may  be  put  into  several  different 
forms  by  means  of  equation  (9)  §  6  above.  The  form  that  is 
useful  to  us  is  the  one  in  which  instead  of  dm/dt  we  have  d(md$)/dt 
as  we  have  already  seen.     Now  d  (d%)/dt  =  —  md%.     Therefore 

d  (md%)/dt  =  (m  -  m2)  dg 

=  (dm/dt  -  SaVm  -  m2)  d%. 
Substituting  from  this  for  dm/dt,  and  then  transforming  by 
equation  (9)  §  6  so  as  to  get  rid  of  the  space  variations  of  m  in- 
volved in  SaVm  we  get 

ffju  (dm/dt)  d%  =  JJfvd  (md%)/dt  +  ///  (urn2  -  mSaVu  -  urn')  d% 
=  fffud  (md%)/dt  -  JffmSaVudS. 

Thus  from  equations  (8)  and  (9) 

47ra=-27rV  .  a-+VJJfSVu(V(TT-m*)dS+Vfffud(md%)/dt ...(10), 

and 

P+v-a'/2^7r)-^JPxd%8Vu{VaT-m<T)+ud(md%)/dt)=H...(ll). 

Sir   Wm.  Thomsons  Theory. 

91.  We  are  now  in  a  position  to  examine  the  two  theories. 
We  first  take  Thomson's,  which  is  considerably  the  simpler,  and 
which  therefore  serves  as  an  introduction  to  the  other.  We  have 
then  m  =  0.     Thus  equations  (10)  and  (11)  become 

4tt6-  =  -  2ttV  .  o-2  +  VfffScnVwte  (12), 

p/D  -  <r8/2  +  {^7r)-\fJjSarVud%  =  H  (13), 

for  in  the  present  theories  we  may  put  v  =  0. 

We  shall  consider  the  two  terms  on  the  right  of  equation  (12) 
separately.  The  second  term  gives  then  an  apparent  force  per 
unit  mass  due  to  a  potential 

-  (^7rT\fffSarVvdS (14). 
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Comparing  this  with  equation  (33)  §  46  above  we  see  that  this 
potential  is  the  same  as  that  of  a  magnetic  system  given  by 

4ttI=  Var    (15). 

Now  this  magnetism  is  zero  where  t  is  zero.  In  other  words 
it  is  only  present  where  the  vortex-atoms  are,  and  therefore  it 
cannot  be  so  distributed  as  to  give  an  apparent  force  of  gravitation, 
for  taking  the  view  of  magnetic  matter  expressed  in  equation  (34) 
§  46  we  see  that  there  must,  in  every  complete  vortex-atom,  be 
a  sum  of  magnetic  matter  exactly  =  0.  Nor  again  is  it  likely  to 
explain  the  phenomena  of  permanent  magnets,  because,  assuming 
that  for  any  given  small  space  including  many  vortex-atoms 
fffVardS  is  not  zero,  the  apparent  force  produced  will  affect  all 
other  parts  of  space  independently  of  whether  this  same  integral 
for  them  is  not  or  is  zero.  But  to  explain  the  phenomena  of 
permanent  magnets  we  must  assume  that  the  effect  takes  place 
only  on  portions  of  space  where  there  is  positive  magnetic  matter. 
This  term  then  gives  us  no  phenomena  analogous  to  physical 
phenomena.  As  a  matter  of  fact  it  probably  has  no  visible  effect 
on  large  groups  of  vortices,  for  there  is  no  reason  to  suppose  that 
the  vector  Var  is  distributed  otherwise  than  at  random. 

92.  Let  us  now  consider  the  other  term  in  equation  (12),  and 
neglecting  the  term  already  considered,  put 

o-  =  _V.o-72 (16). 

The  phenomena  resulting  from  this  are  the  same  as  would 
follow  from  a  stress  in  a  medium,  the  stress  beiug  an  equal  tension 
in  all  directions  =  —  a2/2.  Now  comparing  equation  (62)  §  58  with 
equation  (4)  §  89,  we  see  that  a  depends  on  t/47t  in  exactly  the 
same  way  as  does  H  on  C.  The  question  then  arises — is  the  stress 
we  are  now  considering  equivalent  in  its  mechanical  effects  upon 
the  vortex-atoms  to  the  stress  given  by  equation  (68)  §  60  above, 
which  explains  the  mechanical  effect  of  one  current  on  another  ? 
We  saw  in  §  60  that  this  stress  is  a  tension  —  W/Stt  along  the 
vector  H  and  an  equal  pressure  in  all  directions  at  right  angles. 
The  effects  then  would  be  the  same  only  if  o-  be  at  right  angles  to 
the  surface  of  our  atom.  But  this  is  obviously  not  in  general  the 
case.  From  this  analogy  we  can  see  however  what  approximately 
will  happen  to  our  atom.  For  instead  of  a  being  at  right  angles 
to  the  surface  it  is  in  all  probability  very  nearly  tangential. 
Assuming  that  it  is  actually  tangential  we  see  that  at  the  surface 
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of  the  atom  we  have  a  tension  exactly  corresponding  to  the 
pressure  which  in  the  electric  analogue  will  be  exerted  on  this 
surface.  In  other  words,  the  atoms  will  act  on  each  other  very 
approximately  in  what  may  be  called  a  converse  way  to  the  small 
circuits  in  the  electric  analogue;  i.e.  where,  in  the  electric  analogue 
there  is  an  attraction,  in  the  hydrodynarnic  case  there  will  be  an 
apparent  repulsion,  and  vice  versa. 

Now  each  vortex-atom  forms  a  small  circuit  and  therefore  acts 
in  the  converse  way  to  a  small  magnet.  In  other  words,  each  atom 
acts  upon  each  other  atom  as  if  it  were  charged  with  attracting 
magnetic  matter.  Thus  we  see  that  if  we  could  suppose  certain 
extra  atomic  vortices  to  exist  and  to  be  disposed  throughout  space 
in  what  at  present  must  be  considered  quite  an  artificial  manner 
with  reference  to  the  atomic  vortices  we  could  rear  up  a  fabric 
which  would  explain  gravitation.  This  conception  however  is  of 
very  little  use  for  our  present  purpose. 

From  these  considerations  I  think  we  have  every  reason  to 
believe  that  Thomson's  theory  in  its  native  simplicity  does  not 
promise  to  lead  us  to  the  physical  phenomena  of  matter.  We 
pass  on  therefore  to  Hicks's. 


Prof.  Hicks's  Theory. 

93.  Hicks  in  his  theory,  as  I  understand  him,  assumes  that 
the  bubbles  always  remain  associated  with  the  same  particles  of 
the  fluid.  This  of  course  is  probably  not  the  case.  By  reason  of 
the  variation  of  the  pressure  with  the  time  it  is  probable  that 
evanescent  bubbles  start  into  existence  and  disappear  at  various 
parts  of  the  fluid.     This  requires  some  few  preliminary  remarks. 

The  particles  of  the  fluid  with  which  bubbles  are  permanently 
(i.e.  throughout  the  greater  part  of  each  small  but  not  infinitely 
small  interval  of  time)  associated  are  those  where  the  intensity  of 
spin  is  greatest.  If  the  intensity  of  spin  is  quite  various  at 
different  points  we  shall  thus  have  vortices  where  there  is 
generally  no  bubble,  extra  material  vortices  in  fact.  We  must 
suppose  these  distributed  quite  at  random  till  the  more  exact 
mathematical  treatment  of  our  problem  leads  us  to  suppose  other- 
wise. Now  evanescent  bubbles  will  occur  rather  in  these  vortices 
than  in  parts  of  the  fluid  where  is  no  vortex  at  all  (if  we  may 
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suppose  such  parts  to  exist),  and  of  course  they  will  occur  more 
readily  in  stronger  than  in  weaker  vortices.  At  the  present  stage 
of  the  theory  then  we  may  suppose  evanescent  bubbles  to  occur  in 
all  parts  of  the  fluid.  As  a  first  approximation  to  the  considera- 
tion of  the  effect  of  these  bubbles  we  may  assume  a  part  m!  of  ra 
to  be  continuously  distributed  through  space.     Putting 

m  =  m'  +  M  (17), 

we  must  suppose  M  to  be  present  only  at  the  material  bubbles 
where  it  is  probably  discontinuous,  whereas  m  is  continuous 
throughout  both  the  material  vortices  and  the  rest  of  the  fluid. 

Now  when  on  account  of  variation  of  pressure  m'  and  M  are 
affected — is  it  probable  that  in  the  neighbourhood  of  a  permanent 
bubble  m  and  M  are  of  the  same  or  opposite  signs  ?  To  answer 
this  question  observe  that  what  we  call  m  continuously  distributed 
is  really  a  series  of  discontinuous  values  of  m  scattered  at  random 
through  space  so  that  m'  is  probably  very  small.  A  decrement  of 
pressure  will  cause  an  increment  of  evanescent  bubbles,  i.e.  a 
decrement  in  in'.  An  increment  of  the  permanent  bubbles  will 
also  take  place  the  magnitude  of  which  by  what  has  been  said 
concerning  m'  will  not  be  by  any  means  accounted  for  by  the 
decrement  in  in.  There  will  therefore  also  be  a  decrement  in  M. 
Similarly  for  an  increment  in  the  pressure.  M  and  in  may  there- 
fore be  assumed  to  be  of  the  same  sign. 

After  noticing  that  m  is  continuous  and  therefore  that  there  is 
no  objection  to  introducing  its  space  variations  we  are  furnished 
with  all  the  materials  necessary  for  discussing  our  problem.  The 
equation  we  shall  use  is  (10)  of  §  90  above.  We  divide  its  discus- 
sion into  two  parts  as  follows. 


Consideration  of  all  the  terms  except  — V {a2)/ 2. 

94.  The  reasons  that  we  have  already  seen  in  §  91  for 
neglecting  VjfjScrTVudS  still  hold  good  so  we  put  this  aside. 
This  is  not  the  case  with  —VfffmSaVudS  but  we  can  neglect 
"Vffj'ud  {md%)Jdt  for  the  average  value  of  d  (md%)/dt  for  any  particle 
is  zero.  The  only  term  to  consider  then  is  —Vfffm8o-Vud%, 
Putting  as  in  equation  (17)  m  =  m  +  M  we  have 

-  VfffMSaVudS  -Vfffm' SaVud*. 
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The  first  term  of  this  can  be  neglected  for  the  same  reasons  as 
for  neglecting  that  containing  V<tt.  Applying  equation  (9)  §  6  to 
the  last  term  and  neglecting  the  surface  integral  as  usual  we  get 

V///«£V (ma)  d%  =  VJJJu  [m  (m  +  M)  +  tS<rVm'}d%. . .(18). 

The  last  term  can  probably  be  neglected  though  we  cannot  give 
such  good  reasons  as  for  the  other  terms  we  have  neglected.  At 
any  rate  in  places  not  near  permanent  bubbles  SaVm  is  as  likely 
to  be  positive  as  negative  and  vice  versa,  so  that  such  portions  of 
space  will  on  the  whole  produce  no  effect  on  the  permanent 
bubbles.  If  /S'uVm'  contributes  anything  for  parts  of  space  in  the 
neighbourhood  of  permanent  bubbles  we  must  be  content  at 
present  with  the  assumption  either  that  the  contribution  is  in 
general  positive  or  that  if  it  be  negative  it  is  not  sufficient  to 
cancel  the  effect  of  the  positive  term  m'M.  Remembering  that 
m  is  small  compared  with  M  we  are  left  with  the  positive  term 
m'M.  This  as  can  be  easily  seen  from  the  form  of  equation  (18) 
leads  to  an  apparent  law  of  gravitation  for  our  permanent 
bubbles. 

95.  The  gravitational  mass  which  we  must  on  this  supposition 
assign  to  each  permanent  bubble  varies  as  the  average  value  of 
m'M  for  that  bubble*.  Now  we  saw  in  §  87  that  the  probable 
measure  of  mass  of  a  permanent  bubble  was  proportional  to  its 
average  size.  Do  these  two  results  agree  ?  I  cannot  say,  but  even 
if  they  do  not  these  considerations  would  still  explain  the  motions 
of  the  solar  system,  but  if  the  sun  and  Jupiter  (say)  were  to 
collide  their  subsequent  motion  would  not  be  that  due  to  the 
collision  of  two  bodies  the  ratio  of  whose  masses  is  that  which  is 
accepted  as  the  ratio  of  the  sun's  and  Jupiter's.  As  a  matter 
of  fact  however  I  should  imagine  that  the  average  value  of  m'M 
for  a  permanent  bubble  is  proportional  to  its  average  volume  and 
this  simply  as  a  consequence  of  the  reasoning  in  §  87  above. 

A  conclusion  at  any  rate  to  be  drawn  from  the  above  is  that 
there  is  a  presumption  in  favour  of  Hicks's  theory  explaining 
gravitation. 

*  There  is  one  important  difference  to  be  noticed  between  tbis  and  Hicks's 
explanation  of  gravitation.  His  depends  on  the  synchronous  pulsations  of  distant 
vortices.  I  do  not  wish  to  imply  that  I  do  not  believe  in  the  existence  of  such 
synchronous  pulsations,  but  by  the  above  we  see  that  gravitation  can  probably 
be  explained  independently  of  them.  As  a  matter  of  fact  such  synchronous 
pulsations  probably  actually  occur  on  account  of  the  variation  of  II  with  the  time, 
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Consideration  of  the  term  —  V  .  a2/2. 

96.  In  considering  this  term  we  adopt  the  method  of  §  92 
and  consider  an  electric  analogue.  The  analogue  is  an  electro- 
magnetic field  for  which  in  the  notation  of  §§  46  to  60  above  at 
every  point. 

H  =  <7    (20). 

For  this  field  we  have  at  once 

4ttC=  FVH  =  t  (21). 

The  distribution  of  the  magnetism  in  the  field  is  somewhat  arbi- 
trary, but  in  the  notation  of  equations  (61)  and  (62)  §  82  it  will  be 
found  that  everything  is  satisfied  by  putting 

4nrI  =  -VSq (22). 

This  gives  as  it  should  #V  (H+  4"rrl)  =  0,  which  is  in  fact  the  only 
equation  it  is  necessary  to  satisfy.     We  have  further 

B  =  H  +  4ttI  =  a  -  VSq  =  VVq, 

so  that  A=  Vq  (23). 

Thus  all  the  important  vectors  in  the  analogue  are  determined. 
It  remains  to  compare  the  mechanical  effects  of  the  analogue  with 
the  term  -  V  .  a2/2. 

I,  it  must  be  observed,  is  not  confined  to  the  bubbles,  but  is 
distributed  throughout  space. 

97.  If  we  now  assume  that  bubbles  have  not  always  existed 
in  the  positions  which  we  call  permanent,  there  cannot  at  the 
surface  of  the  bubbles  be  any  circulation  round  them.  This  makes 
the  velocity  at  the  surface  almost  normal  to  it,  so  that  the  stress 
given  in  equation  (6Q)  §  60  reduces  to  a  tension  —  />'H  (2B  —  H)/87r 
over  the  surface,  i.e.  we  have  a  pressure 

-o-'^Tr  +  SaVVq/^Tr. 

Now  on  account  of  the  absence  of  circulation  VF</  is  very  small 
and  may  therefore  be  neglected.  Thus  we  get  a  pressure  —  a'^Sir, 
and  are  thus  led  once  more  to  a  "converse"  of  the  analogue. 
This  at  once*  leads  to  another  reason  for  the  law  of  gravitation  if 


*  [Note  added,  1892.     Because  the  density  of  attracting  magnetic  matter  of  the 
analogue  =  -&'y-HM7r=  -m/4ir.] 
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the  pulsations  are  synchronous.  This  we  have  already  seen  to  be 
probable. 

The  present  consideration  of  the  subject  is  merely  to  point  to 
a  method  of  investigating  the  theory  of  vortex-atoms.  I  therefore 
leave  the  subject  here,  not  attempting  to  force  the  phenomena  we 
have  been  considering  to  tally  with  the  known  phenomena  of 
electricity  and  magnetism.  Nevertheless  I  may  say  that  the 
prospect  of  discussing  these  things  by  means  of  the  present  subject 
can  scarcely  be  considered  as  distant  after  what  has  gone  before. 

To  sum  up,  this  first  application  of  the  method  leads  to  a  pre- 
sumption in  favour  of  Hicks's  theory  leading  to  an  explanation  of 
both  the  important  properties  of  matter — inertia  and  the  law  of 
gravitation — and  there  is  also  reason  from  it  to  hope  that  the 
phenomena  of  electro-magnetism  are  not  unlikely  to  receive  an 
explanation.  Thomson's  theory  on  the  other  hand  would  seem  to 
fail  in  the  first  two  at  any  rate  of  those  endeavours. 

98.  We  close  the  essay  with  the  fulfilment  of  the  promise 
made  towards  the  end  of  §  88.  In  that  section  it  will  be  remem- 
bered we  considered  a  hypothetical  fluid  for  which  D  =  tanh  (p/c), 
and  made  this  do  duty  for  a  liquid  containing  bubbles.  Strictly 
speaking  our  liquid  is  bounded  at  the  bubbles  and  therefore  as  a 
bounded  liquid  should  it  be  treated.  For  such  a  liquid  we  require 
an  equation  corresponding  to  equation  (10)  §  90,  and  if  possible 
equation  (11)  also.  This  last  I  have  been  unable  to  obtain,  and  I 
am  not  sure  that  to  solve  the  problem  explicitly  is  possible. 

Our  problem  only  deals  with  an  incompressible  fluid,  but  as 
the  removal  of  this  restriction  does  not  greatly  complicate  the 
work  we  will  consider  the  general  case  of  a  bounded  compressible 
fluid.     We  have 

4™  =  V*fffu*d%  =  -  VfffVnadZ 

=  v  ///m  (t  +  m)  d%-V  ffudXa, 
whence 

4nrd<r/dt  =  Vfffd  {u  (r  +  m)  d%}jdt  -  V  ffd(ndla)ldt. 

The  justification  of  using  d/dt  on  the  left  and  d/dt  under  the  inte- 
gral sign  will  appear  if  the  increment  (da/dt)  dt  in  a  at  a  given 
point  in  the  time  dt  be  considered.  It  will  be  observed  that  the 
meaning  here  to  be  attached  to  it  will  be  —  SaVu,  as  in  the 
differentiation  d/dt  with  regard  to  the  time  the  origin  of  n  is 
assumed  to  be  fixed. 
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I  shall  now  merely  indicate  the  method  of  procedure.  By  the 
method  exhibited  in  §  89  we  can  prove  that 

d  (urd%)/dt  =  -  WuVa-rdS  -  u<tStM% 
and  that 

d  (umd%)/dt  =  —  mSaVuds  +  ud  {md%)jdt. 

From  this  we  can  deduce  that 

4nrd<r/dt  =  V  VV  fJjuVa-rdS  +  V///  {-  mS<rVud%  +  ud  (md%)/dt} 
-  VffuaSrdX  -  Vffud  (d%<r)/dt  +  V  JfdlaSaVu, 
and  from  this  again  we  get 

&  =  Viu  +  Vtv'  +  Vv (24), 

where  iv  and  w  are  scalars  and  v  a  vector  given  by 

4ttw  =  -  2tto-2  +///  {d%SVu  (  Vo-t  -  m<r)  +  ud  (md%)/dt) . .  .(25), 
so  that  iv  is  in  fact  the  H  —  v  —  P  of  equation  (11)  §  90 

4ttw'  =  ff  (SdSaSaV u  -  udSdta/dt) (26 ), 

4™  =//(-  uaSrdS  +  VdlaSaVu  -  udVdXa/dt)   . .  .(27). 

This  last  equation  may  be  put  in  what  for  our  purposes  is  the 
more  convenient  form 

4™  =  -//  {uaSrd^  +  d  (uVd2<r)/dt}  (28). 

Again  it  may  be  put  in  a  form  free  from  d/dt ;  for  Vd^a  =  0 
because  the  surface  is  a  free  surface  and  d(dX)/dt=^ViSd^.a1*)  as 
can  easily  be  proved  by  considerations  similar  to  those  in  §  83. 
Thus 

4ttv  =//  (-  uaSrdt  +  VdtaSaVu  -  uVV  ^SdZtrJ (29). 

In  the  problem  we  have  to  discuss  m  =0,  so  that  w  gives  only 
terms  which  we  discussed  in  considering  Sir  Wm.  Thomson's 
theory.  Observing  that  if  in  iv  we  change  md%  into  —  Sdlcr  we 
get  for  the  part  of  w  containing  m,  w  ;  we  see  that  w'  only  gives 
terms  that  we  have  virtually  discussed  under  Hicks's  theory.  We 
have  however  entirely  neglected  v.  Are  we  justified  in  this  ?  In 
the  first  place  we  have  seen  that  if  the  bubbles  have  not  always 
been  associated  with  those  parts  of  the  fluid  with  which  they  now 
are  there  is  round  every  bubble  absolutely  no  circulation.     This 

*  [Note  added,  1892.     This  should  be  d  (d2)jdt=  -  mdS  +  Vi^Sc^,  and  there- 
fore equation  (29)  should  be 

iwv  =  J'j  [  VdZaSoyu  -  u  (<r5rd2  -  mVdZa  +  Vv-pSdZaJ  \ . 
This  does  not  affect  our  present  problem  because  m  —  0  in  our  case.] 

M.  8 
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shows  that  for  any  one  bubble  JJVdXcr  =  0*,  and  therefore  we  are 
justified  in  neglecting  the  last  term  in  equation  (28).  We  are 
probably  also  justified  in  neglecting  the  first  terra,  for  probably  t 
is  very  nearly  tangential  to  the  surface  and  therefore  Srd%  =  0. 

[Note  added,  1892.  The  whole  of  this  last  section  is  in  rather 
a  nebulous  stage,  and  since  writing  it  I  have  not  had  sufficient 
leisure  to  return  to  the  matter.  I  hesitated  whether  to  include  it 
in  the  present  issue.  But  since,  notwithstanding  the  absence  of 
any  reliable  results,  it  serves  very  well  to  illustrate  how  investiga- 
tions are  conducted  by  Quaternions,  I  have  thought  it  worth  pub- 
lishing. 

Should  anybody  feel  inclined  to  attempt  to  apply  the  method 
or  an  analogous  one  it  is  well  to  note  that  in  the  Phil.  Mag.  June, 
1892,  p.  490,  I  have  given  the  more  general  result  sought  in  this 
last  section.  As  the  result  is  .not  there  proved  I  give  one  proof 
which  seems  instructive.  It  exhibits  the  great  variety  of  suitable 
quaternion  methods  of  dealing  with  physical  questions.  It  fur- 
nishes incidentally  a  fourth  quaternion  proof  of  the  properties  of 
vortices.  It  also  illustrates  how  special  quaternion  methods 
developed  for  use  in  one  branch  of  Physics  at  once  prove  them- 
selves useful  in  other  branches. 

Adopting  the  notation  and  terminology  of  Phil.  Trans.  1892, 
p.  686,  §§  5 — 7,  let  a  and  t  be  taken  as  an  intensity  and  flux  respec- 
tively, t  still  being  =  Wcr  and  therefore  t'  =  WV.  Thus  a-' 
and  t  are  the  actual  velocity  and  double  spin  respectively  and 
the  equation  of  motion  is 

o-'  =  -  V  (v  +  P). 

Putting  g'  =  \~xa  and  operating  on  the  equation  by  ^' 

*  One  way  out  of  many  of  proving  this  is  as  follows.  Divide  the  bubble  up  into 
a  number  of  infinitely  near  sections  by  planes  perpendicular  to  the  unit  vector  a. 
For  any  one  section  $Sdp<r  =  0.  Consider  d^,  to  be  the  element  of  the  surface  cut 
off  by  the  following  four  planes,  (1)  the  plane  of  section  considered,  (2)  the  con- 
secutive plane  of  section,  (3)  the  two  planes  perpendicular  to  dp,  and  through  the 
extremities  of  the  element  dp.  Thus  if  x  be  the  distance  between  the  two  sections, 
dp  =  x~1Vad1,,  whence 

JJSadSo-=0, 

for  the  surface  of  the  bubble  between  the  two  sections.  But  adding,  we  may  sup- 
pose this  integral  to  extend  over  the  whole  bubble.  Thus  SajjVd"E(r=0  for  the 
whole  bubble ;  therefore  o  being  a  quite  arbitrary  unit  vector  we  have  for  the  whole 
bubble  jjVdZa=0. 
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Now 

X'd  (x'^/dt  .*  =  -  *y-V  =  -  % V  =  V1&r1V  =  V .  </8/2. 

.-.     <5-  +  V.<r/2/2  =  -V(u  +  P) (A), 

which  can  easily  be  deduced  from  or  utilised  to  prove  Lord 
Kelvin's  theorem  concerning  "  flow,"  equation  (23)  §  72  above. 

As  d/dt  is  commutative  with  V,  f  =  0  or  t  is  an  absolute  con- 
stant for  each  element  of  matter.  This  being  interpreted  at  once 
gives  the  well-known  properties  of  vortices  in  their  usual  form. 

If  in  the  equation  4-7T  (v  +  P)  =  S .  V2  fffu  ( v  +  P)  d%  we  carry 
one  V  across  the  integral  sign,  get  rid  of  its  differentiations  which 
affect  u  by  equation  (9)  §  6  above,  and  then  do  the  same  with  the 
other  V  we  get 

4tt  (v  +  P)  =JJ  [0  +  P)  SdXVu  -  uSdXV  (v  +  P)}  +fffu  V2  (v+P)  d% 

At  surfaces  of  discontinuity  in  <r,  v  and  P  will  both  be  continuous, 
so  that  instead  of  JJ(v+P)  Sd!V u  we  may  write  ffb(v+P)  SdXVu. 
In  the  last  equation  substitute  throughout  for  V(?j  +  P)  from 
equation  (A).     Thus 

4tt  (v  +  P)  =//6  (v  +  P)  SdZVu  +JfuSdX&-fffuSV&d% 

+  (JJuSdXV  .  0-/2/2  -/J)W  V2ds/2) 
=  ffb  (v  +  P)  SdZVu  +jfiiSd%&  -jjfuSV&d%  +JffSVuV  (o-'2)  ds/2, 
which  is  equation  (36)  of  the  Phil.  Mag.  paper. 

If  the  standard  position  and  present  position  of  matter  coincide 
it  is  quite  easy  to  prove  that 

dSdt'a/dt  =  Sdl  (&  +  VV^aJ 
d  (SVaW)ldt  =  SV  {&  +  W^aJ  d%. 

Substituting  for  Sd1&,  SV&  from  these  in  the  last  equation  we 
get  equation  (32)  of  the  Phil.  Mag.  paper;  but  this  equation  can 
also  be  proved  directly.  It  should  be  noticed  that  equations  (34) 
and  (35)  of  that  paper  have  been  wrongly  written  down  from 
equation  (32).     In  each  read  +  V  (o-2/2)  for  -  V  (o"72).] 
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